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'with  respect  to  the  torso.  Kinematics  of  relative  motion  between  two  body 

segments  is  presented  in  a  format  suitable  for  utilization  of  four  sonic  emitters 
on  the  moving  body  segment.  ]  The  necessary  computer  software  for  PDP  11/34A  is 
also  developed  to  collect  bc/th  kinematic  and  force  data  on  human  subjects. 

^The  numerical  results  are  presented  from  a  set  of  experiments  conducted  on 
three  male  and  three  female  subjects  to  determine  their  active  muscle  resistance 
against  the  external  forces  to  dislodge  their  limbs  from  the  initial  or  so-called 
I  stowed  position.  \  From  a  second  set  of  experiments,  extensive  numerical  results 
on  determination  of  the  active  muscle  force  capabilities  of  the  subjects,  when 
their  limbs  are  dislodged  from  the  initial  configuration,  are  also  presented  for 
various  orientations  of  the  limbs  with  respect  to  the  torso.  In  addition  to 
these  two  sets  of  experiments,  some  representative  results  of  a  special  set  of 
experiments  dealing  with  the  forced  kinematic  motion  of  the  shoulder  complex  are 
included  in  the  report. 

Vit  is  concluded  that  although  there  are  intra-  and  inter- subject  variations 
for  the  maximum  values  of  the  resistive  muscle  force  and  moment  data  for  the 
limbs,  there  are  some  trends  one  can  establish  for  the  behavior  of  their 
magnitudes .  ^It  is  expected  that  incorporation  of  the  results  of  the  present 
research  intov the  multi-segmented  mathematical  models  of  the  human  body  should 
improve  the  lopg-time  response  capabilities  of  these  models  so  that  they  can 
simulate  more  realistically  the  biodynamic  events  which  take  place  prior  to 
flail  injuries. i 
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INTRODUCTION 


Abandoning  disabled  aircraft  which  can  not  be  flown  may  be 
accomplished  by  either  bailing  out  via  seat  ejection,  or  by  the  ejection 
of  the  encapsulated  crew  station  containing  one  or  more  aircraft 
personnel.  If  the  separation  from  the  disabled  aircraft  occurs  by  seat 
ejection,  the  body  of  the  ejected  crew  member  is  subjected  to  a  high-speed 
air  stream  (the  windblast)  generating  large  magnitudes  of  aerodynamic 
forces  on  the  higher  and  lower  extremities  of  the  body.  Despite  the 
design  improvements  on  the  seats  with  provision  for  the  limbs  to  be 
lodged  in  a  favorable  position  during  ejection,  dislodging  of  the  limbs 
occur  frequently  when  the  body  is  ejected  into  the  windblast. 

Primary  cause  of  the  dislodging  of  the  limbs  is  the  difference 
between  the  magnitudes  of  the  deceleration  forces  acting  on  the  limbs  and 
the  torso.  Note  that  immediately  after  ejection,  the  seat  and  its  occu¬ 
pant  have  the  same  speed  as  the  aircraft;  however,  due  to  shape  (frontal 
area)  and  mass  differences  between  the  limbs  and  the  torso  to  which  the 
mass  of  the  ejection  seat  is  also  attached  the  aerodynamic  drag  forces 
acting  on  these  parts  will  be  different,  higher  on  the  limbs  than  on  the 
torso.  Thus,  there  is  a  strong  tendency  of  dislodging  of  the  limbs  from 
their  so-called  stowed  position.  The  dislodged  limbs  swing  freely  and 
without  control  and  they  are  stopped  only  by  some  part  of  the  seat  struc¬ 
ture  or  by  collision  with  another  part  of  the  body.  The  resulting,  so- 
called  "flail  injuries"  are  the  dislocations  of  the  major  articulating 
joints  and  various  types  of  long  bone  injuries.  The  probability  of  these 
injuries  increase  significantly  with  increasing  aircraft  speed  [1],  It 
was  also  estimated  that  [2]  almost  half  of  all  ejections  taking  place 
under  combat  conditions  result  in  flail  injury  or  death  to  the  pilot. 

Several  years  ago  recognizing  the  seriousness  of  the  flail  injury 
problem,  the  Aerospace  Medical  Research  Laboratory  (AMRL)  of  the  Air  Force 
Systems  Command  has  initiated  both  experimental  and  analytical  programs. 
Some  experimental  wind  tunnel  measurements  have  been  already  conducted 
[3,4,5]  and  a  theoretical  effort  has  been  started  by  Schneck  [6],  In 
Reference  [6]  a  potential  flow  solution  is  presented  for  estimating  the 
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pressure  distribution  around  the  lower  arm  when  the  seated  body  is 
subjected  to  windblast.  In  this  theoretical  study  by  inclusion  of  the 
effects  of  the  flow  separation  it  was  also  possible  to  predict  the 
presence  of  a  pressure  drag  force  which  acts  to  throw  the  lower  arm 
outward,  away  from  the  torso.  Results  of  this  study  can  provide  some 
input  data  in  regard  to  prescribing  time-dependent  forces  on  the  body  seg¬ 
ments  of  the  Articulated  Total  Body  Model  (ATBM)  of  the  AMRL.  The  ATBM 
is  particularly  a  useful  model  to  study  the  pilot  ejection  problem 
provided  appropriate  joint  torques  and  segmental  loading  data  to  simulate 
windblast  forces  are  known  and  incorporated  into  the  model. 

Approximately  five  years  ago,  a  research  program  under  co-sponsorship 
of  the  Mathematics  and  Analysis  Branch  of  the  AMRL  and  the  National 
Highway  Traffic  Safety  Administration  (NHTSA)  was  initiated  at  the 
Engineering  Mechanics  Department  of  The  Ohio  State  University  to  collect 
passive  resistive  force,  moment  and  torque  data  at  the  five  major  joints 
of  the  human  body.  The  research  was  conducted  with  some  obvious  limita¬ 
tions  on  live  human  subjects.  The  joints  which  were  studied,  in  order  of 
descending  priority,  were  the  shoulder,  knee,  hip,  elbow  and  ankle.  The 
final  report  [7]  of  this  forementioned  research  program  along  with  a 
series  of  papers  [8-15]  illustrate  various  aspects  of  the  research  and 
present  experimental  results  on  both  passive  and  active  resistive  force, 
moment  and  torque  data  at  the  major  articulating  joints  of  the  human  body. 
We  will  next  present  very  briefly  some  of  the  highlights  of  the  experi¬ 
mental  phase  of  this  previous  research  program  [7],  The  primary  purpose 
of  this  brief  background  presentation  is  to  illustrate  the  substantial 
changes  and  improvements  that  were  made  on  the  experimental  apparatus  in 
order  to  accomplish  tasks  of  the  present  research  program  sponsored  by 
the  Air  Force  Office  of  Scientific  Research. 

The  major  components  of  the  experimental  apparatus  developed  for  the 
previous  research  program  were  a  subject  restraint  system  (Fig.  1),  a 
global  force  applicator,  GFA  (Fig.  2)  and  the  exoskeletal  device,  ESD 
(Figs.  1  §  4) .  The  subject  restraint  system  with  its  pitch,  and  yaw 
capabilities  allowed  orientation  of  the  subject  to  various  selected 
positions.  This  was  an  important  feature  of  the  subject  restraint  system 


Fig.  1.  Overall  view  of  the  subject  restraint  system  and  data 
collection  equipment.  Subject  is  prepared  for  the  shoulder  force 
and  moment  data  collection,  [7]. 


Fig.  2.  Force  is  being  applied  by  means  of  the  global  force 
applicator  (GFA)  on  the  subject's  arm,  [7]. 
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Fig.  3.  Close-up  view  of  the  force  transducer  and  the  force  cuff,  [7] 


Fig.  4.  Close-up  view  of  the  exoskeletal  device  (ESD) ,  [7] 


which  assisted  in  the  mobilization  of  subject's  appropriate  body  segment 
at  a  constant  elevation  while  the  experimentor  moved  the  adjacent  body 
segment  by  means  of  the  GFA  throughout  its  entire  range  of  motion  as  shown 
for  the  shoulder  in  Figs.  2  §  3.  To  eliminate  the  gravitational  component 
of  the  moment  values  at  the  joint,  the  force  application  to  the  moving 
body  segment  was  made  in  a  horizontal  plane  while  the  elevation  of  the 
moving  body  segment  was  maintained  by  a  support  line  in  such  a  way  that 
in  the  direction  of  the  force  application  a  relaxed  floating  type  of 
motion  of  the  moving  body  segment  was  achieved.  Both  the  security  and  the 
support  characteristics  of  the  restraint  system  allowed  total  relaxation, 
thus  minimizing  active  muscle  effects  during  the  experiment. 

Forces  were  applied  to  the  moving  body  segments  of  the  subject  by 
means  of  the  GFA  (Fig.  2)  which  consists  of  four  links  and  eight  revolute 
joints  containing  high  precision  potentiometers.  The  outputs  of  the 
potentiometers  and  the  length  dimensions  of  the  links  were  utilized  to 
obtain  the  direction  as  well  as  the  point  of  application  of  the  force 
vector.  Note  that  in  Fig.  2  the  last  two  links  of  the  GFA  are  shown.  The 
GFA  is  terminated  by  a  force  transducer  and  a  force  cuff  which  is  free  to 
rotate  about  its  axis.  The  force  transducer  was  designed  and  built  to 
measure  all  three  components  of  the  force  and  moment  vectors.  Of  course, 
the  predominant  force  component  was  the  one  along  the  direction  of  the 
last  link  of  the  GFA  and  all  the  other  force  as  well  as  the  moment  compo¬ 
nents  were  relatively  small  in  magnitude  if  one  maintained  approximately 
perpendicular  force  application  on  the  moving  body  segment. 

The  third  major  component  of  the  experimental  apparatus  was  the  ESD 
as  shown  in  Fig.  4.  The  ESD  also  consists  of  eight  revolute  joints  whose 
rotations  are  monitored  by  means  of  high  precision  potentiometers. 

Although  there  were  major  design  differences  between  the  ESD  and  the  GFA, 
from  the  kinematics  point  of  view  the  ESD  can  be  considered  as  a 
miniaturized  version  of  the  GFA.  Fundamental  design  requirements  of  the 
ESD  were  (a)  capability  of  providing  complete  freedom  of  motion  between 
two  adjacent  body  segments,  and  (b)  capability  of  providing  sufficient  data 
to  determine  kinematics  of  the  relative  motion  between  two  body  segments 
which  are  connected  with  a  complex  anatomical  joint.  Principles  of 


orthopedic  bracing  were  utilized  in  attaching  the  ends  of  the  ESD  to  the 
moving  and  fixed  body  segments.  For  example,  for  the  shoulder  an  upper 
arm  shell  and  the  back  rail  of  the  chair,  where  the  torso  support  shell 
is  attached,  were  used  for  the  ESD  connections  as  shown  in  Figs.  1  §  2. 

For  the  fitting  of  the  upper  arm  shell,  special  care  was  given  to  locate 
the  bony  landmarks  such  as  the  medial  and  lateral  epicondyles  of  the 
humerus  so  that  contacts  can  be  maintained  between  these  landmarks  and  the 
shell.  In  addition,  the  relative  soft  tissue  motion  between  the  upper  arm 
and  the  shell  was  kept  minimal  because  of  the  following  reasons:  a)  light 
weight  and  very  low  friction  construction  of  the  ESD,  b)  balancing  of  the 
weight  of  the  ESD  by  means  of  pulleys  and  small  weights  in  such  a  way  that 
the  subjects  did  not  feel  its  presence,  and  c)  application  of  the  GFA 
force  on  a  portion  of  the  arm  away  from  the  shell  segment  to  which  the 
moving  end  of  the  ESD  was  attached. 

The  development  of  the  ESD,  the  GFA,  the  subject  restraint  system, 
and  the  associated  theoretical  concepts  were  made  to  achieve  at  least  three 
major  tasks.  The  first  one  was  the  quantitative  determination  of  the  range 
of  motion  in  the  joints;  the  second  one  was  quantitative  determination  of 
passive  resistive  force  and  moments  at  the  joints;  the  third  task  was 
quantitative  determination  of  passive  resistive  torques  associated  with  the 
rotational  motion  of  the  body  segments  about  their  long  bone  axes.  In  this 
research  program  three  male  subjects  were  tested  and  the  tests  were 
repeated  at  least  three  times. 

In  this  report,  first,  a  detailed  description  of  the  major  components 
of  the  experimental  apparatus  is  presented.  This  is  followed  by  a 
presentation  of  kinematics  of  relative  motion  between  two  body  segments 
which  is  cast  into  a  format  suitable  for  employment  of  four  sonic  emitters 
on  the  moving  body  segment.  Next,  a  description  of  the  computer  software 
for  the  collection  of  both  kinematic  and  force  data  on  human  subjects  is 
provided. 

The  numerical  results  are  presented  from  experiments  in  which  three 
male  and  three  female  subjects  were  tested.  Two  sets  of  experiments  were 
conducted  for  the  collection  of  active  muscle  force  and  moment  data.  The 
first  set  of  experiments  was  designed  to  determine  the  active  muscle 


resistance  of  subjects  when  there  is  a  strong  tendency  of  dislodging  of 
the  limbs  from  their  stowed  position.  The  second  set  of  experiments  was 
designed  to  determine  the  active  muscle  resistance  of  subjects  when  the 
limbs  are  dislodged  from  their  stowed  position.  In  this  set  of  experiments 
various  combinations  of  0  and  $  angles  for  the  moving  body  segments  were 
considered.  The  angles  0  and  <j>  define  the  orientation  of  the  moving  body 
segment  with  respect  to  the  fixed  body  segment.  In  addition  to  these  two 
sets  of  experiments,  numerical  results  of  a  set  of  experiments  dealing 
with  the  forced  kinematic  motion  of  the  shoulder  complex  are  presented. 

We  will  next  supply  a  description  of  the  major  components  of  the  experi¬ 
mental  setup  developed  for  the  tasks  of  the  present  research  program. 

DESCRIPTION  OF  THE  MAJOR  COMPONENTS  OF  THE  EXPERIMENTAL  SETUP 

Major  changes  were  made  in  almost  every  aspect  of  the  previous  program 
in  order  to  improve  the  accuracy,  repeatability  and  reliability  of  the 
data.  These  changes  are:  a)  design,  manufacturing  and  calibration  of  a 
new  force  and  moment  transducer,  b)  replacement  of  the  GFA  with  a  new 
force  applicator  which  utilizes  three  sonic  emitters,  c)  replacement  of 
the  ESD  with  an  elastic  cuff  which  holds  four  sonic  emitters,  d)  purchase 
of  sonic  digitizing  equipment  manufactured  by  Science  Accessories  Corp. 
(Model  6P6-3D-3) ,  e)  replacement  of  the  previous  obsolete  data  collection 
system  with  the  purchase  of  a  PDP  11/34  mini-computer  with  floating  point 
processor  manufactured  by  Digital  Equipment  Corp.,  and  f)  interfacing  a 
14-bit,  32-channel  analog-digital  converter  manufactured  by  Data 
Translation  Corp.  (Model  DT1712/5714-64)  with  PDP  11/34  computer. 

Further  details  about  the  items,  listed  above  and  developed  by  us  during 
the  research  program  will  next  be  presented. 

THE  FORCE  AND  MOMENT  TRANSDUCER  AND  ITS  CALIBRATION 

The  force  and  moment  transducer,  which  is  referred  to  as  the  force 
transducer  in  the  introduction  section  of  this  report,  is  physically  sim¬ 
ilar  to  the  force  transducer  utilized  in  previous  research  program  [7]  on 
the  measurement  of  resistive  torques  in  major  human  joints.  The  trans¬ 
ducer  which  was  designed,  built  and  calibrated  for  the  present  research 
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program,  however,  employs  a  number  of  significant  improvements  over  its 
predecessor.  These  improvements  are  basically  of  two  types:  1)  modifica¬ 
tions  to  the  physical  design,  and  2)  a  significantly  altered  calibration 
approach.  The  subsequent  discussion  will  describe  these  changes  in  some 
detail  and  indicate  the  benefits  derived  from  their  implementation. 

As  before  the  transducer  is  essentially  a  cylindrical  shell  with 
strain  gages  located  on  the  shell  surface  to  sense  force  and  moment  in 
three  orthogonal  directions  as  shown  in  Fig.  5.  However  certain  physical 
design  alterations  have  been  made  to  increase  the  transducer  sensitivity. 
Since  the  predominant  force  component  which  occur  in  use  is  in  the  axial 
or  z-direction,  design  changes  were  made  specifically  to  increase  the 
load  sensing  capability  in  the  axial  force  measurement  circuit.  This  was 
accomplished  in  two  ways:  1)  by  decreasing  the  shell  wall  thickness 
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Fig.  5  Schematic  drawing  of  the  force  and  moment  transducer 
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(from  .100"  in  the  previous  design  to  .032"  in  the  current  transducer), 
and  2)  by  increasing  the  number  of  T  rosettes  used  to  sense  axial  strain 
from  two  to  four.  The  choice  of  shell  wall  thickness  was  based,  primarily, 
on  a  consideration  of  the  tranducer's  ability  to  withstand  predicted  load 
configurations  while  maintaining  an  adequate  safety  factor  against 
buckling.  Of  course,  cost  and  machining  capability  were  also  considered. 
The  addition  of  two  more  T  rosettes  naturally  increases  the  sensitivity  by 
increasing  the  change  in  resistivity  and  hence  the  change  in  bridge 
voltage  for  a  given  applied  load.  In  this  way,  an  approximate  150-o 
increase  in  axial  sensitivity  over  the  previous  transducer  was  achieved. 

The  change  in  wall  thickness  also  served  to  increase  the  sensitivity  of 
all  the  other  load  measuring  circuits. 

The  force  transducer  currently  in  use  is  also  greater  in  length  than 
the  previous  one.  This  was  done  primarily  to  accommodate  more  gages  (12 
single  gages  and  6  rosettes  as  opposed  to  8  single  gages  and  4  rosettes 
previously)  and  also  to  improve  the  shear  force  measurement  capabilities. 
The  shear  load  measuring  circuits  have  been  changed  from  half  bridges  to 
full  bridges,  and  their  calibration  techniques  have  been  significantly 
altered.  The  calibration  techniques  will  now  be  discussed  in  some  detail. 

One  of  the  basic  problems  encountered  in  transducers  which  measure 
force  and  moment  in  multiple  directions  is  cross-talk  among  the  measuring 
circuits.  The  cross-talk  is  due  to  erroneous  output  of  the  circuits 
belonging  to  the  other  directions  when  a  force  or  moment  is  applied  in  a 
particular  direction.  This  effect  can  result  from  slight  deviations  in 
gage  alignments,  imperfections  occurring  in  the  gage  application,  or 
imperfections  which  occur  due  to  machining  tolerances  which  make  ideal 
load  applications  impossible.  While  reasonable  care  in  the  construction 
of  the  transducer  has  kept  the  cross-talk  to  a  minimum  we  also  account  for 
their  effects  in  the  subsequent  load  analysis. 

Calibration  was  performed  by  applying  known  loads  to  each  channel  in 
turn  and  recording  the  six  outputs  of  the  transducer.  Ideally,  only  the 
channel  corresponding  with  the  applied  load  should  register  the  load. 

Here  a  channel  designates  quantities  observed  on  the  strain  gage  bridges 
corresponding  to  Fx,  F^,  Fz,  M^,  M^,  and  Considering  primary  and 


cross-talk  signals  of  the  transducer  behave  linearly  with  applied  forces 
and  moments  throughout  the  working  range  of  the  transducer,  it  is  then 
possible  to  characterize  the  behavior  of  the  transducer  by  the  relationship 

[S]  =  [Cl  [L]  (1) 


where  the  matrices  [S]  and  [L]  are  the  column  vectors  representing  the 
bridge  output  voltages  (signals)  and  the  input  loads,  respectively.  The 
elements  of  matrix  [C]  are  determined  from  the  calibration  data.  The 
diagonal  elements  represent  the  primary  signals  and  the  non-diagonal 
elements  represent  the  cross-talk  effects. 

For  example,  consider  that  we  are  calibrating  the  axial  mode  of  the 
transducer  and  that  the  axial  bridge  is  our  first  measuring  circuit. 

Then  the  signals,  s^,  that  we  observe  as  the  bridge  outputs  are  as  follows 


If 


=  CnL1 

=  c21lx ' 


=  C31L1 

=  C4iL! 

=  C51L1 
=  C61L1 


(the  primary  signal) 


cross-talk 

(ideally,  these  should  be  zero) 


we  perform  calibrations  of  the  subsequent  modes  we  obtain: 


(2) 


mode  2 

mode  3 

S1  =  C12L2 

S1  =  C13L3 

(primary  signal)  s2  =  C22L2 

S2  =  C23L3 

etc 

s3  =  C32L2 

(primary  signal)  s3  =  C-jjLj 

S6  =  C62L2  s6  =  C63L3 


Note  here  that  the  C^'s  are  the  slopes  of  best  fit  lines  of  voltages 
obtained  at  a  variety  of  loads.  Therefore  for  any  aribitrary  loading,  the 
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total  observed  signals  for  the  transducer  channels,  S^,  arise  as  a  result 
of : 


S1  =  C11L1  +  C12L2  +  C13L3  +  C14L4  +  C15L5  +  C16L6 


S6  =  C61L1  +  C62L2  +  C63L3  +  C64L4  +  C65L5  +  C66L6 


(3) 


or. 


[S]  =  [C]  [L] 


The  coefficient  matrix  obtained  as  a  result  of  the  transducer  calibration 
is  shown  in  Table  1. 

With  these  mathematical  concepts  of  the  load  analysis  in  mind,  let  us 
consider  some  of  the  subtle  physical  aspects  of  the  transducer  calibration 
that  were  employed.  The  calibration  of  the  axial  force  (F„)  and  torque 


Table  1 

Calibration  Crosstalk  Matrix  [C] 
(Units  are  Volts/lb  or  Volts /lb-in  ) 


0.01532 

0.00023 

0.00000 

-0.00016 

0.00090 

0.00108 

-0.00090 

0.03986 

0.00000 

0.00097 

-0.00020 

-0.00390 

0.00037 

0.00012 

0.02987 

-0.00009 

0.00229 

0.00028 

0.00045 

-0.00003 

-0.00028 

0.02924 

0.00033 

0.00707 

0.00047 

-0.00011 

0.00094 

-0.00030 

-0.10616 

0.00164 

-0.00001 

0.00008 

-0.00021 

-0.00149 

-0.00106 

0.10876 

(M  )  modes  are  relatively  simple  due  to  the  relative  ease  with  which  we 
can  apply  these  types  of  loads.  The  application  of  pure  shears  (F  and 
F  )  and  pure  moments  (M  and  M  ),  is  very  difficult,  however,  in  the 

y  x  y 

absence  of  highly  specialized  equipment.  These  modes  were  calibrated  in 
the  following  manner. 

The  transducer  was  fixed  at  one  end  so  as  to  perform  as  a  cantilevered 
beam  and  steadily  increasing  loads  were  applied  in  two  specific  locations 
as  shown  in  Fig.  6.  Of  course  these  types  of  loading  involve  both  a  shear 
load,  =  F,  and  a  bending  moment  with  respect  to  the  moment  bridge 
location,  =  d^xF  and  =  d2xF.  Therefore,  we  have  two  sets  of  signals: 
As  a  result  of  Lj  =  F  and  L2  =  d^xF: 

S1  =  C11L1  +  C12L2 

(4) 

S6  =  C61L1  +  C62L2’ 

and  as  a  result  of  =  F  and  L2  =  d2xF: 

S1  =  C11L1  +  C12L2 

(5) 


s'  =  C  L  +  C  I 
6  L61L1  C62L2 


If  we  denote  b^,  b2,  ....  b^  as  the  slopes  of  the  best  fit  voltage  vs. 
load  lines,  then  we  see  that: 


=  C11  +  dlC 1 2  =  bl 


and 


S» 

T  =  cn  +  d2c  1 2  =  bi 


(6) 


•=—  =  C,,  +  d.C,-.,  =  bc 
F  61  1  62  6 


S’ 

r  =  C61  +  d2C6 2  *  b6 
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Fig.  6  Loading  of  the  transducer  for  the  shear  and  moment  calibration. 


In  Eq.  (6)  there  are  12  equations  in  12  unknowns;  we  solve  for  the  unknowns, 
however,  two  at  a  time,  e.g. 


C11  +  dlC12  =  bl 


C11  +  d2C12  =  bl 


(7) 


This  process  gives  us  two  column  vectors  of  [C]  for  the  corresponding 
moment  and  shear  applications.  A  subsequent  repeat  of  this  procedure  in  a 
direction  orthogonal  to  the  one  just  performed  yields  the  next  two  column 
vectors  of  [C] .  The  fifth  and  sixth  columns  of  [C]  corresponding  to  the 
axial  and  torque  applications  completes  the  calibration  matrix. 

We  have  shown  how  the  loads  were  applied  for  the  shear  and  bending 
moment  calibrations  in  Fig.  6.  Finally,  the  axial  load  and  torque  calibra¬ 
tion  configurations  are  included  as  Figs.  7a  and  b  respectively  for  the 
sake  of  completeness. 

PRINCIPLES  OF  SONIC  DIGITIZING  AND  A  FORCE  APPLICATOR  WITH  SONIC  EMITTERS 

Before  we  discuss  some  details  of  the  force  applicator  we  will 
introduce  the  principles  of  sonic  digitizing.  In  the  current  research,  the 
kinematics  of  both  the  moving  body  segment  and  the  force  applicator  are 
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Fig.  7b  Torque  application  on  the  transducer  for  the  torque  calibration. 


monitored  by  a  Graf/Pen  Sonic  Digitizer  manufactured  by  Science  Accessories 
Corporation.  Sonic  digitizing  is  the  process  of  converting  information 
on  location  or  position  in  one,  two,  or  three  dimensions  to  digital  values 
in  a  form  suitable  for  data  transmission,  storage,  and  processing. 

The  system  used  to  accomplish  this  conversion  consists  of  two  or  three 
microphone/sensor  assemblies  (for  2-D  and  3-D  data  conversions, 
respectively),  an  electronic  control  unit  and  a  generator/multiplexer  unit 
which  is  used  to  select  and  power  eight  sonic  impulse  emitters.  To  explain 
the  digitizer's  operation,  let  us  first  consider  the  2-D  mode  of  operation. 

The  two-dimensional  microphone/sensor  assembly  consists  of  two 
perpendicular  linear  microphones  as  shown  in  Fig.  8.  These  sensors  define 
a  planar  effective  working  area  of  approximately  35  cm  x  35  cm.  The 
Graf/Pen  uses  impulses  generated  at  the  tip  of  the  emitter  to  calculate 
its  position  in  the  working  area.  The  times  required  for  the  sound  waves 
to  reach  the  two  microphone/sensors  are  converted  into  distance  measure¬ 
ments  (x  and  y  coordinates) .  The  measurements  are  then  transmitted  to  the 
computer  in  BCD  Cartesian  form.  The  2-D  mode  is  primarily  used  to  digitize 
joint  surface  curves  which  are  obtained  either  from  X-ray  projections  or 


Fig.  8  The  2-dimensional  microphone/ sensor  assembly. 
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physical  cross-sections  of  the  articulating  joints.  These  coordinate 
curves  are  subsequently  used  in  the  joint  model  analysis.  Note  that 
concurrent  with  this  present  work  a  research  program  on  dynamic  modeling 
of  human  articulating  joints  has  been  conducted  for  the  AMRL  [16,17].  The 
2-D  microphone/sensor  assembly  also  has  a  menu  capability  which  provides 
the  ability  of  simultaneous  alphanumeric  data  entry  to  the  computer  by 
the  digitizer. 

Of  course,  one  advantage  of  using  sound  as  a  ranging  device  is  that 
digitation  need  not  be  confined  to  a  plane.  The  3-D  microphone/ sensor 
unit  used  for  the  kinematic  data  acquisition,  consists  of  four  linear 
microphones  arranged  in  a  planar,  rectangular  manner.  These  microphones, 
in  turn,  define  a  3-dimensional  effective  working  volume  of  approximately 
150  cm  x  75  cm  x  180  cm  along  the  x,  y,  and  z  directions  respectively. 

The  3-D  mode  of  operation  is  similar  to  that  of  the  2-D  mode.  In  the  3-D 
set  up,  however,  the  distances  measured  are  the  slant  ranges  to  each  of 
the  four  coplanar  sensors.  Thus  the  information  generated  by  each  sensor 
represents  the  radius  of  a  circular  arc  which  includes  the  impulse  source 
(the  tip  of  the  sonic  emitter]  and  is  in  a  plane  perpendicular  to  the 
sensor.  Four  sensors  are  used  merely  for  the  purpose  of  accuracy.  In 
actual  operation,  the  digitizer  examines  the  signal  from  all  four  sensors, 
selects  the  three  smallest  signals,  and  disregards  the  fourth.  The 
location  of  the  sonic  emitter  is  then  calculated  as  being  at  the  inter¬ 
section  of  the  three  smallest  arcs.  The  three  slant  ranges  are  easily 
converted  into  Cartesian  x,  y,  and  z  coordinates  by  a  microprocessor  in 
the  control  unit  and  converted  for  transmission.  In  the  3-D  mode  of 
operation,  the  sonic  digitizer  will  monitor  up  to  seven  sonic  emitters  and 
output  their  positional  information  in  serial  form.  One  additional  sonic 
emitter  is  always  reserved  for  2-D  data  or  alphanumeric  input. 

As  stated  previously,  one  of  the  purposes  of  the  sonic  digitizer  is 
to  monitor  the  kinematics  of  the  force  applicator  used  in  this  research. 
Fig.  9  shows  the  picture  of  the  present  force  applicator,  which  is  a 
significant  improvement  over  the  multiple-linkage  force  applicator  (GFA) 
used  in  the  previously  conducted  research  [7].  The  force  transducer  which 
is  used  has  already  been  described  in  a  previous  section  of  this  report. 


SONIC  EMITTERS 


FORCE  TRANSDUCER 


Fig.  9  The  force  applicator. 


The  position  and  direction  of  force  application  are  obtained  by  a  vector 
analysis  which  proceeds  from  the  position  information  of  three  sonic 
emitters  located  on  the  force  applicator.  Note  that  there  are  more  than 
three  sonic  emitters  on  the  force  applicator.  However,  the  coordinate 
information  obtained  from  three  sonic  emitters  is  sufficient  for 
determination  of  the  position  and  direction  of  force  application. 

Aside  from  the  improvements  in  the  force  transducer  itself,  this 
force  applicator  has  several  advantages  over  the  previously  used  multiple 
linkage  system.  It  has  significantly  improved  our  mobility,  and  its  ease 
of  handling  allows  better  control  of  the  magnitude  and  direction  of  force 
application.  In  addition,  data  acquisition  is  now  controlled  by  the 
force  applicator  operator.  An  on-off  switch  has  been  installed  on  the 
front  handle  of  the  force  applicator  which  controls  the  firing  of  the 
sonic  emitters  located  on  both  the  force  applicator  and  on  the  moving 
body  segment.  The  signals  transmitted  from  the  sonic  digitizer  to  the 
computer  are  then,  in  turn,  used  to  trigger  A  to  D  conversions  of  the 
signals  originating  from  the  force  transducer.  The  digital  readings  from 
the  transducer  are  synchronized  with  the  coordinate  data  from  the  sonic 
emitters  for  input  into  data  files.  This  process  continues  until  it  is 
terminated  by  the  person  controlling  the  force  applicator. 


KINEMATICS  BY  MEANS  OF  SONIC  EMITTERS 


In  this  section  theoretical  aspects  of  the  kinematics  of  the  relative 
motion  in  three  dimensional  space  will  be  presented  by  considering  the 
motion  of  the  upper  arm  with  respect  to  the  torso  and  casting  the  equa¬ 
tions  into  a  format  suitable  for  the  utilization  of  sonic  emitters. 

The  quantitative  determination  of  the  nature  of  the  relative  motion 
between  two  body  segments  (e.g.  the  upper  arm  and  the  torso),  which  are 
connected  by  a  complex  anatomical  joint,  is  of  prime  importance  to  the 
biomechanician  as  well  as  to  those  in  medicine.  The  simple  hinge  joint 
with  one  degree  of  freedom  and  the  ball  and  socket  joint  with  three 
degrees  of  freedom  are  the  most  popular  joints  employed  in  multisegmented 
mathematical  models  of  the  human  body.  Under  both  physiological  and 
external  loads,  each  articulating  human  joint  can  be  considered  as  a 
general  joint  which  can  display  up  to  six  degrees  of  freedom,  to  some 
extent.  A  good  example  of  a  general  joint  is  the  shoulder  complex,  which 
exhibits  four  independent  articulations  among  the  humerus,  scapula, 
clavicle,  and  thorax.  Of  course  at  the  shoulder  complex,  the  six  degrees 
of  freedom  refers  to  the  motion  of  the  humerus  relative  to  the  torso.  If 
one  considers  the  total  number  of  degrees  of  freedom  for  the  motions 
executed  by  the  various  bones  of  the  shoulder  complex  one  can  easily  get 
a  number  much  higher  than  six  even  with  the  proper  consideration  of  vari¬ 
ous  constraints  present  in  the  joint  complex. 

The  basic  concept  for  the  study  of  the  relative  motion  between  two 
body  segments  is  quite  fundamental  and  takes  its  impetus  from  the  study 
of  the  general  motion  of  two  rigid  bodies,  say  A  and  B,  in  three- 
dimensional  space.  The  relative  motion  of  body  B  with  respect  to  body  A 
can  be  characterized  by  a  unique  axis  called  the  screw  axis.  The  relative 
displacement  of  body  B  from  one  position  to  another  can  be  defined  in 
terms  of  a  rotation  a  about,  the  translation  s  along  the  screw  axis.  For 
each  incremental  displacement  of  body  B  with  respect  to  body  A,  a  new 
screw  axis  is  defined.  In  fact,  if  the  displacements  taken  by  body  B  are 
made  infinitesimal  in  size,  the  collection  of  the  screw  axes  will  form  a 
ruled  surface  called  an  axode.  There  are  two  unique  axodes,  one  associated 
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with  the  motion  of  body  B  with  respect  to  body  A,  and  another  one  associated 
with  the  motion  of  A  with  respect  to  B.  During  the  relative  motion,  the 
two  axodes  roll  and  slide  relative  to  each  other  along  a  generator  which 
is  momentarily  common  to  both  axode  surfaces  [18].  For  a  ball  and  socket 
joint,  the  axodes  become  two  rolling  cones  with  common  apexes  at  the 
center  of  the  sphere;  for  a  hinge  joint,  they  simply  degenerate  to  a  single 
axis,  i.e.  the  axis  of  the  joint.  Of  course,  in  each  one  of  these 
specialized  cases  the  translation  s  along  the  screw  axis  becomes  zero. 

Let  us  consider  the  relative  motion  of  arm  designated  as  body  segment 
B  in  Fig.  10  with  respect  to  the  torso  which  is  designated  as  body  segment 
A.  For  the  intention  of  studying  the  relative  motion  between  these  two 


Body  Segment  B 


Fig.  10  Motion  of  body  segment  B  with  respect  to  body  segment  A  in  three 
dimensional  space.  Displacement  of  body  segment  B  from  position  1  to 
position  2  is  illustrated. 


body  segments,  let  us  assume  that  one  of  them  e.g.,  body  segment  A,  be 
fixed  and  the  other  one,  body  segment  B,  is  moving  relative  to  body 
segment  A.  Note  that  the  range  of  motion  between  body  segments  A  and  B 
is  controlled  by  the  joint  anatomy  and  ligamentous  as  well  as  muscle 
forces  present  during  motion.  It  is  very  essential  to  point  out  before 
we  proceed  further  that  no  human  body  segment  is  a  rigid  body  in  the 
sense  defined  in  mechanics.  However,  for  the  purposes  of  studying  the 
nature  of  motion  in  a  given  anatomical  joint  the  adjacent  body  segments 
to  this  joint  can  be  assumed  to  be  rigid  bodies  if  certain  precautions  are 
observed.  Let  the  unprimed  xyz  and  the  primed  x'y'z'  cartesian  coordinate 
systems  be  attached  to  the  body  segments  A  and  B,  respectively.  These 
coordinate  systems  can  be  also  referred  to  as  fixed  and  moving  coordinate 
systems  since  body  segment  A  is  assumed  to  be  fixed  and  body  segment  B  is 
considered  to  be  moving. 

Let  P  be  the  point  designating  one  of  the  sonic  emitters  attached  to 
the  body  segment  B  whose  position  is  indicated  by  vectors  p  and  r^  (f^  @ 
position  1;  r2  @  position  2)  in  references  x'y'z'  and  xyz,  respectively. 
Since  point  P  is  fixed  in  a  moving  body  segment,  the  components  of  vector 
p  in  reference  x'y'z'  will  remain  constant,  whereas  the  components  of 
vector  rj  will  change  as  the  body  segment  B  moves  relative  to  body  segment 
A.  In  fact,  the  relative  position  of  the  moving  segment  is  considered 
completely  determined  with  respect  to  the  fixed  segment  if  for  every  point 
in  the  moving  segment  and  its  associated  local  position  vector  p,  the 
corresponding  r^  can  be  found.  From  Fig.  10  the  relationship  between 
vectors  r^,  fQ  and  p  can  be  written  in  compact  matrix  form 

trj]  -  [r0]  +  [T]  [p]  (8) 

where  the  elements  of  matrix  T,t..  are  the  direction  cosines  of  the  O'x', 
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O'y'  and  O'z'  axes  relative  to  the  axes  of  the  xyz  reference.  Thus, 
matrix  T  represents  rotational  orientation  of  x'y'z',  whereas  matrix  [rQ] 
represents  separation  of  x'y'z'  with  respect  to  the  fixed  xyz  reference. 

It  is  more  convenient  to  express  Eq.  (8)  in  terms  of  augmented  vectors 
r,  ,  p  whose  first  component  is  1,  and  a  single  4x4  matrix  T  by 
adding  the  equation  1  =  1  to  the  system  of  equations  contained  in  Eq.  (8). 
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Thus ,  Eq . 


(8)  takes  the  following  expanded  and  compact  forms: 
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Note  that  there  is  a  new  matrix  T  for  each  position  of  body  segment  B, 
however,  at  a  given  position  the  matrix  T&  is  the  same  for  all  points  of 
body  segment  B.  Hence,  determination  of  the  matrix  T  is  sufficient  to 
know  the  position  of  the  moving  body  segment  relative  to  the  fixed  one. 

For  the  total  description  of  the  relative  motion  between  two  body 
segments,  besides  the  knowledge  of  the  instantaneous  positions  of  the  one 
body  segment  relative  to  another,  we  must  have  a  description  of  the  nature 
of  the  displacement  during  the  motion  of  the  moving  segment  from  position 
1  to  position  2.  Needless  to  say,  instantaneous  positions  and  displace¬ 
ments  are  very  closely  related.  Displacement  analysis  results  in 
determination  of  the  set  of  screw  axes;  thus,  it  is  essential  for  the 
accurate  determination  of  the  locations  of  the  joint  centers  (i.e.  centers 
of  rotation) .  As  explained  at  the  beginning  part  of  this  section, 
displacement  of  body  segment  B  from  position  1  to  position  2  can  be  defined 
in  terms  of  a  rotation  a  about,  the  translation  s  along  the  screw  axis. 

Fig.  11.  If  the  two  positions  considered  are  close  enough,  then  the 
displacements  are  incremental  and  refer  to  a  new  screw  axis  each  time. 

Let  us  again  consider  the  point  P  which  defines  the  location  of  one 
of  the  sonic  emitters  fixed  to  the  body  segment  B.  As  the  body  segment 
B  moves  from  position  1  to  position  2,  we  can  define  the  position  vectors 
and  ?2  which  are  expressed  in  terms  of  components  in  xyz  references: 


=  Xji  +  y^j  +  Zjk  and  =  X2i  +  y2^  +  (10) 

The  components  of  x^  will  be  related  to  the  components  of  r^  by  the 
following  special  case  of  an  affine  transformation: 
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Screw  Axis 


Fig.  11  Representation  of  the  motion  of  body  segment  B  from  position  1  to 
position  2  in  terms  of  a  rotation  a  about  and  translation  s  along  the 
screw  axis.  A  pure  rotation  of  magnitude  a  of  a  point  Q  about  an  axis 
parallel  to  the  screw  axis  is  also  displayed. 
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Eq.  (11)  can  be  written  in  augmented  form  as 
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Eqs.  (11)  and  (12)  can  be  expressed  in  a  more  compact  form  as 
[r2]  =  [MHrj]  +  [a]  and  [r2a]  =  [M&]  [rla] , 


(13a, b) 


respectively.  From  a  mathematical  point  of  view,  the  augmented  form 
describes  the  motion  of  a  three-dimensional  set  of  points  contained  within 
body  segment  B  relative  to  a  hyperplane  [19].  In  Eq.  (13a)  matrix  M 
represents  a  pure  rotation  about  the  screw  axis  and  also  contains  direction 
cosines  information  for  the  screw  axis.  The  augmented  matrix  in 
Eqs.  (12)  and  (13b)  has  the  total  information  for  determination  of  the 
displacement  of  the  moving  body. 

To  determine  the  matrix  M  ,  it  is  sufficient  to  know  the  coordinates 

a 

of  four  non-coplanar  points  Pj,  P2,  P3  and  P4  in  the  moving  body  segment 
(or  on  the  extension  of  the  moving  body  segment)  in  both  positions  1  and  2. 
Four  sonic  emitters  attached  to  the  moving  body  segment  represent  these 
points  whose  coordinates  can  be  monitored  continuously.  Let  us  designate 
the  position  vectors  of  these  four  sonic  emitters  by  p . .  and  r. .  relative 
to  the  moving  x'y'z'  and  the  fixed  xyz  references,  respectively.  In  this 
subscripted  notation,  the  first  subscript  i  (i  =  1,2, 3, 4)  is  used  to 
designate  an  emitter  and  the  second  subscript  j  (j  *  1,2)  is  used  for  the 
designation  of  positions  1  and  2.  The  augmented  vectors  r„a  for  each 
emitter  can  be  related  via  fcq.  (13b), 

tri2J  '  ["a^iial  Ci  =  1,2,3,4)  (14) 

Eq.  (14)  represents  four  matrix  equations  which  can  be  written  in  one 

equation  by  defining  a  new  matrix  A  whose  columns  are  the  augmented 

vectors,  r.  .  , 
ija’ 
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Xlj  X2j  X3j  X4j 

ylj  y2j  y3j  y4j 

z,  .  z. .  z„ .  z  . . 

L1  2j  jj  4j 

With  the  definition  given  by  Eq.  (15),  Eq.  (14)  becomes 

[A2]  =  [Ma][Aj]  (16) 

The  matrix  A  has  a  non -vanishing  determinant  since  it  is  formed  by  the 
position  vectors  of  non-coplanar  points  (sonic  emitters);  thus,  it  has  an 
inverse  and  the  desired  matrix  M  is  easily  obtained  from  Eq.  (16), 

[Ma]  =  [A2] [Aj] -1  (17) 

When  the  matrix  M  is  determined  by  means  of  Eq.  (17)  and  data  obtained 
from  four  sonic  emitters  attached  to  the  moving  body  segment,  information 
such  as  direction  cosines  of  the  screw  axis,  rotation  a  about  and  transla¬ 
tion  s  along  the  screw  axis,  and  the  coordinates  of  the  piercing  point  of 
the  screw  axis  with  one  of  the  reference  planes  can  be  obtained  according 
to  the  theory  presented  previously,  [7,10], 

DESCRIPTION  OF  SOFTWARE  DEVELOPMENT 

Software  is  dependent,  to  various  extents,  upon  the  architecture  of 
the  particular  computer  on  which  it  is  to  run.  A  Digital  Equipment 
Corporation  PDP  11/34A  was  selected  for  this  research  project.  The 
software  for  this  machine  included  an  RSX-11M  operating  system,  a 
Fortran-IV  compiler,  a  scientific  subroutine  package  and  interface  and 
applications  programs.  In  this  section  we  will  discuss  in  some  detail, 
the  system  generation  process,  the  force  transducer  calibration  software, 
and  the  data  acquisition  and  analysis  programs. 

SYSTEM  GENERATION 

Generation  of  the  operating  system  required  considerable  planning  and 
consultation  with  Digital  representatives  in  order  to  tailor  the  system  to 
the  project's  needs.  An  important  and  troublesome  component  of  the  system 
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software  is  the  set  of  routines  to  manage  the  operation  of  peripherals, 
including  the  Data  Translation  model  DT  1712/5714-64,  14-bit  analog- 
digital  converter  and  the  Science  Accessories  Corporation  model  GP6-3D-3 
sonic  digitizer.  These  major  components  of  the  experimental  apparatus 
are  considered  "foreign"  equipment  and  could  not  be  installed  with  the 
system  by  the  Digital  Equipment  Corporation  representatives.  This 
necessitated  considerable  time  and  study  to  become  familiar  with  the  entire 
operating  system  and  hardware.  Provided  with  a  basic  outline  of  a  software 
interface  for  an  A-D  board,  a  subprogram  was  developed  to  address  the 
A-D  board,  to  identify  the  channel  to  be  read,  to  supply  any  "software" 
gain  to  the  signal  and  to  check  if  the  A-D  conversion  was  successful. 

Because  of  a  foreign  hardware  interface  involvement,  the  subprogram 
had  to  be  written  in  MACRO-11,  the  assembly  language  of  the  PDP-11.  This 
necessitated  the  learning  of  an  entirely  new  programming  language  since  we 
had  prior  experience  only  in  Fortran.  The  final  requirement  was  that  this 
subprogram  be  callable  from  a  Fortran  program.  The  result  of  this  effort 
is  the  subprogram  entitled  (OSUATD) ,  a  listing  of  which  is  included  in 
both  the  transducer  calibration  program  (CHNCAL) ,  and  the  data  acquisition 
program  (BIGQIO),  found  in  Appendices  A  and  B. 

TRANSDUCER  CALIBRATION  PROGRAM 

The  transducer  calibration  program  (CHNCAL),  a  listing  of  which  is 
contained  in  Appendix  A,  is  short  and  straightforward.  (CHNCAL)  computes 
the  total,  mean,  standard  deviation  and  the  minimum  and  maximum  values  of 
1000  samples  read  from  each  channel  of  the  A-D  board.  All  non-zero  values 
are  analyzed,  totals  are  accumulated  and  the  minimum  and  maximum  values 
are  found.  Following  this,  the  mean  and  standard  deviation  are  calculated. 
The  divisor  for  the  standard  deviation  calculation  is  one  less  than  the 
number  of  samples  taken. 

Output  includes  the  channel  number  being  sampled,  the  mean,  the 
minimum  and  maximum  values  read,  and  a  status  code.  If  errors  are 
encountered,  appropriate  descriptive  messages  are  displayed.  The  program 
may  be  run  stand-alone  for  single  sweeps  of  the  A-D  channel (s),  or  it  may 
be  non  from  an  indirect  command  file.  The  indirect  command  file,  a 
separate,  monitor  control  routine,  provides  prompts  for  test  identification 
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and  executes  (CHNCAL)  each  time  a  calibration  test  is  set  up,  thus 
eliminating  the  need  for  manually  installing  and  executing  the  program  for 
each  test  set-up.  The  command  file  (TRNCAL)  prompts  for  the  name  of  the 
calibration  test  being  performed,  displays  the  date  and  time  of  the  test, 
prompts  for  the  amount  of  the  load  applied  to  the  transducer  and  executes 
(CHNCAL),  repeating  the  date,  time  and  applied  load  prompts  until  the 
testing  sequence  is  completed.  A  listing  of  (TRNCAL)  is  included  in 
Appendix  A. 

Two  subroutines  are  called  in  (CHNCAL).  The  subroutine  (OSUATD) , 
described  briefly  above,  obtains  the  A-D  readings,  checks  for  successful 
conversion,  and  applies  gain  if  desired.  It  is  called  1000  times  for  each 
channel  tested.  (TALLY)  is  the  statistical  subroutine  which  takes  the 
1000  readings  for  a  channel  and  calculates  the  total,  mean,  the  minimum 
and  maximum  values  and  checks  for  null  or  single-value  input.  It  is 
called  once  for  each  channel  tested.  The  program  tests  channels  32  to  37 
of  the  A-D  board,  but  may  be  easily  edited  for  other  sampling 
configurations. 

DATA  ACQUISITION  PROGRAM 

The  data  acquisition  program,  (BIGQI0) ,  is  designed  to  create  an 
unformatted  data  file  containing  digitized  points  in  space  and  concurrent 
transducer  readings.  A  complete  listing  of  this  program  is  contained  in 
Appendix  B.  The  program  queues  from  the  Science  Accessories  sonic 
digitizer.  Each  point  is  described  by  a  set  of  three  coordinates,  relating 
the  position  to  the  origin  of  the  3-D  sonic  tablet.  Seven  individual 
points  are  monitored  throughout  the  testing.  The  first  three  points  are 
obtained  from  sonic  emitters  mounted  on  the  force  applicator.  The  remain¬ 
ing  four  points  are  mounted,  by  use  of  an  elastic  cuff,  on  the  moving 
body  segment  of  the  subject.  Thus,  each  scan  of  the  digitizer  monitors 
the  positions  of  these  seven  points. 

The  sonic  digitizer  can  scan  its  seven  channels  up  to  eight  times  each 
second.  Since  this  is,  in  computer  reference,  a  long  time,  the  sonic 
digitizer  is  considered  to  be  the  limiting  factor  in  the  rate  of  data 
acquisition.  It  is  for  this  reason  that  the  data  acquisition  program 
queues  from  the  digitizer.  During  the  time  between  the  sampling  of  the 
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digitized  points,  the  computer  scans  each  of  the  six  transducer  channels 
100  times  and  takes  the  mean  of  each  channel's  readings,  where  each 
channel  of  the  A-D  board  represents  a  separate  strain  gage  bridge  of  the 
force  transducer. 

(BIGQIO)  prompts  for  a  name  to  be  given  to  the  data  file.  The  format 
for  the  file  name  is  [SIX  CHARACTER  STRING] ,  [THREE  CHARACTER  STRING] , 
[SINGLE  DIGIT],  (e.g.:  SH0LD.DAT;!).  The  file  name  is  used  later  to 
identify  a  particular  test;  configuration  from  the  test  log,  which  is  kept 
separately  and  contains  a  description  of  the  test  position,  the  subject's 
identification,  the  3-D  sonic  tablet  position  relative  to  the  body 
coordinate  system  and  the  data  file  name. 

The  file  structure  is  currently  set  for  one  file  per  test.  Each  file 
may  contain  a  variable  number  of  records.  Each  record  contains  a  complete 
set  of  digitized  points  and  transducer  readings.  (BIGQIO)  prompts  for  the 
number  of  records  to  be  allocated  for  each  test.  The  number  entered  in 
response  to  this  prompt  should  be  approximately  five  times  the  duration  of 
the  test  (in  seconds) .  For  an  average  ten  second  test,  50  should  be 
entered  for  the  number  of  records,  where  each  record  contains  one  sweep  of 
the  seven  sonic  emitters  and  the  mean  of  700  readings  for  each  of  the  six 
transducer  channels. 

The  program  begins  data  collection  whenever  the  first  sonic  emitter 
coordinates  are  received.  (BIGQIO)  will  halt  if  the  number  of  records 
allocated  is  exceeded,  if  the  A-D  board  fails  to  convert  a  signal  or  if  the 
digitizer  sends  a  series  of  zeros.  Since  the  number  of  records  should  be 
generously  estimated,  the  collection  of  data  should  be  completed  before 
the  record  allocation  is  exceeded. 

Two  subroutines  are  utilized  in  (BIGQIO).  The  system  library  routine 
(ASSIGN)  is  used  to  identify  the  data  file,  assigning  a  logical  unit  number 
and  name  to  the  file.  The  second  subroutine  (OSUATD)  has  been  described 
above. 

DATA  ANALYSIS  PROGRAM 

The  program  (TABLE)  is  designed  to  read  data  from  disk  storage  and 
analyze  it,  printing  the  results  for  each  position  during  a  particular 
test.  The  program  is  interactive,  requiring  manual  input  from  the 
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terminal.  The  prompts  for  input  are  self-explanatory.  Each  prompt 
includes  the  number  of  permissible  characters  for  each  entry.  Since  a 
series  of  tests  may  be  performed  for  a  given  chair  position,  prompts  are 
divided  into  general  position  and  identification  prompts,  and  test 
specific  prompts.  The  position  and  identification  prompts  are: 

(a)  (Enter  the  name  of  joint  tested  [ S-9 ] : )  Enter  up  to  nine 
alphanumeric  characters  for  use  in  the  output  heading. 

(b)  (Enter  subject  name  or  number  [S-25]:)  Enter  up  to  25 
alphanumeric  characters  for  identification  on  the  output 
heading. 

(c)  (Enter  the  coordinates  of  the  3-D  board  origin,  with  respect 
to  the  fixed  body  segment  origin:)  Enter  the  X,  Y  and  Z 
coordinates  in  decimal  numbers. 

(d)  (Enter,  in  degrees,  the  rotation  of  the  3-D  board  from  the 
vertical  [N-8]:)  Enter  the  appropriate  decimal  number. 

(e)  (Enter,  in  degrees,  the  rotation  of  the  3-D  board  from  the 
X-Y  plane  [N-8]:)  Enter  a  decimal  number. 

(f)  (Enter  the  coordinates  of  the  fixed  joint  center  with  respect 
to  the  fixed  body  segment  origin:)  Enter  the  appropriate 
coordinates  as  decimal  numbers. 

(g)  (Enter  the  length  to  the  JOINT  CENTER  [N-8] :)  Enter  a  decimal 
number . 

(h)  (Enter  the  tilt  of  the  chair,  in  degrees  [N-8]:)  Enter  a 
decimal  number. 

(i)  (Enter  the  rotation  of  the  chair,  in  degrees  [N-8] :)  Enter  a 
decimal  number. 

For  each  test  under  the  above  board  and  chair  positions  for  the 
specified  joint  and  subject,  the  following  prompts  are  issued: 

(a)  (Enter  any  messages  (Maximum  -  80  characters):)  Enter  up  to 
80  characters  which  describe  the  specific  test. 

(b)  (Enter  data  file  name  [S- 13] :)  Enter  up  to  13  characters  which 
identify  the  file  containing  the  data  for  the  test  (e.g.: 
SH0LD.DAT;  1). 
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(c)  (Enter  the  number  of  records  to  be  read  [N-5] : )  Enter  the 

number  of  records  contained  in  the  data  file  (this  information 
was  output  from  (BIGQIO).) 

Data  analysis  program  outputs  to  the  hardcopy  terminal  and  includes: 
a  header  identifying  the  test  and  subject;  the  calculated  spherical 
coordinates,  9  and  <fi,  of  the  moving  body  segment;  the  force  and  moment 
components  and  magnitudes;  the  calculated  fixed  joint  center  and  the 
distances  between  various  fixed  sonic  emitters.  The  distances  between 
fixed  sonic  emitters,  such  as  those  on  the  force  applicator,  should  remain 
constant  and  thus  serve  as  a  check  on  the  validity  of  the  data. 

All  calculations  are  performed  with  respect  to  the  3-D  board 
coordinate  system.  The  results  are  then  transformed  to  the  fixed  body 
coordinate  system.  Therefore,  the  transformation  matrix  will  be  dealt 
with  last.  We  next  provide  a  brief  description  of  the  subroutines  of  data 
analysis  program. 


Point  of  Force  Application  [SUBROUTINE  FORPT] 

The  first  calculations  are  to  determine  the  point  of  application  of 
the  force.  This  subroutine  is  called  once  for  every  set  of  data  points. 

It  uses  the  coordinates  of  the  first  three  sonic  emitters  to  determine  the 
direction  and  position  of  the  point  of  force  application. 

Moving  Body  Segment  Position  [SUBROUTINE  POSITN] 

The  position  of  the  moving  body  segment  is  necessary  for  calculation 
of  the  joint  center.  This  subroutine  calculates  the  position  from  the  four 
sonic  emitters  located  on  the  moving  body  segment.  The  process  is  similar 
to  locating  the  point  of  force  application. 

Spherical  Coordinates,  9,  <p  (SUBROUTINE  SPHERE] 

The  location  of  the  moving  body  segment  is  located  by  the  spherical 
coordinates,  9  and  <p.  This  subroutine  calculates  these  angles  when  called 
from  (SUBROUTINE  JTCNT) .  The  vector  between  the  positions  of  the  moving 
body  segment  toward  the  joint  center  is  first  transformed  to  the  fixed 
body  segment  coordinate  system  and  9  and  4>  are  calculated  directly. 


Fixed  Joint  Center  [SUBROUTINE  F JTCNT ] 

Given  the  input  of  the  coordinates  of  a  fixed  position  at  the  joint 
center,  this  subroutine  calculates  a  relatively  fixed  joint  center  located 
at  the  end  of  the  long  bone  axis  of  the  moving  body  segment. 

Matrix  Inversion  [SUBROUTINE  MINV] 

This  subroutine  inverts  a  matrix  using  the  standard  Gauss-Jordan 
method.  The  determinant  is  also  calculated.  A  determinant  of  zero 
indicates  that  the  matrix  is  singular. 

Matrix  Multiplication  f SUBROUTINE  GMPRD] 

This  subroutine  multiplies  two  matrices  to  form  a  resultant  general 
matrix. 

Transformation  Matrix  [SUBROUTINE  MATRIX] 

When  all  the  positions,  locations  and  values  have  been  calculated, 
the  reference  coordinate  system  (3-D  board),  must  be  related  to  the  fixed 
body  coordinate  system.  This  subroutine  rotates  and  translates  the  3-D 
coordinate  system  into  the  fixed  body  segment  coordinate  system,  using  the 
angles  of  the  3-D  board  (tilt  and  rotation) ,  the  position  of  the  respective 
origins  and  the  angles  of  the  chair  (tilt  and  rotation) . 

DATA  COLLECTION  ON  HUMAN  SUBJECTS 

The  main  thrust  of  the  research  program  reported  here  is  the 
collection  of  the  active  muscle  force  and  moment  data  on  the  higher  and 
lower  extremities  of  the  human  body.  In  this  section  of  the  report, 
first,  anthropometry  of  the  subjects  and  later  a  detailed  description  of 
the  experiments  are  presented. 

ANTHROPOMETRY  OF  SUBJECTS 

In  this  research  program,  six  healthy  subjects  (3  male,  3  female) 
between  the  ages  of  20-22  were  tested.  Subjects  were  university  students 
with  no  special  training  in  athletics.  Selected  anthropometric  measure¬ 
ments  of  these  six  subjects  are  given  in  Table  2  for  male  subjects  and  in 
Table  3  for  female  subjects.  The  body  segment  circumferences  given  in 
these  two  tables  were  obtained  according  to  the  following  definitions 
[20]: 
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(a)  Shoulder  Circumference:  maximum  horizontal  circumference  over 
the  deltoid  muscle. 

(b)  Chest  Circumference:  horizontal  circumference  at  nipple  level 
during  normal  breathing. 

(c)  Waist  Circumference:  horizontal  circumference  at  level  of  the 
greatest  lateral  indentation  of  trunk. 

(d)  Wrist  Circumference:  minimum  circumference  above  the  wrist 
bones  (styloid  processes  of  radius  and  ulna). 

(e)  Lower  Arm  Circumference:  maximum  circumference  wherever  found, 
with  the  upper  arm  horizontal,  forearm  vertical,  and  the  elbow 
90°,  muscles  tensed. 

(f)  Elbow  Circumference:  circumference  encompassing  the  elbow  tip 
and  the  elbow  crotch,  with  the  upper  arm  horizontal,  forearm 
vertical,  and  the  elbow  at  90°,  muscles  tensed. 

(g)  Biceps  Circumference:  maximum  circumference  with  elbow  bent  at 
90°  and  biceps  maximally  flexed. 

(h)  Arm,  Axillary  Circumference:  horizontal  circumference  at  the 
armpit,  arms  hanging  loosely  at  the  side. 

(i)  Arm,  Scye  Circumference:  circumference  encompassing  the  highest 
skeletal  point  on  the  lateral  edge  of  the  acromion  process  and 
the  highest  point  in  the  axilla  (armpit). 

(j)  Thigh,  Upper  Circumference:  horizontal  circumference  just 
below  the  gluteal  furrow. 

(k)  Thigh,  Lower  Circumference:  horizontal  circumference  just  above 
the  knee. 

(l)  Calf  Circumference:  maximum  horizontal  circumference,  wherever 
found. 

(m)  Ankle  Circumference:  minimum  horizontal  circumference  above  the 
projections  of  the  ankle  bones  (external  and  internal  malleoli). 

The  following  definitions  were  used  for  the  linear  dimensions  given 
in  Tables  2  and  3: 

(a)  Forearm-Wrist  Length:  distance  from  tip  of  elbow  to  the  styloid 
process  of  radius. 


Table  2 


SELECTED  ANTHROPOMETRY  OF  MALE  SUBJECTS  1,  2,  3. 
(All  length  dimensions  are  in  cm) 


1 

2 

3 

663 

756 

757 

Stature 

170 

171.2 

185.5 

Shoulder  Circumference 

113 

112.5 

111.2 

Chest  Circumference 

94.3 

100.2 

97.0 

Waist  Circumf erence 

75.0 

87.4 

81.1 

Wrist  Circumference 

16.6 

17.5 

16.0 

Lower  Arm  Circumference 

30.5 

29.5 

28.6 

Elbow  Circumference 

31.6 

32.5 

31.4 

Biceps  Circumference 

32.1 

33.5 

31.5 

Arm,  Axillary  Circumference 

31.4 

33.8 

29.8 

Arm,  Scye  Circumference 

42.1 

48.6 

45.6 

Thigh,  Upper  Circumference 

54.6 

57.2 

57.8 

Thigh,  Lower  Circumference 

39.0 

42.5 

39.5 

Calf  Circumference 

35.4 

39.6 

38.5 

Ankle  Circumference 

21.5 

24.0 

22.5 

Forearm  -  wrist  length 

25.2 

26.5 

28.5 

Shoulder  -  elbow  length 

34.2 

34.4 

39.4 

Shoulder  -  height,  sitting 

58.0 

58.8 

60.5 

Sitting  height 

90.1 

89.6 

92.0 

Shoulder  breadth 

46.7 

48.5 

47.1 

Chest  breadth 

28.8 

32.0 

30.5 

Chest  depth 

19.0 

21.0 

20.6 

Waist  depth 

18.9 

21.5 

19.6 

Buttock  -  knee  length 

57.1 

59.8 

62.1 

Buttock  -  popliteal  length 

41.1 

48.2 

52.0 

Knee  height,  sitting 

51.8 

51.2 

57.9 

Elbow  -  to  -  elbow  breadth 

46.5 

49.5 

47.2 

Hip  breadth,  sitting 

35.0 

36.5 

35.6 

Knee  -  to  -  knee  breadth,  sitting 

19.1 

22.5 

19.9 

Leg  length  (tibial  point  -  Malleo.) 

37.8 

37.5 

45.0 

Trochanterion  height 

87.0 

88.0 

98.0 
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Table  3 


SELECTED  ANTHROPOMETRY  OF  FEMALE  SUBJECTS  1,  2,  3. 
(All  length  dimensions  are  In  cm) 


1 

2 

3 

Weight  (Newtons) 

578 

592 

645 

Stature 

163.2 

174 

169.5 

Shoulder  Circumference 

101.3 

101.5 

107.6 

Chest  Circumf erence 

84.5 

88.5 

98.0 

Waist  Circumference 

74.3 

70.5 

78.1 

Wrist  Circumference 

15.5 

15.6 

16.4 

Lower  Arm  Circumference 

25.5 

24.3 

26.4 

Elbow  Circumference 

26.7 

27.8 

28.0 

Biceps  Circumference 

29.3 

27.6 

29.4 

Arm,  Axillary  Circumference 

28.7 

28.5 

31.6 

Arm,  Scye  Circumference 

39.1 

39.4 

39.8 

Thigh,  Upper  Circumference 

60.0 

56.5 

62.2 

Thigh,  Lower  Circumference 

42.1 

38.5 

43.1 

Calf  Circumference 

38.2 

37.8 

37.9 

Ankle  Circumference 

22.9 

22.6 

22.6 

Forearm  -  wrist  length 

25.0 

27.5 

25.9 

Shoulder  -  elbow  length 

32.2 

35.3 

34.9 

Shoulder  -  height,  sitting 

61.8 

58.6 

61.5 

Sitting  height 

89.2 

90.2 

90.4 

Shoulder  breadth 

43.7 

43.2 

44.5 

Chest  breadth 

27.2 

28.0 

28.6 

Chest  depth 

19.0 

19.0 

20.2 

Waist  depth 

17.4 

17.3 

18.2 

Buttock  -  knee  length 

57.1 

60.1 

57.6 

48.3 

48.8 

48.0 

Knee  height,  sitting 

49.8 

52.0 

52.1 

Elbow  -  to  -  elbow  breadth 

41.0 

39.6 

42.2 

Hip  breadth,  sitting 

41.6 

40.2 

41.5 

Knee  -  to  -  knee  breadth,  sitting 

17.0 

15.0 

16.8 

Leg  length  (tibial  point  -  Malleo.) 

36.4 

39.6 

39.0 

Trochanterion  height 

85.0 

90.0 

85.5 
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(b)  Shoulder-Elbow  Length:  distance  from  the  top  of  the  acromion 
process  to  the  bottom  of  the  elbow  when  subject  sits  erect, 
upper  arm  vertical  at  side  and  the  forearm  making  a  right  angle 
with  upper  arm. 

(c)  Shoulder  Height,  Sitting:  vertical  distance  from  the  sitting 
surface  to  the  most  lateral  point  on  the  superior  surface  of 
the  acromion  process  of  the  scapula. 

(d)  Sitting  Height:  vertical  distance  from  the  sitting  surface  to 
the  top  of  the  head. 

(e)  Shoulder  Breadth:  maximum  horizontal  distance  across  the 
deltoid  muscles  when  subject  sits  erect,  upper  arms  vertical  at 
sides  and  forearms  extended  horizontally. 

(f)  Chest  Breadth:  horizontal  distance  across  the  chest  at  nipple 
level  when  subject  stands  erect  and  breathing  normally. 

(g)  Chest  Depth:  horizontal  distance  from  front  to  back  of  chest 

at  nipple  level  when  subject  stands  erect  and  breathing  normally. 

(h)  Waist  Depth:  horizontal  distance  between  the  back  and  abdomen 
at  the  level  of  the  greatest  lateral  indentation  at  the  waist 
when  subject  stands  erect  with  relaxed  abdomen. 

(i)  Buttock-Knee  Length:  horizontal  distance  from  the  plane  of  the 
rearmost  point  on  the  buttocks  to  the  front  of  the  knee  when 
subject  sits  erect  with  knees  and  ankles  at  right  angles. 

(j)  Buttock-Popliteal  Length:  horizontal  distance  from  the  plane 
of  the  rearmost  point  on  the  buttocks  to  the  back  of  the  lower 
leg  at  the  knee  when  subject  sits  erect  with  knees  and  ankles 
at  right  angles. 

(k)  Knee  Height,  Sitting:  vertical  distance  from  the  floor  to  the 
uppermost  point  on  the  knee  when  subject  sits  erect  with  knees 
and  ankles  at  right  angles. 

(l)  Elbow- to-elbow  Breadth:  maximum  horizontal  distance  across  the 
lateral  surface  of  the  elbows  when  subject  sits  erect  with  upper 
arms  vertical  and  touching  the  sides,  and  forearms  extended 
horizontally. 


•rr 


(m)  Hip  Breadth,  Sitting:  maximum  horizontal  distance  across  the 
hips  when  subject  sits  erect  with  knees  and  ankles  at  right 
angles. 

(n)  Knee-to-knee  Breadth,  Sitting:  maximum  horizontal  distance 
across  the  lateral  surfaces  of  the  knee  when  subject  sits  erect 
with  knees  at  right  angles  and  pressed  together  lightly. 

(o)  Leg  Length:  distance  between  the  tibial  point  and  malleolar 
point. 

(p)  Trochanteric  Height:  vertical  distance  between  the  trochanteric 
point  and  floor  when  subject  stands  erect. 

[hiring  anthropometric  measurements  and  subsequent  experiments  male 
subjects  wore  swimming  trunks,  female  subjects  wore  gymnast's  leotards. 

DESCRIPTION  OF  EXPERIMENTS 

Two  sets  of  experiments  on  the  collection  of  active  force  and  moment 
data  were  considered.  The  first  set  of  experiments  was  designed  to 
determine  the  active  muscle  resistances  and  corresponding  moment  values 
at  the  shoulder  and  hip  joints  of  subjects  when  there  is  a  strong  tendency 
of  dislodging  of  the  upper  and  lower  extremities  from  their  initial 
(stowed)  positions.  These  experiments  included  seat  armrest  initiation 
force  simulation,  D-ring  force  simulation,  and  face  curtain  force  simulation. 
Brief  description  of  these  force  applications  are  given  next  by  the  aid 
of  a  coordinate  system  for  which  positive  x-direction  refers  to  anterior 
direction,  positive  y-direction  to  right  lateral  and,  naturally,  positive 
z-direction  refers  to  inferior  direction  if  the  origin  of  the  cartesian 
coordinate  system  is  fixed  in  the  mid-torso  location.  In  Fig.  12  this 
coordinate  system  is  shown  as  (xt,yt,zt).  The  angles  and  ^2>^2 

shown  in  Fig.  12  define  the  angular  orientation  of  the  subject  restraint 
system  and  the  microphone  assembly,  respectively,  with  respect  to  the 
laboratory  coordinate  system  (X,Y,Z). 

In  the  experiment  entitled  "seat  armrest  initiation  force  simulation" 
we  determined  the  magnitudes  of  forces  necessary  to  separate  the  lower 
arm  from  the  initial  seat  armrest  position.  For  the  first  part  of  the 
experiment  the  force  cuff  of  the  force  applicator  was  placed  along  the 
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Fig.  12  Schematic  drawing  of  experimental  setup. 

y-direction  on  the  lower  arm  and  on  the  wrist,  as  shown  in  Figs.  14  fj  15, 
thus,  subjects  applied  an  isometric  force,  F^,  along  the  negative 
y-direction,  i.e.  medially,  on  the  force  transducer.  In  the  second  part 
of  the  experiment  the  force  application  by  the  subjects  was  along  the 
positive  z-direction  as  shown  in  Figs.  16  $  17,  and  these  forces  are 
tabulated  as  in  results  section.  In  the  "D-ring  force  simulation" 
experiment  the  subjects  as  shown  in  Figs.  18  8  19,  with  extended  arm  and 
hand  near  the  crotch  applied  an  isometric  force  along  the  negative 
y-direction  on  the  force  transducer.  In  the  "face  curtain  force 
simulation"  experiment  the  subjects  with  a  flexed  elbow  brought  their 
hands  toward  the  lateral  side  of  their  face  and  at  this  position  they 
applied  an  isometric  force  along  the  negative  y-direction  on  the  force  cuff 


Fig.  14.  Force  application  on  the  wrist  along  the  negative  y  direction 
by  one  of  the  male  subjects  (MSI). 
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Fig.  16.  Force  application  along  the  positive  z  direction  by  one  of 
the  male  subjects  (MSI) — Seat  armrest  initiation  force. 
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Fig.  18.  Force  application  along  the  negative  y  direction  by  one  of 
the  male  subjects  (MSI) --D-nng  force  simulation. 
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Fig.  19.  Force  application  along  the  negative  y  direction  by  one  of 
the  female  subjects  (FS2) --D-ring  force  simulation. 


Fig.  20.  Force  application  on  the  ankle  along  the  positive  y  direction 
by  one  of  the  male  subjects  (MSI). 


Fig.  24.  Force  application  on  the  wrist  by  one  of  the  male  subjects 
(MSI)  when  forearm  is  60°  laterally  rotated. 
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Fig.  27.  Force  application  on  the  elbow  by  one  of  the  male  subjects 
(MSI)  with  fully  extended  elbow  and  arm  orientation  at  <}>  =  0°  5  0  = 


Fig.  28.  Force  application  on  the  wrist  by  one  of  the  male  subjects 
(MSI)  with  fully  extended  elbow  and  arm  orientation  at  4  =  0°  5  0  = 


Fig.  29.  Force  application  on  the  elbow  by  one  of  the  male  subjects 
(MSI)  with  fully  extended  elbow  and  arm  orientation  at  <j>  =  0°  5  9  =  60° 


Fig.  30.  Force  application  on  the  wrist  by  one  of  the  male  subjects 
(MSI)  with  fully  extended  elbow  and  arm  orientation  at  <£  *  0°  §  0  =  120 
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Fig.  31.  Force  application  on  the  elbow  by  one  of  the  male  subjects 
(MSI)  with  fully  extended  elbow  and  arm  orientation  at  <j>  =  0°  5  6  * 


Fig.  32.  Force  application  on  the  wrist  by  one  of  the  male  subjects 
(MSI)  with  fully  extended  elbow  and  arm  orientation  at  $  =  30°  §  6 


Fig.  33.  Force  application  on  the  wrist  by  one  of  the  male  subjects 
(MSI)  with  fully  extended  elbow  and  arm  orientation  at  <J>  =  60°  &  6  =  60°. 


Fig.  34.  Force  application  on  the  wrist  by  one  of  the  male  subjects 
(MSI)  with  fully  extended  elbow  and  arm  orientation  at  ^  *  60°  §  e  =  120°. 
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Fig.  35.  Force  application  on  the  wrist  by  one  of  the  male  subjects 
(MSI)  with  fully  extended  elbow  and  arm  orientation  at  4>  =  90°  5  6  =  30°. 


Fig.  36.  Force  application  on  the  wrist  by  one  of  the  male  subjects 
(MSI)  with  fully  extended  elbow  and  arm  orientation  at  $  =  90®  6  0  «  60°. 
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Fig.  37.  Force  application  on  the  elbow  by  one  of  the  male  subjects 
(MSI)  with  fully  extended  elbow  and  arm  orientation  at  4>  «  90°  5  0  =  60° 


Fig.  38.  Force  application  on  the  wrist  by  one  of  the  male  subjects 
(MSI)  with  fully  extended  elbow  and  arm  orientation  at  $  =  90°  50=  90° 
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Fig.  42.  Force  application  on  the  elbow  by  one  of  the  female  subjects 
(FS2)  with  fully  extended  elbow  and  arm  orientation  at  $  =  0°  5  8  =  3 


Fig.  45.  Force  application  on  the  elbow  by  one  of  the  female  subjects 
(FS2)  with  fully  extended  elbow  and  arm  orientation  at  $  =  0°  §  9  =  90° 


Fig.  46.  Force  application  on  the  wrist  by  one  of  the  female  subjects 
(FS2)  with  fully  extended  elbow  and  arm  orientation  at  <j>  -  0°  5  0  *  120 


Fig.  50.  Force  application  on  the  wrist  by  one  of  the  female  subjects 
(FS2)  with  fully  extended  elbow  and  arm  orientation  at  $  =  90°  59= 


Fig.  51.  Force  application  on  the  wrist  by  one  of  the  female  subjects 
(FS2)  with  fully  extended  elbow  and  arm  orientation  at  $  =  90°  50= 


Fig.  55.  Force  application  on  the  knee  by  one  of  the  male  subjects 
(MSI)  with  fully  extended  knee  and  leg  is  30°  abducted. 


Fig.  56.  Force  application  on  the  ankle  by  one  of  the  male  subjects 
(MSI)  with  fully  extended  knee  and  leg  is  60°  abducted. 


Fig.  57.  Force  application  on  the  ankle  by  one  of  the  female  subjects 
(FS3)  with  knee  at  90°  flexion  and  upper  leg  is  30°  abducted. 


Fig.  58.  Force  application  on  the  knee  by  one  of  the  female  subjects 
(FS3)  with  knee  at  90°  flexion  and  upper  leg  is  30°  abducted. 


Fig.  61  Force  application  on  the  ankle  by  one  of  the  female  subjects 
(FS3)  with  fully  extended  knee  and  leg  is  60°  abducted. 

of  the  force  applicator  which  was  placed  on  the  lower  arm  very  close  to 
the  elbow  joint.  The  same  procedure  was  repeated  for  the  force  applica¬ 
tion  on  the  wrist.  These  forces  were  called  and  F^  for  the  elbow  and 
the  wrist  application  modes,  respectively.  Figures  20  6  21  show  the 
experiment  for  determination  of  the  resistance  of  subjects  against 
separation  of  lower  leg  from  the  seated  configuration. 

The  second  set  of  experiments  was  designed  to  determine  the  active 
muscle  resistance  of  subjects  when  their  extremities  were  dislodged  from 
their  initial  positions.  If  we  use  the  spherical  coordinate  angles  0 
and  <j>  to  define  the  orientation  of  the  upper  arm  or  upper  leg  with  respect 
to  the  torso  we  can  consider  various  combinations  of  8  and  <t>  in  this  set 
of  experiments.  Definition  of  8  and  $  are  given  by  taking  the  upper  arm 
and  the  torso  as  an  example  in  Fig.  13.  Figs.  22-51  show  the  photographs 
of  the  second  set  of  experiments  for  the  upper  extremities  of  both  male 
and  female  subjects.  The  photographs  in  Figs.  52-61  contain  views  for 
the  second  set  of  experiments  conducted  on  the  lower  extremities  of  the 


subjects.  In  all  these  experiments  the  isometric  force  application  the 
subjects  was  approximately  three  seconds  and  the  tests  were  repeated  twice 
in  a  test  session  for  each  subject. 

ACTIVE  MUSCLE  FORCE  AND  MOMENT  RESULTS 
FOR  THE  UPPER  EXTREMITIES 

Active  muscle  force  and  moment  results  for  the  arm  and  shoulder  are 
tabulated  in  Tables  I-XIV  of  Appendix  D.  In  these  tables  the  numerical 
values  for  the  maximum  magnitudes  of  the  active  muscle  forces  of  the  arms 
and  the  corresponding  moments  at  the  shoulder  joints  of  the  subjects  are 
given  for  two  tests.  It  should  be  emphasized  that  the  numerical  values 
givon  in  these  tables  are  the  magnitudes  of  the  force  and  moment  vectors 
having  all  three  components;  although,  depending  upon  the  direction  of  the 
force  application  one  component  of  the  force  vector  predominates  the  other 
two.  It  should  be  also  recognized  that  due  to  static  equilibrium,  the 
magnitude  of  the  force  vector  measured  by  the  force  transducer  actii  g  on 
the  arm  also  represents  the  magnitude  of  the  resistive  force  vector  in  the 
shoulder  complex. 

Coordinate  information  obtained  from  the  sonic  emitters  on  the  upper 
arm  cuff  was  utilized  to  define  the  shoulder  joint  location  in  the 
following  manner.  The  position  of  the  cuff  was  first  defined  with  respect 
to  the  bony  landmarks  such  as  the  medial  and  lateral  epicondyles  of 
humerus;  then,  the  long-bone  axis  of  humerus  was  determined  by  using 
coordinate  data  of  the  sonic  emitters.  The  shoulder  joint  point  about 
which  the  moments  were  calculated  was  assumed  to  be  at  a  prescribed  dis¬ 
tance  on  this  analytically  determined  long-bone  axis.  Note  that  the 
definition  of  shoulder  joint  is  somewhat  arbitrary  since  one  can  choose  an 
infinite  number  of  points  in  the  "shoulder  complex"  space. 

The  term  "shoulder  complex"  refers  to  the  combination  of  the  shoulder 
joint  (frequently  identified  as  the  glenohumeral  joint)  and  the  shoulder 
girdle  which  includes  clavicle  and  scapula  and  their  articulations.  Thus, 
the  usage  of  the  term  shoulder  joint  in  place  of  the  shoulder  complex  is 
somewhat  misleading  if  one  does  not  bear  in  mind  the  various  components 
of  the  region  and  their  articulations.  The  most  striking  feature  of  the 


shoulder  complex  is  the  presence  of  four  independent  articulations  among 
the  bones  of  the  shoulder  complex,  i.e.  the  clavicle,  scapula,  humerus, 
and  the  thorax.  The  sternoclavicular  joint  where  the  clavicle  articulates 
with  respect  to  the  manubrium  of  sternum,  and  the  acromioclavicular  joint 
where  the  clavicle  meets  the  acromion  process  of  the  scapula  are  the  two 
articulations  of  the  clavicle.  The  glenohumeral  joint  is  the  ball  and 
socket  joint  where  the  humerus  articulates  with  the  glenoid  cavity  of  the 
scapula.  The  fourth  articulation  of  the  shoulder  complex  is  the  scapulo- 
thoracic  joint  where  the  scapula  rotates  on  the  thorax.  In  the  true  sense 
the  scapulothoracic  articulation  is  not  a  joint,  but  this  definition  is 
of  some  value  when  describing  the  movements  of  the  scapula  over  the  thorax 
Figure  62  shows  a  typical  resistive  force  application  by  subjects  as 
a  function  of  time.  Since  the  subjects  were  instructed  not  to  suddenly 


Fig.  62  Resistive  muscle  force  vs.  time  curves  for  the  force  Fj  by  the 
male  (MSI,  MS2,  MS3)  and  the  female  (FS1,  FS2,  FS3)  subjects. 
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make  forceful  or  jerky  effort  during  force  application,  tho  maximum  val.es 
of  force  occurred  at  a  later  time  than  it  was  anticipated.  The  force 
values  tabulated  in  tables  are  in  units  of  newtons  and  they  represent  the 
maximum  values  of  curves  similar  to  the  ones  shown  in  Fig.  62.  The  contents 
of  Tables  I  §  II  are  self-exp Ian. tory.  The  tabulated  results  in  Tables 
III  §  IV  were  obtained  from  the  tests  in  which  the  lower  arms  of  the 
subjects  were  dislodged  from  the  seat  armrest  along  the  lateral  direction 
while  the  upper  arm  orientation  was  kept  along  the  torso.  Fig.  63  shows 
results  from  Tables  III  §  IV. 

In  the  second  set  of  experiments,  in  which  both  the  lower  and  upper 
arms  of  the  subjects  were  dislodged  from  the  initial  position,  the 
results  are  quite  extensive  because  of  various  combinations  of  0  and  p 
angles  which  define  the  orientation  of  the  arm,  with  fully  extended  elbow, 
with  respect  to  the  torso.  The  9  angle  refers  to  the  angle  between  the 
z-axis  of  the  torso  and  the  long-bone  axis  of  the  upper  arm  (this  angle 
also  defines  the  shoulder  flexion-extension  in  the  sagittal  plane).  The 
$  angle  refers  to  the  angle  between  the  projection  of  the  long-bone  axis 
of  the  upper  arm  on  the  xy-plane  and  the  x-axis  (the  positive  and 
negative  values  of  this  angle  also  define  the  shoulder  abduction  and 
adduction,  respectively).  The  numerical  results  for  the  second  set  of 
experiments  are  tabulated  in  Tables  V-XIV  and  associated  plots  are  given 
in  Figs.  64-69.  In  Figs.  64-69  various  arm  positions  as  a  function  of  0 
are  considered  in  the  planes  of  $  =  0°  (the  sagittal  plane),  p  =  30°, 
p  =  60°,  ip  =  90°  (the  frontal  plane)  and  p  =  120°.  In  Figs.  64,  65,  5  66 
the  angle  -9  along  with  <j>  =  0°  defines  the  shoulder  extension  positions  in 
the  sagittal  plane.  Figures  64-69  can  be  put  in  three  categories.  In  the 
first  category,  the  maximum  values  of  resistive  muscle  force,  for  the 
external  force  applications  on  the  elbow,  are  plotted  for  various  arm 
positions  for  both  male  and  female  subjects  (Figs.  64  §  67).  In  the 
second  category,  the  maximum  values  of  resistive  muscle  force  by  the  male 
subjects  against  the  external  force  applications  on  the  elbow  and  the  wrist 
at  various  arm  positions  are  plotted  for  various  arm  positions  (Figs.  65, 

6  68).  In  the  third  category,  the  maximum  values  of  resistive  muscle 
force  by  the  female  subjects  against  the  external  force  applications  on 
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Fig.  63  Maximum  values  of  resistive  muscle  force  by  the  male  (MSI,  MS2,  MS3) 
and  the  female  (FS1,  FS2,  FS3)  subjects  at  various  lower  arm  positions  while 
the  upper  arm  is  rotated  by  keeping  its  orientation  along  the  torso. 


UPPER  ARM  POSITION, -6° 

Fig.  64  Maximum  values  of  resistive  muscle  force  by  the  male  (MSI,  MS2,  MS3) 
and  the  female  (FS1,  FS2,  FS3)  subjects  at  various  arm  positions  in  <J>  =  0° 
pla 


77 


0 -SO' 


320  r 

I 

280  j- 

r 

240  j- 

f- 

200  I- 


160 


80 


40 

0 


.J _ _ l _ I _ 

30°  60°  90°  120° 

UPPER  ARM  POSITION,  Q° 


Ll.  CD 

3- 

lS 

2  c: 

50 
=>  — 

Su, 
><  > 


■  cn 


UPPER  ARM  POSITION, 


UPPER  ARM  POSITION,  0° 


Fig.  67  Maximum  values  of  resistive  muscle  force  by  the  male  (MSI,  MS2 
MS3)  and  the  female  (FS1,  FS2,  FS3)  subjects  against  the  external  force 
applications  on  the  elbow  at  various  arm  positions  in  <f>  =  SO”,  60°,  90° 
and  120°  planes. 
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the  elbow  and  the  wrist  at  various  arm  positions  are  plotted  (Figs.  66, 
§  69). 


ACTIVE  MUSCLE  FORCE  AND  MOMENT  RESULTS 
FOR  THE  LOWER  EXTREMITIES 

The  maximum  magnitudes  of  the  active  muscle  forces  of  the  leg  and  the 
corresponding  moments  at  the  hip  joint  are  tabulated  in  Tables  XV-XXVIII 
of  Appendix  E.  Tests  were  conducted  for  various  combinations  of  the 
spherical  coordinates  8  and  $,  describing  the  position  of  the  upper  leg 
in  relation  to  the  torso  with  the  knee  flexed  and  the  knee  locked.  Tests 
were  repeated  twice  for  each  subject.  The  numerical  values  in  the  tables 
are  the  maximum  magnitudes  of  the  active  muscle  forces  of  the  leg  and  the 
corresponding  hip  joint  moments  for  the  two  tests.  Force  applications  by 
the  subjects  were  at  the  knee  and  just  above  the  ankle. 

The  orientation  of  the  upper  leg  with  respect  to  the  torso  was 
determined  in  an  analogous  manner  to  the  upper  extremity  orientation,  as 
discussed  above.  The  9  angle  refers  to  the  angle  between  the  z-axis  of 
the  torso  and  the  femoral  axis  (9  =  90°  for  seated  position;  6=0°  for 
standing  position).  The  angle  <j>  measures  the  projection  of  the  femoral 
axis  on  the  xy  plane  and  the  x-axis. 

The  moments  were  calculated  about  a  point  in  the  hip  which  was  defined 
similarly  to  the  moment  center  of  the  shoulder  complex.  The  femoral  axis 
was  determined  by  obtaining  the  coordinates  of  two  bony  landmarks  of  the 
femur:  the  lateral  margin  of  the  femoral  condyle  and  the  greater 
trochanter.  Using  the  coordinates  from  the  elastic  cuff,  an  axis  was 
defined.  The  moment  center  was  assumed  to  be  a  measured  distance  from 
the  elastic  cuff  on  this  axis.  This  gives  an  arbitrary  point  since  any 
one  of  an  infinite  number  of  points  could  be  chosen. 

Figure  70  shows  the  resistive  force  applications  by  both  male  and 
female  subjects,  for  the  knee  flexed  and  the  knee  locked  configurations. 

On  the  abscissa  of  Fig.  70  the  positions  start  from  an  initial  seated 
position  (9  =  90°,  <j>  =  0°)  and  include  various  leg  positions  corresponding 
to  dislodging  the  leg  along  the  lateral  direction.  Maximum  dislodged 
position  is  taken  as  the  subject's  indicated  limit.  In  addition  to  the 
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seated  configuration,  tests  were  conducted  for  various  positions  of  0, 
with  varying  4>  angles  corresponding  to  dislodging  of  the  leg  along  the 
lateral  direction.  Figures  71  and  72  show  the  maximum  magnitudes  of 
resistive  muscle  force  for  constant  9  values  of  120°  and  40°  respectively. 
Finally,  Fig.  73  depicts  the  maximum  magnitudes  of  resistive  muscle 
force  for  two  laterally  dislodged  positions  of  the  leg  for  9  equal  to  zero. 

RESULTS  ON  THE  FORCED  KINEMATIC  MOTION 
OF  THE  SHOULDER  COMPLEX 

In  addition  to  the  active  force/moment  tests  which  have  been 
discussed  previously,  a  number  of  tests  have  been  performed  which  involve 
the  passive  resistance  of  the  shoulder  complex  to  forced,  drawer-type 
motions.  In  these  tests  the  humerus  is  forced  axially  for  a  variety  of 
upper  arm  orientations  and  the  resistive  force  versus  the  joint  center 
displacement  is  monitored.  Representative  results  for  a  male  subject  are 
presented  in  Figs.  77-80.  However,  before  we  present  the  test  results, 
we  will  first  discuss  the  testing  procedure. 

Although  the  forced  drawer  tests  are  similar  to  the  type  of  testing 
done  previously,  a  number  of  changes  have  been  made  in  the  experimental 
set-up  and  in  the  data  analysis  and  acquisition  schemes.  These  changes 
were  made  for  the  following  reasons: 

1)  To  provide  redundant  data  for  uhe  kinematic  analysis.  (More 
emitters  are  used  to  provide  additional  data  in  cases  where 
certain  emitters  are  blocked  from  the  microphone/ sensor  board.) 

2)  To  improve  the  accuracy  of  the  data  acquired. 

3)  To  increase  the  ease  of  data  acquisition. 

An  overall  discussion  of  the  experimental  procedure  will  serve  to  elaborate 
on  the  effect  of  the  hardware  additions  and  software  revisions  on  these 
tests. 

As  previously  explained,  the  microphone/ sensor  board  axis  system  and 
the  fixed  body  axis  system  were  always  described  relative  to  a  laboratory 
axis  system  by  the  physical  measurement  of  the  angular  orientation  of  the 
restraint  chair  and  the  microphone  board  with  respect  to  the  laboratory 
(inertial)  axis  system.  These  data  were  tedious  to  collect,  cumbersome  to 
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Fig.  73  Maximum  values  of  resistive  muscle  force  by  the  male  (MSI,  MS2, 
MS3)  and  the  female  (FS1,  FS2,  FS3)  subjects  for  the  locked  knee  configu¬ 
rations  and  for  two  lateral  positions  of  the  upper  leg  with  9=0°. 


input  in  the  data  analysis  program,  and  subject  to  a  degree  of  "judgment" 
by  the  experimentor  during  the  actual  measurement. 

In  order  to  increase  both  accuracy  and  the  ease  of  positional  data 
acquisition,  a  relative  axes  locator  device  (RALD)  was  developed.  The 
RALD  is  shown  in  Fig.  74.  The  RALD  uses  four  sonic  emitters  arranged  in 
a  prismatic  fashion  to  describe  an  axis  system  and  coordinate  location  in 
space  relative  to  the  microphone  board  axis  system.  During  any  test  set¬ 
up,  the  RALD  is  aligned  with  the  fixed  body  axis  system  and  the  sonic 
emitter  data  are  used  to  calculate  the  direction  cosine  matrix  which 
describes  the  fixed  body  axis  system  relative  to  the  microphone  board 
axis  system.  Inversely,  the  board  axis  system  can  always  be  determined 
with  respect  to  the  fixed  body  axis  system  from  these  data.  Therefore, 
there  is  no  need  to  measure  any  gross  board  or%  chair  angular  orientations. 
This  approach  is  obviously  easier  and  more  flexible  from  a  data  acquisition 
point  of  view.  In  addition,  the  accuracy  of  the  calculated  data  is 


Fig.  75  Upper  arm  cuff  fitted  with  six  sonic  emitters. 
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increased  since  multiple  physical  length  and  angle  measurements  are 
eliminated  along  with  their  more  numerous  transformation  matrix 
calculations.  The  laboratory  axis  system  is  eliminated  and  the  fixed 
body  system  is  always  established  as  the  "inertial"  axis  system.  Physical 
dimensions,  which  are  input  in  the  data  analysis  program,  transfer  the 
location  of  the  fixed  body  center  from  the  RALD  origin  to  its  actual 
position. 

Another  device  which  was  designed  to  increase  the  data  acquisition 
capabilities  was  a  new  arm  cuff  which  has  six  emitters  instead  of  the 
original  four.  This  was  made  possible  by  the  procurement  of  an  additional 
digitizer  multiplexor  unit  which  expanded  the  previous  capability  to  15 
emitters.  The  new  cuff,  which  is  shown  in  Fig.  75,  eliminates  many 
problematic  situations  in  which  emitters  are  blocked  from  sensors  on  the 
board  due  to  the  kinematics  of  the  test  involved.  Since  only  three 


Fig.  76  Subject  in  position  for  superior-inferior  drawer  tests. 
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Fig.  77  Resistive  force  versus  drawer  displacement 
for  9  =  90°  ,  <{>  =  0°  upper  arm  orientation. 
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Fig.  78  Resistive  force  versus  drawer  displacement 
for  9  =  90'’,  <p  ~  30°  upper  arm  orientation. 


MAGNITUDE  OF  RESISTIVE  FORCE  (N) 


DRAWER  DISPLACEMENT  FROM 
NEUTRAL  POSITION  (cm  ) 


(a) 


DRAWER  DISPLACEMENT  FROM 
NEUTRAL  POSITION  (cm) 

(b) 


DRAWER  DISPLACEMENT  FROM 
NEUTRAL  POSITION  (cm ) 


(c) 


Fig.  79  Resistive  force  versus  drawer  displacement  for  a)  0  =  90" 
4>  *  60%  b)  9  =  90°,  4>  =  90°,  and  c)  9  =  90°,  <j>  =  115°  upper  arm 
orientation. 
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three  emitters  can  be  blocked  from  microphone  "view"  without  losing  needed 
kinematic  data.  Of  course,  the  appropriate  computer  software  revisions 
were  made  to  accommodate  these  hardware  additions. 

For  the  actual  tests,  the  subject  was  secured  in  the  chair  restraint 
system  which  allows  unimpeded  movement  of  the  shoulder  complex.  With  the 
subject's  elbow  flexed  to  90°,  the  force  applicator  was  fixed  to  the  arm 
by  means  of  a  rigid  cuff  extension  so  that  the  force  could  be  applied 
along  the  humeral  axis.  Fig.  76  shows  the  subject  in  position  for  the 
0=0°  drawer  tests.  In  addition,  the  rigid  cuff  placed  the  shoulder  in 
a  position  of  approximately  25°  of  medial  rotation.  This  position  is  near 
the  midpoint  of  the  shoulder  medial-lateral  range  and  is  a  comfortable 
position  for  the  subject.  The  test  results  shown  in  Figs.  77-79  were  for 
drawer  motions  with  humeral  axis  orientations  of  0  =  90°  and  $  =  0°,  30°, 
60°,  90°,  and  115°.  Fig.  80  shows  the  results  of  a  0  =  0°  test.  For  each 
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Fig.  80  Resistive  force  versus  drawer  displacement  for  the  superior- 
inferior  drawer  (0  =  0°). 
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drawer  orientation,  the  subject  was  asked  to  allow  his  shoulder  complex  to 
assume  its  most  comfortable,  "natural"  position.  This  established  the 
neutral  position  for  each  drawer.  The  test  results  are  plots  of  force 
magnitude  along  the  force  applicator  axis  versus  displacement  of  the 
humeral  head  along  this  same  axis.  In  each  case,  the  subject  was 
instructed  to  avoid  muscular  resistance,  and  each  test  was  terminated 
when  the  subject  expressed  discomfort.  The  rate  of  displacement  for  each 
test  was  maintained  between  2-3  centimeters  per  second. 

CONCLUDING  REMARKS 

Based  on  the  numerical  results  presented  in  this  report,  several 
concluding  remarks  can  be  made  on  the  active  muscle  force  response  of  the 
upper  and  lower  extremities  of  the  human  body  to  the  externally  applied 
forces : 

1.  On  the  average  for  the  male  subjects  the  maximum  magnitudes  of 
forces  necessary  to  separate  the  lower  arm,  from  the  initial  seat  armrest 
position,  along  the  mediolateral  direction  (F^)  is  higher  than  that  of 
inferior- superior  direction  (F2);  whereas,  for  the  female  subjects  the 
results  indicate  F2  is  larger  than  F^.  In  fact,  average  to  F2  ratios 
are  found  to  be  1.35  for  males  and  .84  for  females.  In  addition,  the 
average  value  of  the  maximum  magnitudes  of  F^  and  F2  for  male  subjects 
are  2.58  and  1.60,  respectively,  times  those  of  the  female  subjects. 

2.  For  the  D-ring  force  simulation  case  the  average  value  of  the 
maximum  magnitudes  of  the  resistive  forces  for  the  male  subjects  is  about 
2.69  times  that  of  the  female  subjects.  Similar  male  to  female  subject 
comparison  gives  the  ratios  of  2.54  and  2.64  for  the  forces  F_  and  F^, 
respectively,  of  the  face  curtain  force  simulation  experiment. 

3.  In  the  experiments  for  which  the  lower  arm  of  the  subject  was 
dislodged  from  the  seat  armrest  along  the  lateral  direction  while  the 
upper  arm  was  kept  along  the  torso,  the  results  for  both  male  and  female 
subjects  show  only  small  variations  on  the  average  value  of  the  maximum 
magnitudes  of  force  F^,  as  a  function  of  the  laterally  rotated  lower  arm 
position.  The  highest  F^  values  for  all  subjects  occur  for  the  initial 
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lower  arm  position.  The  male  vs.  female  subject  comparison  ratios  for  the 
Fj  forces  range  from  2.34  to  2.56  in  this  case. 

4.  For  the  shoulder  extension  positions  in  the  sagittal  plane,  the 
average  values  of  the  maximum  magnitudes  of  the  resistive  muscle  forces, 
against  the  external  force  applications  on  the  elbow,  steadily  increase 
from  229  N  at  6  *  0  to  327  N  at  the  maximum  extension  position  for  the 
male  subjects.  However,  for  the  female  subjects,  the  same  resistive  force 
starts  with  120  N  at  9  =  0  and  increases  to  160  N  at  6  =  -30°  and  decreases 
to  149  N  at  the  maximum  extension  limit.  In  addition,  for  the  same 
shoulder  extension  positions  in  the  sagittal  plane,  the  ratios  of  the 
average  values  of  the  maximum  magnitudes  of  the  resistive  forces  against 
the  force  application  on  the  elbow  to  that  of  the  wrist  range  from  1.89  to 
2.08  for  the  male  subjects,  and  1.83  to  2.08  for  the  female  subjects. 

5.  The  average  values  of  the  maximum  magnitudes  of  the  resistive 

forces,  against  the  external  force  applications  on  the  elbow,  for  the 
dislodged  positions  of  the  arm  with  a  fixed  value  of  <j>  =  30°,  start  with 

213  N  at  0  =  30°  and  decrease  steadily  to  179  N  at  9  =  120°  for  the  male 

subjects;  for  the  female  subjects,  the  same  quantity  starts  with  115  N  at 
9  =  30°,  exhibits  a  gradual  decrease  to  98  N  at  9  =  90°,  and  subsequently 
increase  to  146  N  at  0  =  120°.  The  ratios  of  the  average  values  of  the 
maximum  magnitudes  of  the  resistive  forces  against  the  external  force 
applications  on  the  elbow  to  that  of  the  wrist  range  from  1.58  to  1.68 
for  the  male  subjects,  and  1.96  to  2.33  for  the  female  subjects. 

6.  The  average  values  of  the  maximum  magnitudes  of  the  resistive 

forces,  against  the  external  force  applications  on  the  elbow,  for  the 
dislodged  positions  of  the  arm  with  a  fixed  value  of  $  =  60°,  start  with 

218  N  at  0  =  30°  and  shows  a  gradual  increase  to  240  N  at  9  =  90°  and 

subsequent  drop  to  182  N  at  9  =  120°  for  the  male  subjects;  for  the  female 
subjects,  the  same  quantity  starts  with  107  N  at  9  =  30°,  decreases 
gradually  to  82  N  at  9  =  90°  and  shows  an  increase  to  106  N  at  9  =  120°. 

The  ratios  of  the  average  values  of  the  maximum  magnitudes  of  the 
resistive  forces  against  the  external  force  applications  on  the  elbow  to 
that  of  the  wrist  range  from  1.71  to  2.12  for  the  male  subjects,  and  1.76 
to  2.39  for  the  female  subjects. 
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7.  The  behavior  of  the  average  values  of  the  maximum  magnitudes  of 
the  resistive  forces,  against  the  external  force  applications  on  the 
elbow,  for  the  dislodged  positions  of  the  arm  with  a  fixed  value  of  <j>  = 

90°  (in  the  frontal  plane)  is  very  similar  for  both  male  and  female 
subjects.  This  quantity  starts  with  231  N  for  the  male  subjects  (101  N 
for  the  female  subjects)  at  0  =  30°  and  after  a  slight  increase  at 

0  =  60°,  it  decreases  gradually  to  169  N  (77  N  for  the  female  subjects)  at 
0  =  120°.  The  corresponding  elbow  to  wrist  resistive  force  ratios  range 
from  1.56  to  1.98  for  the  male  subjects,  and  1.85  to  2.26  for  the  female 
subjects. 

8.  The  behavior  of  the  average  values  of  the  maximum  magnitudes  of 
the  resistive  forces,  against  the  external  force  applications  on  the 
elbow,  for  the  dislodged  positions  of  the  arm  with  a  fixed  value  of  = 
120°  is  very  similar  for  both  male  and  female  subjects.  This  quantity 
starts  with  245  N  for  the  male  subjects  (129  N  for  the  female  subjects)  at 
9  =  30°  and  gradually  decreases  to  210  N  (101  N  for  the  female  subjects) 
at  0  *  120°.  The  corresponding  elbow  to  wrist  resistive  force  ratios 
range  from  1.77  to  2.01  for  the  male  subjects,  and  1.77  to  1.96  for  the 
female  subjects. 

9.  For  the  initial  seated  position  (9  =  90°,  <J>  =  0°),  the  maximum 
magnitude  of  force  necessary  to  laterally  dislodge  the  leg  is  greater  when 
the  knee  is  flexed  than  when  the  knee  is  locked,  greater  for  force  applied 
to  the  knee  than  for  force  applied  to  the  ankle,  and  greater  for  the  male 
subjects  than  for  the  female  subjects.  In  terms  of  ratios,  we  find  that 
the  ratio  of  the  average  values  of  the  maximum  magnitudes  of  the  resistive 
forces  with  the  knee  flexed  to  the  force  with  the  knee  locked  is  1.61  for 
the  male  subjects,  and  1.67  for  the  female  subjects,  when  the  force 
application  is  on  the  ankle.  When  the  force  is  applied  at  the  knee,  the 
required  force  to  dislodge  the  leg  is  1.47  times  the  force  when  applied 

on  the  ankle  for  the  male  subjects,  and  for  the  female  subjects,  this 
ratio  is  1.93,  where,  for  both  male  and  female  subjects,  the  knee  is 
locked.  When  the  knee  is  flexed,  the  ratios  of  force  on  the  knee  to  force 
on  the  ankle  are  1.24  for  males,  and  1.55  for  females. 
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10.  The  male  vs.  female  subject  comparison  ratios  for  the  initial 
seated  position  show  that  to  dislodge  the  leg  the  male  subjects  require 
1.51  times  the  force  needed  for  the  female  subjects  for  force  applied  on 
the  knee  with  the  knee  locked,  and  1.98  for  force  application  on  the 
ankle  with  the  knee  locked  position.  When  the  knee  is  flexed,  the  male 
subjects'  response  is  found  to  be  1.52  times  the  female  subjects' 
response,  for  force  applied  at  the  knee,  and  1.90  times  the  female 
subjects'  response  when  the  force  is  applied  at  the  ankle. 

11.  The  maximum  magnitude  of  force  required  to  dislodge  the  leg  is 
generally  greater  for  the  knee  flexed  than  for  the  knee  locked.  The 
ratios  of  these  forces  tend  to  decrease  as  the  subject's  limit  is 
approached  in  lateral  dislodgement  ($  increasing)  for  any  given  position 
of  e.  When  the  angle  of  dislodgement  (<j>)  is  held  constant  this  force 
ratio  for  the  knee  flexed  to  the  knee  locked  configurations  tends  to 
increase  with  increasing  9.  Representative  ratios  are  cited  here  for 

9  =  120°  with  force  application  on  the  knee,  and  dislodgement  O)  going 
from  0°  to  the  subject's  limit.  For  the  male  subjects  the  range  is  1.11 
to  1.03;  for  the  female  subjects,  the  range  is  1.02  to  0.90.  For  the 
case  of  a  constant  4>  with  force  application  on  the  ankle  and  9  going  from 
0°  to  120°,  we  see  a  range  of  0.75  to  1.37,  for  the  male  subjects,  and  a 
range  of  1.22  to  1.45  for  the  female  subjects. 

12.  An  interesting  observation  'an  also  be  made  by  examining  the 
ratios  of  the  required  dislodging  force  applied  to  the  knee  and  the 
required  dislodging  force  applied  to  the  ankle.  For  the  male  subjects, 
this  ratio  tends  to  increase  as  the  leg  is  dislodged  (<f>  increasing)  for 
all  values  of  9  with  either  the  knee  locked  or  flexed  positions.  This 
range,  for  example,  is  from  1.40  to  1.85,  for  9  =  120°,  and  <J>  going  from 
0°  to  the  subject's  lateral  limit,  with  the  knee  flexed  positions.  On 
the  other  hand,  for  the  female  subjects,  this  ratio  decreases  for  the 
same  position  sequences.  For  example,  with  0  =  120°,  and  $  going  from  0° 
to  the  subject's  limit,  the  range  decreases  from  2.72  to  1.77,  for  the 
knee  locked  case. 

Following  concluding  remarks  are  appropriate  for  the  results  of  the 
forced  kinematic  motion  (drawer  tests)  of  the  shoulder  complex.  These 


preliminary  tests  provided  quantitative  data  on  the  variability  of  the 
"stiffness"  of  the  shoulder  complex  dependent  upon  the  upper  arm 
orientation.  The  shoulder  complex  appears  least  stiff  for  the  9  =  90°, 

4>  =  30°  drawer  for  the  tests  which  were  performed.  This  is  very  close  to 
the  midpoint  of  possible  $  angle  orientations  with  6  =  90°.  In  addition,  as 
we  increase  or  decrease  the  $  orientation  from  $  =  30°  the  stiffness  for  a 
given  displacement  increases.  A  number  of  the  curves  exhibit  a  rapid 
increase  in  resistance  initially  followed  by  a  region  of  less  rapid 
increase.  It  is  believed  that  this  phenomenon  is  an  artifact  of  the 
testing  procedure.  It  is  possibly  due  to  a  reflex  resistance  caused  by 
the  activation  of  stretch  receptors  in  the  muscles  of  the  shoulder  complex. 

The  results  presented  in  this  report  were  obtained  from  three  male 
and  three  female  subjects  who  are  young,  healthy  and  with  no  special 
training  to  increase  muscular  strength.  It  is  reasonable  to  expect  that 
the  values  of  the  resistive  muscle  forces  and  their  corresponding  moments 
at  the  joints  depend  not  only  on  sex  and  anthropometry  of  the  subjects  as 
shown  in  this  research  but  also  on  age,  physical  fitness,  and,  particularly, 
degree  of  training  for  muscular  strength.  We  can  also  point  out  that, 
although  there  are  both  intra-  and  inter-subject  variations  for  the 
maximum  values  of  the  active  muscle  force  and  moment  data  for  the 
extremities,  there  are  some  trends  one  can  establish  for  the  behavior  of 
their  magnitudes. 

It  is  expected  that  incorporation  of  the  active  muscle  force  and 
moment  data  into  the  multi- segmented  mathematical  models  of  the  human  body 
should  improve  the  long-time  response  capabilities  of  these  models  so  that 
they  can  simulate  more  realistically  the  biodyrr.mic  events  which  take 
place  prior  to  flail  injuries.  Improvement  •  mathematical  models 

with  incorporation  of  the  active  muscle  force  an^.  .loment  data  results 
from  the  notion  that  ca'culated  deceleration  profiles  of  the  limbs  with 
respect  to  the  torso  when  the  body  is  ejected  into  the  windblast  will  be 
substantially  different  than  those  calculated  from  a  model  including  only 
the  passive  resistive  force  and  moments  in  major  articulating  joints. 
Recognizing  that  dislocations  of  the  major  articulating  joints  along  with 
associated  soft  and  hard  tissue  injuries  are  closely  related  to  the 


magnitude  of  the  deceleration  forces,  one  can  see  the  importance  of 
incorporating  active  muscle  force  and  moment  data  into  the  multi- segmented 
mathematical  models  of  the  human  body. 
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APPENDIX  A 


TRANSDUCER  CALIBRATION  PROGRAM 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


PROGRAM  CHNCAl 

PURPOSE:  TO  SAMPLE  A-D  CHANNELS  1000  TIMES  AND  CALCULATE  THE 
MEAN*  STANDARD  DEVIATION*  MAXIMA*  MINIMA* 

FOR  USE  IN  CALIBRATION  OF  THE  FORCE  TRANSDUCER 

subroutines:  osuatd 

TALLY 

[SEE  COMMENTS  IN  SUBROUTINES! 


INTEGER  CHAN  *  DAT  A ( 1 000 ) *  ST AT 
REAL  VTG( 1000) *TOTAL( 1 ) 

DO  10  J=32»43 
CHAN=J 

DO  30  1=1*1000 

DATA(I)=0 

30  CONTINUE 

DO  40  1=1*1000 

CALL  OSUATD<CHAN*0»DATA<I)»STAT) 

VTG ( I ) =0 . 00030571578*DATA ( I > 

40  CONTINUE 

CALL  TALLY  ( VTG » VTG  *  TOTAL  *  A  VER  *  SD  *  VMIN  *  VMAX » 1 000  *  1  *  IER ) 
URITE(S*SO) 

IF (IER  .EQ.  1)WRITE<5*90) 

IF(IER  .EQ.  2)WRITE(5* 100) 

WRITE(5*60)  CHAN  *  AVER  *  SD  *  VMIN  *  VMAX  *  I ER 
10  CONTINUE 

50  FORMAT < ' - ' *  1 IX  * ' CHANNEL  NUMBER' *10X* 'MEAN' *10X* 

♦ 'STANDARD  DEVIATION ' , 10X  > ' MINIMA ' * 10X  * ' MAXIMA ' *  1 OX  * 

♦'ERROR  CODE') 

60  FORMAT('0'»T15*I5*T30*F12.?*T50*F16.9*T74*F12.?*T?0>F12.9* 

♦  Tlll*I5/'-'> 

90  FORMAT ( '  '* 'INPUT  IS  NULL.  MINIMA  =  1.E37*  MEAN  =  STANDARD 

♦  DEVIATION  *  0.0  .') 
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100  FORMAT ( '  '» 'INPUT  HAS  ONLY  ONE  NON-ZERO  READING. 
(MINIMA*  STANDARD  DEVIATION  =0.0 
STOP 
END 


MAXIMA  = 


SUBROUTINE  TALLY 
PURPOSE 

CALCULATE  TOTAL.  MEAN.  STANDARD  DEVIATION.  MINIMUM.  MAXIMUM 
FOR  EACH  VARIABLE  IN  A  SET  (OR  A  SUBSET)  OF  OBSERVATIONS 

USAGE 

CALL  T ALL Y ( A . S . TOTAL , AVER . SD  >  VM I N . VMAX . NO . N V . I ER ) 

DESCRIPTION  OF  PARAMETERS 

A  -  OBSERVATION  MATRIX.  NO  BY  NV 
S  -  INPUT  VECTOR  INDICATING  SUBSET  OF  A.  ONLY  THOSE 
OBSERVATIONS  MITH  A  NON-ZERO  S(J)  ARE  CONSIDERED. 
VECTOR  LENGTH  IS  NO. 

TOTAL  -  OUTPUT  VECTOR  OF  TOTALS  OF  EACH  VARIABLE.  VECTOR 
LENGTH  IS  NV. 

AVER  -  OUTPUT  VECTOR  OF  AVERAGES  OF  EACH  VARIABLE.  VECTOR- 
LENGTH  IS  NV. 

SD  -  OUTPUT  VECTOR  OF  STANDARD  DEVIATIONS  OF  EACH 
VARIABLE.  VECTOR  LENGTH  IS  NV. 

VMIN  -  OUTPUT  VECTOR  OF  MINIMA  OF  EACH  VARIABLE.  VECTOR 
LENGTH  IS  NV. 

VMAX  -  OUTPUT  VECTOR  OF  MAXIMA  OF  EACH  VARIABLE.  VECTOR 
LENGTH  IS  NV. 

NO  -  NUMBER  OF  OBSERVATIONS 

NV  -  NUMBER  OF  VARIABLES  FOR  EACH  OBSERVATION 

IER  -  ZERO.  IF  NO  ERROR. 

-  l.IF  S  IS  NULL.  VMIN=1 . E37 . VMAX=-1 . E37 , . SD=AVER=0 

-  2.  IF  S  HAS  ONLY  ONE  NON-ZERO  ELEMENT.  VMIN=VMAX. 
SD=0 . 0 

REMARKS 

NONE 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
NONE 

METHOD 

ALL  OBSERVATIONS  CORRESPONDING  TO  A  NON-ZERO  ELEMENT  IN  S 


99 


ooo  non  non  non  oooooooo 


VECTOR  ARE  ANALYZED  FOR  EACH  VARIABLE  IN  HATRIX  A. 

TOTALS  ARE  ACCUMULATED  AND  MINIMUM  AND  MAXIMUM  VALUES  ARE 
FOUND.  FOLLOWING  THIS.  MEANS  AND  STANDARD  DEVIATIONS  ARE 
CALCULATED.  THE  DIVISOR  FOR  STANDARD  DEVIATION  IS  ONE  LESS 
THAN  THE  NUMBER  OF  OBSERVATIONS  USED. 


SUBROUTINE  TALLY (A»S» TOTAL , AVER > SD . VHIN, VMAX »N0 .NV.IER) 
DIMENSION  A(1).S(1) .TOTAL(l) .AVER(l) .SD(1) .VMIN(l) .VMAX(l) 

CLEAR  OUTPUT  VECTORS  AND  INITIALIZE  VMIN.VMAX 

IER=0 

DO  1  K=1»NV 

T0TAL(K)=0.0 

AVER<K)=0.0 

SD<K)=0.0 

VMIN(K)=1.0E37 

1  VMAX(K)=-1.0E37 

TEST  SUBSET  VECTOR 

SCNT=0.0 
DO  7  J=l»NO 
I J=J-N0 

IF(S( J) )  2.7,2 

2  SCNT=SCNT+1 .0 

CALCULATE  TOTAL,  MINIMA,  MAXIMA 

DO  6  1=1, NV 
IJ=IJ+NO 
X=A(IJ) 

TOTAL  d)=TQTALd)+X 
IF(X-VMINd))  3,4,4 

3  VMIN( I)=X 

4  IF(X-VMAX(I))  6,6,5 

5  VMAX(I)=X 

6  SDd  )=SDd)+X#X 

7  CONTINUE 

CALCULATE  MEANS  AND  STANDARD  DEVIATIONS 

IF  (SCNT)8,8,9 

8  IER=1 

GO  TO  15 

9  DO  10  1=1,  NV 

10  AVER  d ) =TQTAL  d ) /SCNT 
IF  (SCNT-1,0)  13,11,13 


rnmtmimmiimm 
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DISPLACEMENT  FUNCTION 
IN  LIST 


0 

2 

4 

6 

8. 


outputs: 

If  ISTAT  =  1  then  CONVERSION  was  successful  l  DATA  contains  the  value 
If  ISTAT  =  -1  then  CONVERSION  was  not  successful 


Nuaber  of  paraaeters 
Address  of  CHANNEL 
Address  of  GAIN 
Address  of  DATA 
Address  of  ST AT 


notes: 


Four  arguaents  MUST  be  specified!!  If  not>  control  will  return 
imaediatelu  &  no  operation  will  be  performed. 


.PSECT  osuatd»ru,i»gbl>rel»con 


ADCSR 

= 

170000 

r Define 

the 

A-to-D  CSR 

ADBUF 

- 

ADCSR+2 

J Define 

the 

A-to-B  Buffer 

DACSR 

S 

ADCSR+16 

»Def ine 

the 

D/A  CSR 

DAXBUF 

r 

ADCSR+20 

fDefine 

D/A 

X-data  register 

DAYBUF 

= 

ADCSR+22 

» Define 

D/A 

Y-data  register 

DMACSR 

5 

ADCSR+32 

> Define 

DMA 

CSR 

DMAUCR 

= 

ADCSR+34 

! Define 

DMA 

word  count  register 

DMAADR 

= 

ABCSR+36 

» Define 

the 

DMA  current  address 

register 


osuatd:: 


fEntra  point 


MOV  tADCSR  >  R3  > Point  RO  to  A  to  D  CSR 

MOV  R5»R4  JPoint  R4  to  the  PARM  list 
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NOV 

(R4)+»R0 

•  Move  the  ARG  COUNT  into  RO 

TSTB 

RO 

•  Is  the  nueber  of  ARGs  >  0? 

BNI 

RTN 

»N0»  Br. 

CMPB 

R0f#4< 

ils  the  nunber  of  ARGs  =4.? 

BNE 

RTN 

•  NQ»  Br.  (Unable  to  service) 

MOV 

@(R4)+»R1 

•  Move  the  CHANNEL  NUMBER  into  R1 

CMP 

Rl»#177 

•  Is  the  CHANNEL  NUMBER  .LE.  Max 

BHI 

140* 

>N0r  Br.  (Bad  Channel) 

SWAB 

R1 

•Position  the  CHANNEL  NUMBER  in  the  hidh  byte 

MOV 

0(R4H»R2 

•  Set  R4  =  IGAIN 

CMP 

R2»#3 

•  Is  the  IGAIN  .LE.  3? 

BLOS 

160* 

»YESf  Br.  (Good  sain  supplied) 

140*J 

MOV 

♦-1.R0 

•  Set  ISTAT  to  return  =  -1 

BR 

RTNSTA 

•  Br.  to  RETURN  with  status 

160*: 

ASL 

R2 

•Position  the  IGAIN  bits 

ASL 

R2 

»  ... 

BIS 

R2»R1 

} Construct  the  BIT  MASK  in  R1 

MOV 

R1»0R3 

•Set  the  CSR  for  the  convert 

INC 

0R3 

•Start  the  conversion 

200*: 

TSTB 

0R3 

•Is  the  conversion  DONE? 

BPL 

200* 

>N0»  Br.  (Loop  to  recheck) 

MOV 

ADBUF-ADCSR(R3)f0(R4)+ 

•Move  the  converted  value  to  the  IDATA  slot 

HOV 

*1»R0 

•Set  RO  »  1  to  i»Ply  success 

rtnsta: 

MOV 

RO>08, (R5) 

•Set  the  ISTAT 

rtn: 

RTS 

PC 

•Return  to  the  CALLER 

END 


fThat's  all  folks 


INDIRECT  COMMAND  FILE  <TRNCAL> 


PROMPTS  FOR  NAME  OR  TYPE  OF  CALIBRATION  TEST 
PROMPTS  FOR  AMOUNT  OF  APPLIED  LOAD 


INSTALLS  AND  RUNS  <CHNCAL>  WHENEVER 
A  CARRIAGE  RETURN  IS  ENTERED  AFTER  THE 
AMOUNT  OF  APPLIED  LOAD. 

STOP  BY  TYPING  CTRL  Z 


TIM 

.ENABLE  SUBSTITUTION 

.ASKS  XXXXXX  WHAT  IS  THIS  TEST  <AXIAL»HOMENT >ETC. >? 


TIM 

.NEW5  .ASKS  XXXXXX  WHAT  IS  THE  AMOUNT  OF  APPLIED  LOAD? 
TIM 


RUN  CHNCAL 
TIM 

.GOTO  NEW 

/ 


> 
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APPENDIX  B 


DATA  ACQUISITION  PROGRAM 


PROGRAM  BIGQIO 

DIMENSION  BIGBUF ( 224 , 40 ) , TEMP <  28  > » OUTPUT (34 ) 
L0GICAU1  BIGBUF  ,  TRANS  <  224 ) 

L0GICAU1  R£CDAT<224»2) 

L0GICAU1  FNAMEU3) 

DIMENSION  IPARAM<6,2) ,I0SB<2,2>  »ATD0UT<6) 
DIMENSION  IPARM(4),I0ST0P<2) 

DIMENSION  SUM(6,2)»ATDDAT(6,40) 

INTEGER  COLUMN » TEST  < 1 ) »  CUT  , CHECK 

DATA  TEST/-1/CHECK/O/COLUMN/1/NOUNT/O/CUT/0/ 

DATA  I0RAL/'1010/IREC/1/ 

DATA  M0DE/1/LM0DE/2/ 

C 

C  CREATE  i  OPEN  OUTPUT  FILE 

C 

WRITE(5»5) 

READ(5»10)  (FNAME(I)»I=1»13) 

WRITE(5, 15) 

READ  (5,20)  NREC 
MREC=NREC*2 

CALL  ASSIGN  (1,FNAME,13) 

DEFINE  FILE  1  <MR£C,68,U, IREC) 

GET  THE  BUFFER  ADDRESSES 

CALL  GETADR( IPARAM (1,1) jRECDAT ( 1 , 1  > ) 

CALL  GET ADR (IP ARAN ( 1 , 2 >  > RECDAT (1,2)) 
IPARAM(2,1)=224 
IPARAM(2>2)=224 
CALL  GETADR( IPARM(1 ) , TEST ( 1 ) ) 

IPARM(2)=1 
C 

C  QUEUE  THE  FIRST  I/O 

C 

CALL  QIOdORAL , 4 , 10, 10,  IOSB ( 1 , 1 ) ,  IPARAM(  1,1)) 
DO  25  1=1,150 
DO  25  J=32,37 
K=J-31 

CALL  OSUATD (J,0»IDATA,ISTAT) 

DATA* I DAT A*0 .00030571578 
SUM(K,1)=SUM(K,1)+DATA 

25  CONTINUE 

C 
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C  INITIALIZE  THE  NUMBER  OF  RECORDS  TRANSFERRED 

C 

100  NMQOE=MODE 

MODE=LMODE 
LMODE=NMODE 
C 

C  WAIT  FOR  THE  BUFFER  TO  FILL 

C 

CALL  HAITFRUO) 

CALL  QIO ( I ORAL » 4 , 1 0 , 1 0  *  I QSB ( 1 > MODE ) » I PARAM ( 1 » MODE ) ) 

1  DO  90  1=1,224 

BIGBUF < I , COLUMN ) =RECDAT ( I , LMODE ) 

90  CONTINUE 

DO  30  1=1,150 
DO  30  J=32»37 
K=J-31 

CALL  OSUATD< J,0,IDATA, ISTAT) 

DATA=IDATA*0. 00030571578 
SUN(K, MODE ) =SUM <  K » MODE ) +DAT A 

30  CONTINUE 

2  DO  35  1=1,4 

ATDDAT ( I , COLUMN )=SUM(I> LMODE ) 
SUN(I,LMODE>=0.0 

35  CONTINUE 

IF  (CHECK  .EG.  NREOGOTO  1200 

C 

C  INCREMENT  THE  NUMBER  OF  RECORDS 

C 

COLUMN=CGLUMN+ 1 
CHECK=COLUMN- 1 
GOTO  100 

1200  UR I TE< 5, 45) CHECK 

DO  999  K=l, CHECK 

DO  998  1=1,224 

TRANS ( I ) =BIGBUF ( I ,  K  > 

998  CONTINUE 

DECODE  (224, 530, TRANS)  ( OUTPUT (J) , J=l,28> 

DO  997  1=1,6 

ATDQUT  < 1 ) = ATDD AT ( I » K ) 

OUTPUT  ( 28+ 1 )  *A  TDOUT  <1X0. 006666667 

997  CONTINUE 

DO  996  1=1,25,4 

IND=IFIX(OUTPUT ( I ) ) 

INDEX=IND+( IND-I >*3 
TEMP ( INDEX ) =OUTPUT ( I ) 

IF(TEMP< INDEX)  .EQ.  0,0)G0T0  9?2 
TEMP ( INDEX+i ) =OUTPUT ( I + 1 ) 
IF(TEMP(INDEX+1)  .EG.  0.0)  GOTO  992 
TEMP  < INDEX+2 )=OUTPUT  < 1+2) 
IF(TEMP(INDEX+2)  .EQ.  0.0)  GOTO  992 
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TEMP ( INDEX+3) “OUTPUT  < I +3 ) 

IF(TEMP(INDEX*3>  ,E0.  0.0)  GOTO  992 

996  CONTINUE 

DO  995  1=1 » 28 

OUTPUT ( I ) “TEMP ( I ) 

995  CONTINUE 

IF (MOUNT  ,6T.  0)  GOTO  993 
994  WRITE(5*535)  (FNAME(I)»I=1*13> 

WRITE (5*540)  (0UTPUT(I)»I=1»16) 

WRITE<5*545)  (0UTPUT(I)*I=17.28) 

WRITE<5*550)  (OUTPUT (I), 1=29 *34) 

993  WRITE(riREC)  (0UTPUT(I>*I=1*34) 

GOTO  991 

992  WRITE(5*560>  MOUNT 

CUT=CUT+1 

991  DO  899  1=1*34 

OUTPUT ( I >=0.0 

899  CONTINUE 

DO  888  1=1*6 

ATD0UT(I)=0.0 

888  CONTINUE 

DO  887  1=1*224 

TRANS<I)='  ' 

887  CONTINUE 

DO  886  1=1*28 

TEMF‘<  I)=0.0 

886  CONTINUE 

KOUNT=KOUNT+1 
999  CONTINUE 

1500  CLOSE  (UNIT=1) 

URITE(5»555)  ( FN AME  < I ) *  I = 1  *  1 3  >  *  KOUNT-CUT 
CALL  CLREF(IO) 

5  FORMAT( '$'» 'Enter  the  naae  to  be  given  to  the  data  file  CS-1305  ') 

10  FORMAT ( 13A1 ) 

15  FORMAT ( '4' * 'Enter  the  nuaber  of  records  (digitizer  sweeps)  to  b 

$e  allocated  to  the  data  file  CN-53!  ') 

20  FORMAT (15) 

45  FORMAT ( 'O'* 'SUCCESS. '*16*'  SWEEPS  RECORDED.') 

530  F0RMAT(7(1X,F4.0*3F8.2*3X) ) 

535  F0RMAT( 'O'* 'PROCESSED  DATA  FOR  FILE.  '*13A1) 

540  F0RMAT('0'*4(F4.0*3F8.2>) 

545  F0RMAT('0'*3(F4,0*3F8.2>) 

550  F0RMAT('O'*6F17.9) 

555  F0RMAT('O'*'DATA  WRITTEN  TO  DISK.  '*13A1* 'CONTAINS' *15* '  RECORD 

*S.'> 

560  FORMAT ( '0' * 'RECORD  NUMBER!  '*15*'  CONTAINED  ZERO  VALUES  AND 

i  HAS  BEEN  DELETED.') 

2000  STOP 

END 
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{Enable  LOWER  CASE 


OSUATD 

The  function  of  this  routine  is  to  control  a  DATA  TRANSLATION 
A-to-D.  This  routine  is  callable  from  FORTRAN*  This  routine 
operates  in  a  'SINGLE  SHOT'  MODE  only.  In  addition*  this  routine 
operates  in  a  programmed  reouest  mode  (no  interrupt). 


FORTRAN  CALL! 

CALL  OSUATD ( ICHAN* IGAIN. IDATA. ISTAT) 
where 


* 


ICHAN  -  Is  an  INTEGER*2  variable  seecifing  the  CHANNEL 
to  read.  Must  be  .GE.  0  .AND,  ,LE.  127. 

IGAIN  -  Is  an  INTEGER*2  variable  specifing  the  GAIN  to 
to  use. 

if  *  0  means  gain  of  1 
if  -  1  means  gain  of  2 
if  =  2  means  gain  of  4 
if  =  3  means  gain  of  8. 

IDATA  -  Is  an  INTEGERS  variable  for  the  DATA  read 

ISTAT  -  Is  an  INTEGERS  variable  for  a  STATUS  word, 
(see  OUTPUTS  S  NOTES  for  meaning) 


INPUTS! 


R5  Points  to  a  PARAMETER  LIST 

DISPLACEMENT  FUNCTION 
IN  LIST 


0 

2 

4 

6 

8. 


Number  of  parameters 
Address  of  CHANNEL 
Address  of  GAIN 
Address  of  DATA 
Address  of  ST AT 


108 


OUTPUTS 


If  ISTAT  =  1  then  CONVERSION  was  successful  &  DATA  contains  the  value 
If  ISTAT  a  -1  then  CONVERSION  was  not  successful 


notes: 


Four  arguments  MUST  be  specified!!  If  not»  control  will  return 
immediately  &  no  operation  will  be  performed* 


.PSECT  OSUATD  >  RW  > I » GBL » REL  » CON 


ADCSR 

= 

170000 

{Define 

the 

A-to-D  CSR 

ABBUF 

= 

ADCSR+2 

{ Define 

the 

A-to-B  Buffer 

DACSR 

= 

ADCSR+16 

.Define 

the 

D/A  CSR 

DAXBUF 

= 

ADCSR+20 

{Define 

D/A 

X-data  register 

DAYBUF 

= 

ADCSR+22 

{Define 

D/A 

Y-data  register 

DMACSR 

s 

ADCSR+32 

» Define 

DMA 

CSR 

DMA NCR 

s 

ADCSRF34 

{Define 

DMA 

word  count  register 

DMAADR 

s 

ADCSR+36 

{Define 

the 

DMA  current  address 

register 


osuatd:: 


.Entry  point 


MOV 

♦ADCSR.R3 

{Point  RO 

to 

A  to  D  CSR 

MOV 

R5.R4 

{Point  R4 

to 

the  FARM  list 

MOV 

(R4H.R0 

{Move  the  ARG  COUNT  into  RO 

TSTB 

RO 

{Is  the  number  of  ARGs  >  0? 

BMI 

RTN 

{NO.  Br. 

CMPB 

R0»#4. 

.Is  the  number  of  ARGs  *  4.? 

BNE 

RTN 

>N0»  Br,  (Unable  to  service) 

MOV 

8(R4)+»R1 

{Move  the  CHANNEL  NUMBER  into  R1 

CMP 

R1 r  #177 

{Is  the  CHANNEL  NUMBER  ,LE.  Max 

BHI 

140t 

{NO.  Br,  (Bad  Channel) 
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A 


swab 

R1 

{Position  the  CHANNEL  NUMBER  in  the  hia'n  bate 

MOV 

8(R4H.R2 

.Set  R4  =  I GAIN 

CMP 

R2,*3 

ils  the  I6A1N  .LE.  3? 

BLOS 

160$ 

>YES*  Br.  (Good  aain  supplied) 

140$: 

MOV 

#-l»RO 

(Set  ISTAT  to  return  =  -1 

BR 

RTNSTA 

(Br.  to  RETURN  with  status 

160$: 

ASL 

R2 

(Position  the  I6AIN  bits 

ASL 

R2 

?  »  »  » 

BIS 

R2»R1 

i Construct  the  BIT  MASK  in  R1 

MOV 

R1.0R3 

(Set  the  CSR  for  the  convert 

INC 

0R3 

(Start  the  conversion 

200$: 

TSTB 

8R3 

(Is  the  conversion  DONE? 

BPL 

200$ 

(N0>  Br.  (Loop  to  recheck) 

MOV 

ADBUF- ADCSR  ( R3 ) ( 8  <  R4 )  + 

(Move  the  converted  value  to  the  IDATA  slot 

MOV 

*1»R0 

ISet  RO  *  1  to  ii»ply  success 

rtnsta: 

MOV 

ROf0B.(R5) 

{Set  the  ISTAT 

rtn: 

RTS 

PC 

(Return  to  the  CALLER 

.END 

(That's  all  folks 

> 
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APPENDIX  C 


DATA  PROCESSING  PROGRAM 


C 

PROGRAM  TABLE 
C 

C  DECLARE  $  TYPE  VARIABLES}  DIMENSION  ARRAYS!  INITIALIZE  CONSTANTS 
C 

DIMENSION  P0SITK3)  >RECDAT<34)  .CAL(6>6) 

DIMENSION  STORE (40. 10 > 

DIMENSION  X(6) »QUTPUT<3> »F1<6) >61(6) 

DIMENSION  PTAXIS(3>3) 

DIMENSION  PT1PT2(3) >PT2PT3<3) >PT1PT3(3) >PT5PT7(3) 

DIMENSION  PT4PT7<3> »PT6PT7<3> 

LOGICAL* 1  FNAME ( 1 3 ) >  JTNAME ( 9 ) > SNAME  <  25 ) >  MESS  <  80 ) 

L0GICALI1  DAY (9) >H0UR(8) 

REAL  NORMAL i LENG1 
REAL  MOMNTX  >  MOMNT Y  >  MOMNTZ 
INTEGER  ANS> Y>N>BODYSD 
DATA  Y/'Y'/N/'N'/ 

COMMON  /AC/  PT4PT5<3) 

COMMON  /BC/  P0INTK3)  >P0INT2(3) »P0INT3(3) 

COMMON  /DC/  P0INT4(3) >P0INT5(3) >P0INT6(3) .P0INT7(3) 

COMMON  /EC/  PTAPP(3) >FIXCNT<3) »F0RCE(6) 

COMMON  /FC/  SONI CX  >  SON ICY. SON I CZ  >  T I LT  >  BDROTO  >  CHT ILT  >  CHROT 0 » 

*  BODYSD  >  XSUB  >  TRNS (3.3) 

DATA  IREC/1/B0DYSD/1/F0RCE/640 . 0/ 

DATA  MOUNT/ 1/ 

C 

C  PROMPT  FOR  DATA  AND  OUTPUT  INFORMATION 

C 

URITE(S.l) 

505  URITE(5>5) 

READ  <  5 » 1 0 » ERR=505 )  ( JTN AME  < I ) > I = 1 >  9 ) 

510  WRITE(5>15) 

READ  <  5 » 20 » ERR=5 10 )  ( SNAME  < I ) . I * 1  *  25 ) 

515  WRITE<5>25) 

READ  <5.30>ERR=515)  (MESS(I) > I=l>80) 

520  URITE(5»35> 

READ  <5>40>ERR=520)  (FNAME< I >.1=1.13) 

525  MRITE(5»45) 

READ  (5.50.ERR=525)  NREC 
C 

C  LOCATE  THE  ORIGIN  OF  THE  3-D  BOARD  WITH  RESPECT  TO  THE 
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FIXED  BODY  SEGMENT 


C 
C 

530  WRITE<5»55> 

535  WRITE<5.60) 

READ  <  5 . 65 » ERR=535 )  SONICX 
540  WRITE(5.70) 

READ ( 5 • 65 * ERR=540 )  SONICY 
545  WRITE(5»75) 

READ<  5  >  65 . ERR=545 )  SONICZ 
WRITE<5»3) 

550  WRITE<5»4) 

READ (5.65.  ERR=550 )  XSUB 
555  WRITE(5»80> 

READ ( 5 » 65 » ERR=555 )  TILT 
560  WRITE (5. 85) 

READ (5.65.  ERR=560 )  BDROTO 
565  WRITEC5.87) 

READ ( 5  r  65 . ERR=565 )  LENG1 
570  WRITE(5»275) 

READ(5»65»ERR=570>  CHTILT 
575  WRITE (5. 280) 

RE AD  <  5 » 65 » ERR=575 )  CHROTO 
C 

C  CONVERT  THE  INPUT  ANGLES  TO  RADIANS 

C 

OUTROT=CHROTO 
OUTTIL=CHTILT 
TILT  =  TILTtO. 0174532925 
BDROTO  =  BDRQT040. 01 74532925 
CHTILT  =  CHTILT4-0 . 01 74532925 
CHROTO  =  CHRQTQ4-0. 0174532925 
C 

C  CHECK  WHICH  SIDE  OF  THE  BODY  IS  BEING  TESTED 

C 

IF (SONICX  ,GE.  0.0)  B0DYSD=-1 
C 

C  LOCATE.  IDENTIFY  AND  ACCESS  THE  DATA  FILE 

C 

2000  CONTINUE 

CALL  ASSIGN  (1.FNAME.13) 

DEFINE  FILE  1  (NREC.68.U.IREC) 

C 

C  READ  ONE  RECORD 

C 

500  READ  (1 ' IREC»ERR=4000)  <RECDAT(I> .1*1.34) 

C 

C  ASSIGN  DATA  TO  VARIABLES 

C 

DO  499  1*1.3 

P0INT1<I)*RECDAT(I+1> 
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P0INT2<  I)*RECDAT  < 1+5) 
P0INT3<I)*RECDAT<I'.») 
PQINT4<I)=R£CDAT<IH3> 
P0INT5(I)=RECDAT (1+17) 

P0INT4( I )=RECDAT < 1+21 ) 

P0INT7  < I ) *RECDAT  < I +25 ) 

499  CONTINUE 

501  FORCE ( 1 ) *RECDAT ( 33 ) 

FORCE ( 2 ) *RECDAT ( 34 ) 

FORCE ( 3 ) =RECDAT  <  29  > 

FORCE ( 4 ) “RECDAT (31) 

F0RC£<5)*R£CBAT(32) 

FORCE  <  6 ) =RECDAT ( 30 ) 

C 

C  AFTER  INITIAL  RECORD »  GO  ONi  OTHERWISE 

C 

IF(IREC  ,GE.  3) GOTO  605 

FILL  THE  CALIBRATION  MATRIX 

CAL(l»l)*-0. 10616 
CAL(1»2)=0. 00164 
CAL<1*3)=0. 00047 
CAL<1*4)=0. 00094 
C AL  < 1 » 5 ) =-0 • 00030 
CAL(lr6)=-0. 00011 
CAL(2»l)*-0. 00106 
CAL(2»2)=0. 10876 
CAL <2» 3) =-0.00001 
CAL<2,4)=-0. 00021 
CAL<2»5)=-0. 00149 
CAL(2»6)=0. 00008 
CAL(3*1)=0. 00090 
CAL(3»2)=0. 00108 
CAL<3»3)=0. 01532 
CAL(3»4)=0.0 
CAL<3»5)=-0. 00016 
CAL ( 3  >  6 )  =0 . 00023 
CAL<4»1)»0. 00229 
CAL ( 4 » 2 ) =0 . 00028 
CAL<4»3)=0. 00037 
CAL(4»4)=0. 02987 
CAL  <  4  ?  5 )  =-0 . 00009 
CAL(4»6)=0. 00012 
CAL (5.1) =0.00033 
CAL(5»2)*0. 00707 
CAL(5*3)«0. 00045 
CAL<5»4)=-0.00028 
CAL<5*5>*0. 02924 
CAL<5»6)=0. 00003 
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CAL  ( 6 » 1 ) =-0 . 00020 
CAL (6. 2) =-0.00390 
CAL ( 6 » 3)=-0 . 00090 
CAL(6»4)=0.0 
CAL(6»5)=0. 00097 
CAL (6. 6 ) =0.03986 
C 

C  INVERT  THE  CALIBRATION  MATRIX 

C 

CALL  MINV(CAL»6.D.F1.G1> 

C 

C  OUTPUT  HEADER  INFORMATION 

C 

CALL  DATE (DAY) 

CALL  TIME(HOUR) 

WRITE  (5.200) 

URITE(S.IOO)  (JTNAME(I) .1=1.9) 

WRITEC5.205) 

WRI TE ( 5 » 105 )  DAY . HOUR , OUTROT . ( SNAME ( I ) . I = 1 . 25 ) . OUTT IL 
URITE(S.llO)  (FNAME< I) *  1=1.13) .NREC-1 » (MESS( I) .1=1 .80) 
URITE(5»205) 

URITE(5.320) 

URITE(5. 325) 

605  CONTINUE 

600  CALL  FORPT ( PT APP » PTAX IS » 80DYSD ) 

C 

C  MULTIPLY  THE  TRANSDUCER  VALUES  BY  THE  CALIBRATION  MATRIX 

C  TO  GET  THE  FORCES 

C 

CALL  GMPRD(CAL, FORCE. X,6. 6.1) 

DO  492  1=1.6 

FORCE(I)=X(I) 

492  CONTINUE 

C 

C  CALCULATE  THE  POSITION  OF  THE  MOVING  BODY  SEGMENT 

C 

CALL  P0SITN(P0SIT1»PT4PT5> 

C 

C  CALCULATE  THETA  AND  PHI  FOR  THIS  POSITION 

C 

CALL  SPHERE(THETA.PHI) 

RECDAT(1)=THETA 

RECDAT(2)=PHI 

C 

C  CALCULATE  THE  JOINT  CENTER  AND  PRINT  RESULTS 

C 

CALL  FJTCNT (P0SIT1.LEN01»FI XCNT ) 

CALL  MATRIX(FIXCNT .OUTPUT) 

RECDAT ( 16 ) “OUTPUT ( 1 ) 

RECDAT (17) ^OUTPUT ( 2 ) 
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RECOAT  (18) =OUTPUT ( 3 ) 

CALL  RESULT < PTAXIS* TOTFOX , TOTFOY » TOTFOZ > TOTMQX » 

♦TOTMOY » TOTHGZ  *  2 » RECDAT ) 

RECDAT  <  3 ) =T0TF0X*-4 . 44898 
RECDAT ( 4  > =TOTFO Y*-4 . 44898 
RECDAT <  5 ) =T0TF0Z*-4 . 44898 

RECDAT ( 6 ) =SQRT ( RECDAT ( 3 ) M2+RECDAT (4 ) **2+RECDAT ( 5 ) **2 ) 

RECDAT  ( 13 )  =SQRT  <  RECDAT  ( 7 )  «2*RECDAT  ( 9  >  **2+RECDAT  ( 1 1 )  **2 ) 

RECDAT (14) =SQRT ( RECDAT ( 8 ) **2+RECDAT ( 10 ) W2+RECDAT  < 12 ) **2 ) 

RECDAT ( 15 ) =SGRT ( ( RECDAT ( 7 ) +RECDAT  < 8 ) ) **2+ ( RECDAT ( 9 ) +RECDAT ( 10 ) ) 
$**2+ ( RECDAT  (11)  +RECDAT  (1 2 ) )  M2  > 

DO  705  1*1 » 3 

PT1PT2( I >=P0INT1 ( 1 >-P0INT2< I ) 

PT2PT3< I )=P0INT2< I )-P0INT3( I ) 

PT1PT3( I >=P0INT1 < I >-P0INT3( I ) 

PT4PT5< I )=P0INT4 < I ) -P0INT5< I ) 

PT5PT7( I ) =P0INT5 < I ) -PO INT7 ( I > 

PT4PT7(I)=P0INT4< I)-PQINT7( I ) 

PT6PT7( I )=P0INT6( D-P0INT7 ( I ) 

CONTINUE 

STORE ( KOUNT . 4 ) =SQRT ( PT 1PT2  < 1 ) M2+PT1PT2 ( 2 ) **2+PT 1PT2 ( 3 ) **2 ) 
STORE  ( KOUNT ,  5 )  *SQRT  ( PT2PT3  ( 1 )  M2+PT2PT3  ( 2 )  **2(-PT2PT3  ( 3 )  **2 ) 
STORE ( KOUNT  >  6 ) =SQRT ( PT 1 PT3 (1 ) **2+PTlPT3 ( 2 ) **2+PT 1 PT3 ( 3 ) **2 ) 
STORE ( KOUNT » 7 ) =SQRT ( PT4PT5 < 1 ) «2tPT4PT5( 2 ) **2+PT4PT5 ( 3 ) **2 ) 
STORE ( KOUNT  f  8 ) =SQRT ( PT5PT7 ( 1 ) **2+PT5PT7 ( 2 ) **2+PT5PT7  <  3 ) **2 ) 
STORE ( KOUNT » 9 ) =SQRT  <  PT4PT7 ( 1 ) **2*PT4PT7 ( 2 ) **2+PT4PT7 ( 3 ) U2 ) 
STORE ( KOUNT , 10 ) =SQRT (PT6PT7 ( 1 ) **2+PT6PT7 ( 2 ) **2+PT6PT7 ( 3 ) «2 ) 

WRITE  THE  OUTPUT  IN  COMPACT  FORM 

DO  702  1=1 i 3 

STORE ( KOUNT . I ) =RECDAT (1*15) 

CONTINUE 

WRITE(5»330)  (RECDAT(I) r 1=1 >15) 

IF  THERE  ARE  ANY  MORE  RECORDS ,  60  6ET  THEM! 

KOUNT =K0UNT+1 

IF(IREC  .NE.  NREOGOTO  500 
WRITE<5»335) 

DO  703  1=1 » KOUNT- 1 

WRITE(5»350)  (ST0RE(I»J).J=1»10) 

CONTINUE 

FORMAT  STATEMENTS  FOR  PROMPTS  AND  RESULTS 


FORMAT ('l'r' Enter  reauested  LENGTHS  in  CENTIMETRES  end  reoueste 
♦d  ANGLES  in  DEGREES:  '//> 

FORMAT ('♦'» 'Enter  name  of  Joint  tested  CS-93t  ') 

FORMAT (9A1) 


15  FORMAT*'!'* 'Enter  subject  name  or  number  CS-253*.  ') 

20  FORMAT (25A1) 

25  FORMAT* '0' * 'Enter  3  description  of  the  test  CS-803  ') 

30  FORMAT (80A1) 

35  FORMAT*'!' .'Enter  data  file  name  CS-133!  ') 

40  FORMAT* 13A1) 

45  FORMAT*'!'* 'Enter  number  of  records  to  be  read  CN-53!  ') 

50  FORMAT (15) 

55  FORMAT* 'O' » 'Enter  the  coordinates  of  the  fixed  body  segment  ori 
tain*  with  respect  to  the  3-D  board  origin!') 

60  FORMAT* '!'.T15. 'Enter  the  X-COORDINATE  CN-83!  ') 

65  F0RMAT*F8.2> 

70  FORMAT* '!' *T15* 'Enter  the  Y-COORDINATE  CN-83!  ') 

75  FORMAT* '!' *T15* 'Enter  the  Z-COORDINATE  CN-83!  ') 

3  FORMAT* 'O'. 'Enter  distance  from  FIXED  BODY  ORIGIN  to  SONIC  POI 
iNT  COORDINATE  in  CENTIMETRES!  ') 

4  FORMAT* '!'*T15. 'Enter  the  X-DISTANCE  with  respect  to  the 

!  fixed  body  coordinate  system  CN-83!  ') 

80  FORMAT* '!'* 'Enter *  in  degrees*  the  rotation  of  the  3-D  board  fr 
tom  the  vertical  CN-83!  ') 

85  FORMAT* '$' * 'Enter*  in  degrees*  the  rotation  of  the  3-D  board  fr 
tom  the  X-Z  plane  CN-83!  ') 

87  FORMAT* '!'  * 'Enter  the  length  to  the  fixed  Joint  center  CN-83! 
t'> 

100  FORMAT* '0' *T78*9A1* 'JOINT' ) 

105  FORMAT* '0'*T5. 'DATE!  ',9A1*T138* 'CHAIR  POSITION! ' */*T5, 'TIME!  ' 
t»8Al,T143* 'ROTATION! '»F8.2*'  degrees'* ATS. 'SUBJECT  NAME  AND  NUMB 
tER!  ' .25A1.T143* 'TILT!  ',F8.2*'  degrees') 

110  FORMAT*'  ' ,T5* 'DATA  FILE  NAME!  '.13A1./.T5. 'NUMBER  OF  RECORDS! 

t' ,IS, //,T5* 'DESCRIPTION!  '»80A1) 

200  FORMAT* '1' *170* '-')/) 

205  FORMAT* '0' *170* '-')//) 

206  FORMAT*'  '.170*'-')) 

207  FORMAT* 'O'. 170*'.') > 

275  FORMAT* 't' * 'Enter  the  tilt  of  the  chair*  in  degrees*  from  the  v 

tertical  CN-83!  ') 

280  FORMAT* 't'* 'Enter  the  rotation  of  the  chair*  in  degrees  CN-83! 
t  ') 

300  F0RMAT('0'.T30. 'ERROR  ON  ATTEMPT  TO  READ  NEXT  RECORD') 

320  FORMAT* '0'*T3. 'THETA'. T15. 'PHI' *T44* 'FORCES  (newtons) ' *T88* 

t' MOMENTS  (newton  metres)') 

325  FORMAT*'  '.T7.'(deSrees)'*T26*'FX'.T37*'FY'*T48*'FZ'*T55. 

t'MAGNITUDE' *T69* 'X!rxF' ,T79. 'X!pure' *T91, 'YlrxF' ,T102* 
t'Y5pure'*T114,'Z!rxF'.T124,'Z!pure'*T134.'MAG!rxF'* 
tT144. 'MAG {pure' *T158* 'TOTAL ' ) 

330  FORMAT* 'O'* 15(F7»2*4X)) 

335  FORMAT* 'O'. 'FIXED  JOINT  CENTERS! '.T88. 'SPARK  GAP  DISTANCES!'* 

t/T 47 * ' 1  to  2'.8X.'2  to  3'.8X.'1  to  3'.8X*'4  to  5'*8X.'5  to  7'*8X* 
$'4  to  7'.8X.'6  to  7') 

340  FORMAT* '0' */'!'* 'Are  there  other  files  to  be  processed? 
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♦CY/N3J  ') 

345  F0RMAT(A4) 

350  FORMAT < '  ' >10<F9.2*5X) ) 

C 

C  CLOSE  UP  DATA  FILE  i  THAT'S  ALL  FOLKS! 

C 

WRIT£<5»207) 

2001  URITE(5»340) 

READ(5»345)ANS 

IF( ANS  *EQ.  'N')SOTO  5000 

CL0SE(UNIT=1) 

WRIT£<5»35) 

READ<5»40)  (FNAME(I)>I=1*13) 
HRITE(5»45) 

READ(5»50)  NREC 

NREC=NREC+1 

WRIT£<5»25) 

READ(5»30)  (MESS<I)fI=l*80) 

IREC=1 

K0UNT=1 

DO  2222  1=1 » 34 

RECDAT ( I ) =0 . 0 

2222  CONTINUE 

GOTO  2000 

4000  URITE(5»205) 

WRITE(5»300) 

GOTO  2001 

5000  WRITE(5»205) 

CLQSE(UNIT=1) 

STOP 

END 


SUBROUTINE  FORPT ( PT APP >  PT AXIS » BODYSD ) 

C 

C  SUBROUTINE  TO  CALCULATE  THE  POINT  OF  FORCE  APPLICATION 

C  AND  THE  AXIS  SYSTEM  OF  THE  FORCE  APPLICATOR 

C 

COMMON  /BC/  P0INT1<3)»P0INT2<3) >P0INT3(3) 

DIMENSION  PTAPP<3)tN0RMAL<3)fPTlPT2<3>rPT2PT3(3> 
DIMENSION  PTAXIS(3»3)»X(3)»Y(3) 

REAL  NORMAL »NORLEN 
INTEGER  BODYSD 
DO  10  1=1 »3 

PT1PT2< I )=P0INT2(I ) -POINT 1 ( I ) 

PT2PT3< I )=P0INT3( I ) -P0INT2( I ) 

10  CONTINUE 
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CALL  CRSPRD ( PT1PT2 » PT2PT3 . NORMAL ) 

IF(BODYSD  .GT.  0)G0T0  IS 
NORMAL ' 1 ) =- 1 . 04N0RMAL < 1 ) 

NORMAL (2)*-l. OtNORMAL ( 2 ) 

NORMAL  <2)=-l. OtNORMAL ( 3 ) 

15  NORLEN*SQRT ( PT2PT3 ( 1 ) «2+PT2PT3 ( 2 ) W2+PT2PT3 ( 3 ) *%2 ) *0 . 5 

CALL  UN1TVR( NORMAL) 

DO  20  1=1 >3 

NORMAL < I ) “NORMAL ( I ) CNORLEN 
PT2PT3(I)=PT2PT3<I)*0.5+P0INT2(I> 

PT APP ( I ) “NORMAL ( I ) +PT2PT3 ( I ) 

X< I)=P0INT3< I)-PTAPP< I ) 

Y(I)=P0INT2<I)-PTAPP(I) 

20  CONTINUE 

CALL  UNITVR(PT1PT2) 

CALL  UNITVR(X) 

CALL  UNITVR(Y) 

DO  30  1=1.3 

PTAPP  < I ) =PTAPP ( I ) +PT1PT2 (I >*28. 1305 
PTAXIS(1»I)=X(I> 

PTAXIS(2.I)=Y<I) 

PTAXIS(3. I)=PT1PT2< I ) 

30  CONTINUE 

RETURN 
END 


SUBROUTINE  F JTCNT (P0SIT1 i LENG1 >  FIXCNT ) 

CALCULATES  A  FIXED  JOINT  CENTER 

COMMON  /AC/  PT4PT5<3> 

DIMENSION  P0SITK3)  »FIXCNT<3) 

REAL  LENG1 
DO  5  1*1.3 

FIXCNT<I)=-1,0¥PT4PT5<I> 

5  CONTINUE 

CALL  UNITVR(FIXCNT) 

DO  10  1*1.3 

FIXCNT< I )=FIXCNT< I >*LENG1+P0SIT1 (I> 
10  CONTINUE 

RETURN 
END 


SUBROUTINE  P0SITN(P0SIT1.PT4PT5) 
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10 

15 


25 


SUBROUTINE  TO  CALULATE  THE  POSITION  OF  THE  MOVING  BODY  SEGMENT 

COMMON  /DC/  POINT4<3)»POINT5(3).POINT6(3).POINT7<3) 

DIMENSI ON  PT4PT5  <  3  > » PT5PT7  <  3  > » POS I T 1 ( 3  > » NORMAL ( 3 ) 

REAL  NORMAL 
DO  10  1=1.3 

PT4PT5  < I ) =P0 1 NTS  < I ) -PO I NT4 (I ) 

PT5PT7 ( I ) =P0INT7  < I ) -P0INT5 ( I ) 

CONTINUE 

CALL  CRSPRD ( PT4PT5 . PT5PT7 . NORMAL ) 

CALL  UNITVR<NQRMAL) 

DIST=SQRT  <  PT5PT7  <  1  >  **2m5PT7  ( 2 )  M2+PT5PT7  ( 3 )  **2  >  *0 . 5 
DO  25  1=1.3 

NORMAL  ( I )  =NORMAL  (IXDIST 
PT5PT7(I)=P0INT5<mPT5PT7<I)*0.5 
POSIT 1 ( I ) =NORMAL ( I ) +PT5PT7 ( I ) 

CONTINUE 

RETURN 

END 


SUBROUTINE  SPHERE ( THET A » PH I ) 

C 

C  SUBROUTINE  TO  CALCULATE  THE  SPHERICAL  COORDINATES  OF  THE 

C  MOVING  BODY  SEGMENT 

C 

DIMENSION  A(3).B(3) 

COMMON  /AC/  PT4PT5C3) 

DATA  PI/3.141592654/ 

A<1)=PT4PT5<1) 

A(2)=PT4PT5(2) 

A(3)=PT4PT5(3) 

CALL  UNITVR(A) 

CALL  MATRIX ( A. B) 

THETA=SQRT  ( B  <  1 )  M2+B  ( 2 )  M2 )  /B  ( 3 ) 

THETA=ATAN< THETA) 

IF(B(3) )  5.10.15 
5  THETA=PI+THETA 

GOTO  15 

10  THETA=PI*0.5 

15  PHI=ATAN2(B(2) »B<1) ) 

IF (B< 1 )  .NE.  0.0) GOTO  20 
IF<B(2)  .GT.  0,0)PHI=PI*0,5 
GOTO  1000 

IF (B(2)  .LT.  0.0)PHI=PI*-0.5 
GOTO  1000 
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20 

IF(PHI)  25»30,35 

25 

IF(B(2>  .GE.  0.0)PHI=PHPHI 

GOTO  1000 

30 

PHI=0.0 

IF <B( 1 )  .LT.  0.0>PHI=PI 

GOTO  1000 

35 

IF (B(2)  .LT.  0.0)PHI=-1.0*PI+PHI 

1000 

PHI=PHI*180.0/PI 

THETA=THETA*180.0/PI 

RETURN 

END 


SUBROUTINE  MATRIX(X.OUTPUT) 

SUBROUTINE  TO  DETERMINE  THE  ROTATION  AND 
TRANSLATION  OF  THE  3-D  BOARD  TO  THE  FIXED  BODY  SEGMENT 
COORDINATE  SYSTEM 

COMMON  /FC/  SONICX»SONICY»SONICZ>TILT»BDRQTO»CHTILT»CHROTO> 
*  BODYSD  t XSUB . TRNS (3.3) 

DIMENSION  TMAT<6i3) f Tl<3»3) »T2<3»3>  »T3(3»3) »T32(3»3> 
DIMENSION  X(3) rOUTPUT(3) .F5(3) »G5<3) »XX(3) 

INTEGER  BODYSD 

DATA  K0UNT/l/Tl/9*0.0/ 

IFLAG=0 

XXU)=X<1> 

XX(2)*X<2) 

XX(3)=X(3) 

ARGV=X < 1 ) **2+X ( 2 ) **2+X  <  3 ) **2 
IF (ARGO  .LE.  1 ,2)  GOTO  5 
I FLAG*1 

XX  < 1 ) =X ( 1 ) -SON ICX 
XX(2)=X(2)-S0NICY 
XX ( 3 ) SX ( 3) -SONICZ 
IF(KOUNT  .GT »  1)  GOTO  15 
Tl(l»l)«. 9702957 
T1U, 3)*. 2419219 
Tl<2»2)=1.0 
T1<3»1>*-. 2419219 
Tl(3»3)=. 9702957 
DO  10  I»l»4f3 

IF(I.EQ.l)  PHI =CHROTG+ 1.570798 
IF(I.EQ.l)  THETA=CHTILT 
IF<I,EG>.4>  PHI=BDROTO 
IF(I.EQ.4>  THETA=TILT 
THAT  < I > 1 ) =COS  <  PH I ) 
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TMAT(I.2)=SIN<PHI) 

TNAT<I»3)=0.0 

THAT ( 1+1 i 1)=-SIN( PHI) <COS( THETA) 

THAT ( I 4 1 » 2 ) =COS ( PHI ) *CGS ( THETA ) 

TMAT< 141 .3)=SIN( THETA) 

TMAT< 142. 1)=SIN(PHI)*SIN< THETA) 
TMAT<I42.2)=-C0S<PHI)*SIN(THETA) 

THAT  ( 1 4-2 » 3) =COS  ( THET  A ) 

10  CONTINUE 

CALL  RCUT<THAT.4.T2.T3»<S,3.0) 

CALL  GMPRD(T3»T2»T32»3»3»3) 

CALL  GMPRD ( T  32  >  T 1 » TRNS .3.3.3) 

K0UNT=KQUNT41 
15  DO  20  1=1.3 

0UTPUT(I)=XX(1)*TRNS(1»1)4XX<2)*TRNS<2»I)4XX(3)*TRNS(3>I) 
20  CONTINUE 

IF< IFLAG  ,EQ.  0)  GOTO  25 
OUTPUT ( 1 ) =OUTPUT ( 1 ) 4XSUB 
25  RETURN 

END 


SUBROUTINE  RESULT(PTAXIS.TQTFOX»TOTFOr.TQTFOZ»TOTHOX» 
♦TOTMOY , TOTMOZ . I , RECDAT) 

C 

C  ASSOCIATES  FORCE  S  MOMENT  MAGNITUDES  WITH  THE  POINT  OF 

FORCE  APPLICATION  VECTOR.  CALCULATES  THE  MOMENTS  <r  x  F). 
PREPARES  THE  TOTALS  FOR  OUTPUT 

COMHON  /EC/  PTAPP<3) »FIXCNT(3) >F0RCE<6) 

COHMON  /FC/  SONICX » SON ICY .SONICZ. TILT .BDRQTO.CHTILT.CHROTO. 
t  BODYSD . XSUB » TRNS (3.3) 

DIMENSION  Y(3)»PTAXIS<3»3)»XM0(3) » YMQ(3) .ZMD(3) .WORK (6. 3) 
DIMENSION  WORKX ( 3 ) » WORK Y ( 3 ) . UORKZ ( 3 ) 

DIMENSION  RECDAT (34) »CPTAPP<3) .X(3) »Z(3> 

INTEGER  BODYSD 
CPTAPP(l)=PTAPP(l) 

CPTAPP(2)=PTAPP(2) 

CPTAPP(3)=PTAPP(3) 

CALL  MATRIX(CPTAPP.Y) 

DO  10  J=l»3 

CPTAPP( J)=Y( J) 

10  CONTINUE 

CALL  MATRIX (FIXCNT.Y) 

DO  15  J=1 .3 

CPTAPP( J)=CPTAPP< J)-Y< J) 

15  CONTINUE 
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DO  7  J=l»3 

X(J)=FORCE(l)*PTAXIS(l»J> 
Y(J)=F0RCE(2)*PTAXIS(2.J> 
Z(J)=F0RCE(3)*PTAXIS<3»J) 
W0RKX(J)=F0RCE(4)*PTAXIS(1»J) 
W0RKY(J)=F0RCE(5)*PTAXIS(2»J> 
WGRKZ(J)=F0RCE(6>*PTAXIS<3.J> 

7  CONTINUE 

FRCXB=X ( 1 ) +Y ( 1 ) +Z  ( 1 ) 

FRCYB=X(2)+Y(2)+Z(2) 

FRCZB=X(3)+Y(3)+Z(3> 

URKXB-UORKX ( 1 ) +WORKY ( 1 ) 4W0RKZ ( 1 ) 

WRKYB=WORKX ( 2 ) +WORKY ( 2 ) +WORKZ ( 2 ) 

WRKZBsWORKX  ( 3 )  4-WORKY  ( 3 )  4-WORKZ  ( 3  > 

DO  8  J=l»3 

WORK ( 1 . J ) =FRC  XBtTRNS ( 1 .  J ) 

WORK <2 » J ) =FRC YB*  TRNS ( 2 » J  > 

WORK <3 » J)=FRCZB*TRNS (  3  >  J ) 
WORK(4.J)=WRKXB*TRNS(l»J> 
W0RK(5»J)=WRKYB*TRNS<2.J) 

WORK ( 6 . J ) =WRKZB *TF:NS  <  3 .  J  > 

8  CONTINUE 

RECDAT ( 8  >  =  < WORK (4.1) +WORK (5.1) + WORK ( 6 » 1 ) ) *-0 . 1 1 298 
RECDAT  < 10)=<W0RK(4»2)+W0RK(5»2)H)0RK(6.2) )  *-0.11298 
RECDAT(12)=(W0RK<4,3)+W0RK(5,3)+W0RK(6»3)) *-0.11298 
DO  20  J=l»3 

X(J)=W0RK(1.J) 

Y( J)=U0RK(2» J) 

Z( J)=W0RK<3. J) 

20  CONTINUE 

CALL  CRSPRD ( CPT APP » X . XMO ) 

CALL  CRSPRD (CPTAPP.Y.YMO) 

CALL  CRSPRD (CPTAPP.Z.ZMO) 

DO  30  J=l»3 

U0RK<4» J)=W0RK(4» J)+XMQ( J) 
W0RK(5»J)=W0RK(5»JWM0(J> 
U0RK(6»J)=W0RK(6.J)+ZM0(J> 

30  CONTINUE 

TOTFOX=UORK (1.1) +WORK (2.1) +WORK (3.1) 

TOTFOY=WORK (1.2) +WORK (2.2) +WORK (3.2) 

TOTFOZ=WORK (1.3) f  WORK (2.3) +W0RK (3.3) 

TOTMOX=WORK (4.1) IWORK ( 5 . 1 ) +WORK ( 6 . 1 ) 

TOTMOY*WORK (4.2) +WORK ( 5 » 2 ) +WORK ( 6 » 2) 

TOTMOZ=WORK (4.3) +WORK (5.3) +WORK (6.3) 

RECDAT ( 7 ) = <  XMO ( 1 ) + YMO ( 1 ) +ZMQ ( 1 )) *-0 . 0444898 
RECDAT(9)=(XMO(2WMO(2)+ZMO(2)  >*-0.0444898 
RECDAT(U>*<XM0(3>+YH0(3)+ZM0(3> >*-0.0444898 
RETURN 
END 
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SUBROUTINE  GMPRD 


C 

c 
c 
c 
c 

C  PURPOSE 

C  MULTIPLY  TWO  GENERAL  MATRICES  TO  FORM  A  RESULTANT  GENERAL 

C  MATRIX 

C 

C  USAGE 

C  CALL  GMPRD(A»B*RiN»M»L> 

C 

C  DESCRIPTION  OF  PARAMETERS 

C  A  -  NAME  OF  FIRST  INPUT  MATRIX 

C  B  -  NAME  OF  SECOND  INPUT  MATRIX 

C  R  -  NAME  OF  OUTPUT  MATRIX 

C  N  -  NUMBER  OF  ROUS  IN  A 

C  M  -  NUMBER  OF  COLUMNS  IN  A  AND  ROUS  IN  B 

C  L  -  NUMBER  OF  COLUMNS  IN  B 

C 

C  REHARKS 

C  ALL  MATRICES  MUST  BE  STORED  AS  GENERAL  MATRICES 

C  MATRIX  R  CANNOT  BE  IN  THE  SAME  LOCATION  AS  MATRIX  A 

C  MATRIX  R  CANNOT  BE  IN  THE  SAME  LOCATION  AS  MATRIX  B 

C  NUMBER  OF  COLUMNS  OF  MATRIX  A  MUST  BE  EQUAL  TO  NUMBER  OF  ROW 

C  OF  MATRIX  B 

C 

C  SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 

C  NONE 

C 

C  METHOD 

C  THE  M  BY  L  MATRIX  B  IS  PREMULTIPLIED  BY  THE  N  BY  M  MATRIX  A 

C  AND  THE  RESULT  IS  STORED  IN  THE  N  BY  L  MATRIX  R. 

C 

C  . 

C 

SUBROUTINE  GMPRD<AiB*R»N»M»L) 

DIMENSION  A(1)'B<1)»RU) 

C 

IR=0 

IK=-M 

DO  10  K=1»L 
IK=IK+M 
DO  10  J*1»N 
IR=IR+1 
JI=J-N 
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IB=IK 
R(IR)=0 
DO  10  I=1»M 
JI=JI+N 
I8=IB+1 

10  R<IR)=R<IR)+AUI)*B<IB) 
RETURN 
END 


C 

C  . . . . 

C 

C  SUBROUTINE  HINV 

C 

C  PURPOSE 

C  INVERT  A  HATRIX 

C 

C  USAGE 

C  CALL  MINV<A»N»BiL»M) 

C 

C  DESCRIPTION  OF  PARAMETERS 

C  A  -  INPUT  MATRIX.  DESTROYED  IN  COMPUTATION  AND  REPLACED  BY 

C  RESULTANT  INVERSE. 

C  N  -  ORDER  OF  MATRIX  A 

C  D  -  RESULTANT  DETERMINANT 

C  L  -  WORK  VECTOR  OF  LENGTH  N 

C  M  -  WORK  VECTOR  OF  LENGTH  N 

C 

C  REMARKS 

C  HATRIX  A  MUSI  BE  A  GENERAL  MATRIX 

C 

C  SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 

C  NONE 

C 

C  METHOD 

C  THE  STANDARD  GAUSS-JORDAN  METHOD  IS  USED.  THE  DETERMINANT 

C  IS  ALSO  CALCULATED.  A  DETERMINANT  OF  ZERO  INDICATES  THAT 

C  THE  MATRIX  IS  SINGULAR. 

C 

C  . . . 

C 

SUBROUTINE  MINV(A.N.D.L.M) 

DIMENSION  A(l>fL(l) .M(l) 

C 


IF  A  DOUBLE  PRECISION  VERSION  OF  THIS  ROUTINE  IS  DESIRED.  THE 
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C  C  IN  COLUMN  1  SHOULD  BE  REMOVED  FROM  THE  DOUBLE  PRECISION 

C  STATEMENT  WHICH  FOLLOWS. 

C 

C  DOUBLE  PRECISION  ArD?BIGA»HOLD*DABS 
C 

C  THE  C  MUST  ALSO  BE  REMOVED  FROM  DOUBLE  PRECISION  STATEMENTS 

C  APPEARING  IN  OTHER  ROUTINES  USED  IN  CONJUNCTION  WITH  THIS 

C  ROUTINE. 

C 

C  THE  DOUBLE  PRECISION  VERSION  OF  THIS  SUBROUTINE  MUST  ALSO 

C  CONTAIN  DOUBLE  PRECISION  FORTRAN  FUNCTIONS.  ABS  IN  STATEMENT 

C  10  MUST  BE  CHANGED  TO  DABS. 

C 

C  . 

SEARCH  FOR  LARGEST  ELEMENT 

D=1.0 

NK=-N 

DO  80  K=1»N 
NK=NK+N 
L(K)=K 
M(K)=K 
KK=NK+K 
BIGA=A(KK) 

DO  20  J=K » M 
IZ=N*(J-1) 

DO  20  I=K»N 
IJ=IZ+I 

10  IF <  ABS(BISA)-  ABS(A(IJ) ) )  15»20f20 
15  BIGA-A(IJ) 

L(K)=I 
M(K)=J 
20  CONTINUE 

INTERCHANGE  ROWS 

J=L(K) 

IF(J-K)  35f 35i25 
25  KI=K-N 
DO  30  1*1 >N 
KI*KI+N 
HOLD*-A(KI) 

JI=KI-KFJ 
A(KI)=A( JI) 

30  A(JI)  =HOLD 
C 

C  INTERCHANGE  COLUMNS 

C 

35  I»M(K) 
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IF(I-K)  45,45,38 
38  JP=N*(I-1> 

DO  40  J*1»N 
JK=NK+J 
J1=JP+J 
HQLD=-A(JK) 

A< JK)=A< JI) 

40  A< JI)  =H0LD 


DIVIDE  COLUMN  BY  MINUS  PIVOT  <VALUE  OF  PIVOT  ELEMENT  IS 
CONTAINED  IN  BIGA) 


45  IF(BIGA)  48,46.48 

46  D=0*0 
RETURN 

48  DO  55  1=1, N 
IF < I-K)  50,55,50 
50  IK=NK+I 

A(IK)=A(IK)/(-BIGA) 
55  CONTINUE 


REDUCE  MATRIX 

DO  65  1=1, N 

IK=NK+I 

H0LD=A(IK) 

IJ=I-N 
DO  65  J=1,N 
IJ=IJ+N 

IF(I-K)  60,65,60 
60  IF(J-K)  62,65,62 
62  KJ=IJ-I+K 

A(IJ)=HOLD*A(KJ)+A(IJ) 
65  CONTINUE 

DIVIDE  ROU  BY  PIVOT 


KJ=K-N 
DO  75  J=1,N 
KJ=KJ+N 

IF(J-K)  70,75,70 
70  A<KJ)=A(KJ)/BIGA 
75  CONTINUE 

PRODUCT  OF  PIVOTS 

D=D*BIGA 

REPLACE  PIVOT  BY  RECIPROCAL 
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A<KK)=1.0/BIGA 
80  CONTINUE 
C 

C  FINAL  ROW  AND  COLUMN  INTERCHANGE 

C 

K=N 

100  K=(K-1) 

IF(K)  150*150.105 
105  I=L(K) 

IF(I-K)  120*120*108 
108  JQ=N*(K-1) 

JR=N*(I-1) 

DO  110  J=1,N 
JK=JQ+J 
HOLD=A( JK) 

JI=JR+J 

A<JK)=-A(JI) 

110  A<JI)  =HOLD 
120  J=M(K) 

IF(J-K)  100.100,125 
125  KI=K-N 

DO  130  1=1 rN 

KI=KI+N 

HOLD=A(KI) 

JI=KI-KKJ 

A<KI)=-A<JI) 

130  A(JI>  =HOLD 
GO  TO  100 
150  RETURN 
END 


C 

C  . 

C 

C  SUBROUTINE  RCUT 

C 

C  PURPOSE 

C  PARTITION  A  MATRIX  BETWEEN  SPECIFIED  ROWS  TO  FORM  TWO 

C  RESULTANT  MATRICES 

C 

C  USAGE 

C  CALL  RCUT  (A*L*R*S*N*M,MS) 

C 

C  DESCRIPTION  OF  PARAMETERS 

C  A  -  NAME  OF  INPUT  MATRIX 

C  L  -  ROW  OF  A  ABOVE  WHICH  PARTITIONING  TAKES  PLACE 
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C  R 

C  S 

C  N 

C  M 

C  MS 

C 
C 
C 
C 

C  REMARKS 

C  MATRIX  R  CANNOT  BE  IN  SAME  LOCATION  AS  MATRIX  A 

C  MATRIX  S  CANNOT  BE  IN  SAME  LOCATION  AS  MATRIX  A 

C  MATRIX  R  CANNOT  BE  IN  SAME  LOCATION  AS  MATRIX  S 

C  MATRIX  R  AND  MATRIX  S  ARE  ALWAYS  GENERAL  MATRICES 

C 

C  SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 

C  LOC 

C 

C  METHOD 

C  ELEMENTS  OF  MATRIX  A  ABOVE  ROW  L  ARE  MOVED  TO  FORM  MATRIX  R 

OF  l-l  ROWS  AND  M  COLUMNS.  ELEMENTS  OF  MATRIX  A  IN  ROW  L 

AND  BELOW  ARE  MOVED  TO  FORM  MATRIX  S  OF  N-L+l  ROWS  AND  M 
COLUMNS 


-  NAME  OF  MATRIX  TO  BE  FORMED  FROM  UPPER  PORTION  OF  A 

-  NAME  OF  MATRIX  TO  BE  FORMED  FROM  LOWER  PORTION  OF  A 

-  NUMBER  OF  ROWS  IN  A 

-  NUMBER  OF  COLUMNS  IN  A 

-  ONE  DIGIT  NUMBER  FOR  STORAGE  MODE  OF  MATRIX  A 
0  -  GENERAL 

1  -  SYMMETRIC 

2  -  DIAGONAL 


SUBROUTINE  RCUTIA.L.R.S.N.M.MS) 
DIMENSION  A<l)fR<l)fSU) 


IR=0 

IS=0 

DO  70  J=1»M 
DO  70  1=1 »N 

FIND  LOCATION  IN  OUTPUT  MATRIX  AND  SET  TO  ZERO 

IF(I-L)  20.10.10 
10  IS=IS+1 
S(IS)=0.0 
GO  TO  30 
20  IR=IR+1 
R(IR)=0.0 

LOCATE  ELEMENT  FOR  ANY  MATRIX  STORAGE  MODE 
30  CALL  LOC(I. J.IJ.N.M.MS) 

TEST  FOR  ZERO  ELEMENT  IN  DIAGONAL  MATRIX 
IF< I J)  40.70.40 
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C  DETERMINE  WHETHER  ABOVE  OR  BELOW  L 

C 

40  IF(I-L)  60.50*50 
50  S(IS)=A(IJ) 

GO  TO  70 
60  R(IR)=A(IJ) 

70  CONTINUE 
RETURN 
END 
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TABLE  I 


MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  5  MOMENT  VECTORS 
AT  THE  SHOULDER  JOINTS  OF  THREE  MALE  SUBJECTS 
DURING  DISLODGING  OF  THE  ARM  FROM  THE  INDICATED  POSITIONS 


ARM  POSITION 

MAGNITUDE 

SUBJECT 
MS  1 

SUBJECT 

MS  2 

SUBJECT 
MS  3 

POSITION  #1 

SEAT  ARMREST 

FORCE 

(N) 

225.33 

243.62 

270.60 

271.90 

189.11 

201.03 

INITIATION  FORCE 

MOMENT 

105.18 

108.00 

99.88 

F1 

(N-m) 

107.56 

109.70 

105.88 

POSITION  #1 

SEAT  ARMREST 

FORCE 

(N) 

174.92 

181.18 

217.95 

235.39 

111.81 

113.62 

INITIATION  FORCE 

MOMENT 

58.08 

54.79 

37.54 

F2 

(N-m) 

73.19 

65.78 

38.68 

FORCE 

144.12 

147.04 

122.71 

POSITION  #2 
D-RING  FORCE 

(N) 

174.33 

153.83 

155.23 

SIMULATION 

MOMENT 

69.96 

70.79 

69.98 

(N-m) 

81.33 

75.45 

87.46 

POSITION  #3 

FACE  CURTAIN 

FORCE 

(N) 

216.91 

269.93 

22S.98 

244.42 

189.63 

206.04 

FORCE  SIMULATION 

MOMENT 

65.81 

S5.08 

83.55 

F3 

(N-m) 

85.06 

58.41 

89.95 

POSITION  #3 

FACE  CURTAIN 

FORCE 

CN) 

252.41 

277.27 

225.42 

250.85 

195.55 

210.67 

FORCE  SIMULATION 

MOMENT 

86.99 

68.70 

83.21 

F4 

(N-m) 

89.51 

78.64 

89.72 

131 

J 


TABLE  II 


MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  5  MOMENT  VECTORS  AT  THE 
SHOULDER  JOINTS  OF  THREE  FEMALE  SUBJECTS  DURING  DISLODGING  OF 
THE  ARM  FROM  THE  INDICATED  POSITIONS 


ARM  POSITION 

MAGNITUDE 

SUBJECT 

FS  1 

SUBJECT 

FS  2 

SUBJECT 
FS  3 

POSITION  n 

SEAT  ARMREST 

FORCE 

(N) 

79.34 

88.31 

95.05 

111.62 

83.47 

84.80 

INITIATION  FORCE 

c 

MOMENT 

36.12 

38.14 

35.36 

F1 

(N-m) 

37.00 

46.25 

38.06 

POSITION  #1 

SEAT  ARMREST 

FORCE 

00 

86.03 

88.27 

95.74 

98.81 

139.88 

140.05 

INITIATION  FORCE 

c 

MOMENT 

21.58 

30.07 

51.12 

F2 

(N-m) 

22.35 

31.52 

63.81 

FORCE 

56.56 

57.19 

45.54 

POSITION  *2 

D-RING  FORCE 

(N) 

63.56 

59.18 

51.27 

SIMULATION 

MOMENT 

25.89 

28.76 

22.00 

(N-m) 

28.82 

29.14 

24.57 

POSITION  #3 

FACE  CURTAIN 

FORCE 

(N) 

86.86 

92.99 

96.52 

120.79 

67.10 

68.88 

FORCE  SIMULATION 

c 

MOMENT 

30.78 

30.42 

22.24 

r 

(N-m) 

34.95 

33.35 

23.62 

POSITION  #  3 

FACE  CURTAIN 

FORCE 

(N) 

68.52 

73.43 

92.95 

100.71 

87.01 

103.87 

FORCE  SIMULATION 

c 

MOMENT 

30.74 

32.60 

30.74 

F4 

(N-m) 

31. 2S 

35.22 

37.45 

TABLE  III 

MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  §  MOMENT  VECTORS 
AT  THE  SHOULDER  JOINTS  OF  THREE  MALE  SUBJECTS  AT  THE 
DISLODGED  POSITIONS  OF  THE  LOWER  ARM  FROM  THE  SEAT  ARMREST 


LOWER  ARM  POSITION 

MAGNITUDE 

SUBJECT 
MS  1 

SUBJECT 
MS  2 

SUBJECT 
MS  3 

30°  LATERALLY 

ROTATED 

FORCE 

(N) 

190.08 

223.62 

186.70 

205.05 

150.63 

170.97 

F1 

MOMENT 

81.06 

70.25 

74.57 

(N-m) 

89.61 

75.36 

84.05 

60°  LATERALLY 

ROTATED 

FORCE 

(N) 

216.58 

224.26 

182.32 

195.10 

126.24 

130.89 

F1 

MOMENT 

81.10 

65.45 

62.66 

(N-m) 

83.33 

69.20 

62.95 

ROTATED  TO  THE 

LATERAL  LIMIT 

FORCE 

(N) 

201.02 

244.17 

184.65 

207.39 

120.30 

124.43 

F1 

MOMENT 

79.27 

64 . 45 

55.40 

(N-m) 

97.62 

72.95 

58.68 
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TABLE  IV 

MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  §  MOMENT  VECTORS 
AT  THE  SHOULDER  JOINTS  OF  THREE  FEMALE  SUBJECTS  AT  THE 
DISLODGED  POSITIONS  OF  THE  LOWER  ARM  FROM  THE  SEAT  ARMREST 


LOWER  ARM  POSITION  MAGNITUDE  SUBJECT  SUJHE(rT  SU^JECT 

r  o  1  r  o  Z  ro  j 


30°  LATERALLY 

ROTATED 

FORCE 

(N) 

68.06 

70.83 

77.18 

80.63 

68.14 

74.57 

F1 

MOMENT 

27.47 

30.72 

25.57 

( N-m ) 

29.12 

33.88 

31.06 

60°  LATERALLY 

ROTATED 

FORCE 

(N) 

77.61 

80.17 

72.34 

75.20 

69.20 

85.87 

F1 

MOMENT 

30.56 

28.26 

29.06 

(N-m) 

32.98 

29.30 

31.58 

ROTATED  TO  THE 
LATERAL  LIMIT 

FORCE 

(N) 

74.61 

77.89 

63.34 

69.37 

79.15 

85.15 

F1 

MOMENT 

26.96 

23.29 

33.04 

(N-m) 

27.29 

26.95 

39.40 

TABLE  V 
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MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  §  MOMENT  VECTORS 
AT  THE  SHOULDER  JOINTS  OF  THREE  MALE  SUBJECTS 
FOR  THE  INDICATED  POSITIONS  OF  THE  ARM  WITH  <f>  =  0° 


ARM  POSITION 

MAGNITUDE 

SUBJECT 
MS  1 

SUBJECT 
MS  2 

SUBJECT 
MS  3 

FORCE  APPLICATION 

FORCE 

(N) 

195.43 

257.96 

278.34 

304.91 

165.49 

169.84 

UN  Int  CLoUn 

0  =  0° 

MOMENT 

(N-m) 

55.74 

74.80 

77.34 

89.41 

59.80 

60.87 

FORCE  APPLICATION 

ON  THE  WRIST 

0  =  0° 

FORCE 

(N) 

113.93 

131.51 

163.43 

170.81 

70.81 

74.76 

MOMENT 

(N-m) 

59.20 

65.49 

76.15 

81.71 

38.18 

38.99 

FORCE  APPLICATION 

ON  THE  ELBOW 

0  =  -30° 

FORCE 

(N) 

336.28 

348.20 

349.22 

358.60 

220.20 

228.40 

MOMENT 

(N-m) 

109.09 

112.64 

100.82 

103.64 

77.38 

82.46 

FORCE  APPLICATION 

ON  THE  WRIST 

0  =  -30° 

FORCE 

(N) 

187.47 

192.11 

158.38 

162.46 

92.30 

94.25 

MOMENT 

(N-ra) 

86.84 

89.12 

77.80 

79.78 

49.50 

51.71 

FORCE  APPLICATION 

ON  THE  ELBOW 

0  =  -LIMIT 

FORCE 

(N) 

343.56 

369.63 

368.62 

403.54 

234.70 

240.85 

MOMENT 

(N-m) 

113.54 

117.80 

105.08 

117.95 

81.49 

82.10 

FORCE  APPLICATION 

ON  THE  WRIST 

0  =  -LIMIT 

FORCE 

(N) 

175.52 

227.63 

167.15 

174.98 

97.32 

102.26 

MOMENT 

(N-m) 

78.85 

103.15 

80.69 

87.54 

50. /5 
52.86 
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TABLE  VI 


MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  5  MOMENT  VECTORS 
AT  THE  SHOULDER  JOINTS  OF  THREE  FEMALE  SUBJECTS 
FOR  THE  INDICATED  POSITIONS  OF  THE  ARM  WITH  $  =  06 


ARM  POSITION 

MAGNITUDE 

SUBJECT 
FS  1 

SUBJECT 
FS  2 

SUBJECT 
FS  3 

FORCE  APPLICATION 

ON  THE  ELBOW 

6  =  0° 

FORCE 

(N) 

110.41 

127.22 

111.21 

131.74 

117.04 

119.67 

MOMENT 

(N-m) 

32.82 

34.07 

36.84 

43.54 

30.59 

33.87 

FORCE  APPLICATION 

ON  THE  WRIST 

9  =  0° 

FORCE 

(N) 

59.51 

54.22 

58.47 

58.95 

57.35 

60.04 

MOMENT 

(N-m) 

30.65 

28.50 

31.11 

31.46 

28.64 

31.28 

FORCE  APPLICATION 

ON  THE  ELBOW 

9  *  -30* 

FORCE 

(N) 

183.55 

193.77 

112.05 

123.24 

171.31 

175.78 

MOMENT 

(N-m) 

55.32 

56.63 

55.91 

38.26 

50.44 

50.96 

FORCE  APPLICATION 

ON  THE  WRIST 

9  =  -30° 

FORCE 

(N) 

90.12 

97.93 

53.50 

57.25 

90.17 

94.04 

MOMENT 

(N-m) 

40.39 

45.56 

27.30 

28.58 

41.66 

42.15 

FORCE  APPLICATION 

ON  THE  ELBOW 

9  =  -LIMIT 

FORCE 

(N) 

106.57 

111.54 

141.38 

142.33 

194.15 

197.65 

MOMENT 

(N-m) 

34.09 

34.75 

41.49 

43.72 

57.72 

59.42 

FORCE  APPLICATION 

ON  THE  WRIST 

9  =  -LIMIT 

FORCE 

(N) 

75.40 

87.19 

64.17 

66.94 

90.99 

103.36 

MOMENT 

(N-m) 

35.02 

42.43 

31.96 

32.30 

44.22 

54.76 
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TABLE  VII 


MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  S  MOMENT  VECTORS 
AT  THE  SHOULDER  JOINTS  OF  THREE  MALE  SUBJECTS 
FOR  THE  DISLODGED  POSITIONS  OF  THE  ARM  WITH  $  =  30° 


ARM  POSITION 

MAGNITUDE 

SUBJECT 
MS  1 

SUBJECT 
MS  2 

SUBJECT 
MS  3 

FORCE  APPLICATION 
ON  THE  ELBOW 

8  =  30° 

FORCE 

(N) 

239.41 

249.87 

262.51 

284.39 

110.97 

128.48 

MOMENT 

(N-m) 

75.62 

81.66 

59.03 

61.76 

42.89 

49.41 

FORCE  APPLICATION 
ON  THE  WRIST 

0  =  30° 

FORCE 

(N) 

142.44 

166.42 

155.99 

163.92 

88.05 

92.52 

MOMENT 

(N-m) 

69.91 

76.09 

73.59 

80.01 

53.26 

56.95 

FORCE  APPLICATION 
ON  THE  ELBOW 

8  =  60° 

FORCE 

(N) 

211.71 

261.64 

183.68 

189.12 

144.58 

210.57 

MOMENT 

(N-m) 

71.94 

86.92 

51.93 

61.35 

56.58 

80.81 

FORCE  APPLICATION 
ON  THE  WRIST 

6  =  60° 

FORCE 

(N) 

130.69 

132.63 

123.70 

132.28 

96.44 

101.70 

MOMENT 

(N-m) 

65.85 

67.27 

52.56 

54.52 

61.05 

63.55 

FORCE  APPLICATION 
ON  THE  ELBOW 

9  =  90° 

FORCE 

(N) 

175.82 

186.68 

199.97 

208.88 

156.30 

183.83 

MOMENT 

(N-m) 

60.30 

62.51 

53.20 

56.44 

56.34 

67.10 

FORCE  APPLICATION 
ON  THE  WRIST 

0  =  90° 

FORCE 

(N) 

117.19 

155.61 

104.66 

120.30 

79.50 

85.31 

MOMENT 

(N-m) 

59.31 

69.14 

48.47 

56.75 

44.99 

47.91 

FORCE  APPLICATION 
ON  THE  ELBOW 

0  =  120° 

FORCE 

(N) 

209.91 

233.27 

156.46 

159.02 

153.83 

161.54 

MOMENT 

(N-m) 

63.20 

70.80 

46.82 

46.88 

60.77 

63.14 

FORCE  APPLICATION 
ON  THE  WRIST 

0  =  120° 

FORCE 

(N) 

95.74 

127.94 

125.07 

134.72 

81.19 

82.21 

MOMENT 

(N-m) 

57.02 

66.70 

62.68 

64.66 

54.98 

56.21 

TABLE  VIII 


MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  §  MOMENT  VECTORS 
AT  THE  SHOULDER  JOINTS  OF  THREE  FEMALE  SUBJECTS 
FOR  THE  DISLODGED  POSITIONS  OF  THE  ARM  WITH  $  =  30° 


ARM  POSITION 

MAGNITUDE 

SUBJECT 
FS  1 

SUBJECT 
FS  2 

SUBJECT 
FS  3 

FORCE  APPLICATION 

ON  THE  ELBOW 

0  =  30° 

FORCE 
(N) ; 

111.90 

121.68 

126.84 

128.23 

92.03 

107.61 

MOMENT 

(N-m) 

35.90 

39.59 

41.55 

43.16 

25.99 

31.59 

FORCE  APPLICATION 

ON  THE  WRIST 

0  =  30° 

FORCE 

(N) 

44.83 

57.75 

53.40 

56.89 

40.65 

51.59 

MOMENT 

(N-m) 

22.12 

27.48 

28.92 

30.92 

21.41 

26.82 

FORCE  APPLICATION 

ON  THE  ELBOW 

0  =  60° 

FORCE 

(N) 

107.70 

110.93 

124.52 

130.20 

97.55 

105.38 

MOMENT 

(N-m) 

35.08 

35.45 

57.75 

42.19 

30.17 

32.89 

FORCE  APPLICATION 

ON  THE  WRIST 

0  =  60° 

FORCE 

(N) 

61.40 

62.09 

48.58 

53.23 

45.36 

48.47 

MOMENT 

(N-m) 

29.11 

31.52 

24.79 

26.86 

22.33 

25.13 

FORCE  APPLICATION 

ON  THE  ELBOW 

0  =  90° 

FORCE 

(N) 

109.36 

113.13 

83.32 

110.23 

83.65 

84.73 

MOMENT 

(N-m) 

39.09 

40.81 

22.26 

38.31 

23.40 

24.15 

FORCE  APPLICATION 

FORCE 

(N) 

37.74 

48.51 

58.95 

71.67 

41.07 

42.49 

UN  lnc  nmol 

0  >  90° 

MOMENT 

(N-m) 

20.92 

22.54 

27.84 

34.94 

20.21 

21.07 

FORCE  APPLICATION 
n\i  tup  pt 

FORCE 

(N) 

160.39 

173.74 

109.67 

134.43 

97.65 

99.79 

0  =  120° 

MOMENT 

(N-m) 

50.73 

53.14 

31.07 

35.57 

29.80 

30.98 

FORCE  APPLICATION 
ON  THE  WRIST 

0  =  120° 

FORCE 

(N) 

65.96 

73.15 

60.57 

69.78 

53.09 

54.58 

MOMENT 

(N-m) 

36.87 

41.51 

29.00 

35.30 

27.75 

27.97 

TABLE  IX 


MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  §  MOMENT  VECTORS 
AT  THE  SHOULDER  JOINTS  OF  THREE  MALE  SUBJECTS 
FOR  THE  DISLODGED  POSITIONS  OF  THE  ARM  WITH  $  =  60° 


ARM  POSITION 

MAGNITUDE 

SUBJECT 
MS  1 

SUBJECT 
MS  2 

SUBJECT 
MS  3 

FORCE  APPLICATION 

ON  THE  ELBOW 

9  =  30° 

FORCE 

(N) 

231.59 

265.38 

273.04 

288.02 

122.72 

128.31 

MOMENT 

(N-m) 

71.94 

86.33 

72.03 

75.89 

46.68 

47.69 

FORCE  APPLICATION 

ON  THE  WRIST 

9  =  30° 

FORCE 

(N) 

135.39 

142.05 

166.15 

175.66 

70.74 

76.16 

MOMENT 

(N-m) 

67.88 

68.20 

82.25 

83.40 

43.22 

45.11 

FORCE  APPLICATION 

ON  THE  ELBOW 

9  =  60° 

FORCE 

(N) 

221.19 

297.04 

255.48 

294.77 

128.97 

162.52 

MOMENT 

(N-m) 

94.27 

97.53 

58.04 

61.36 

50.44 

63.22 

FORCE  APPLICATION 

ON  THE  WRIST 

9  =  60° 

FORCE 

(N) 

153.05 

162.77 

132.97 

134.47 

61.22 

69.89 

MOMENT 

(N-m) 

79.13 

82.29 

58.35 

61.35 

38.46 

42.74 

FORCE  APPLICATION 

ON  THE  ELBOW 

9  =  90° 

FORCE 

(N) 

279.48 

281.73 

266.66 

289.51 

154.85 

165.20 

MOMENT 

(N-m) 

82.88 

84.16 

57.42 

62.15 

52.06 

54.70 

FORCE  APPLICATION 

ON  THE  WRIST 

9  =  90° 

FORCE 

(N) 

135.38 

154.43 

115.52 

135.58 

68.37 

69.18 

MOMENT 

(N-m) 

65.30 

75.98 

53.44 

60.25 

42.75 

43.69 

FORCE  APPLICATION 

ON  THE  ELBOW 

9  =  120° 

FORCE 

(N) 

191.87 

209.96 

189.28 

204.40 

144.49 

148.61 

MOMENT 

(N-m) 

55.41 

63.19 

53.12 

56.10 

49.67 

57.76 

FORCE  APPLICATION 

ON  THE  WRIST 

FORCE 

(N) 

106.95 

114.70 

113.14 

117.04 

59.89 

63.32 

TABLE  X 


MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  5  MOMENT  VECTORS 
AT  'll iE  SHOULDER  JOINTS  OF  THREE  FEMALE  SUBJECTS 
FOR  THE  DISLODGED  POSITIONS  OF  THE  ARM  WITH  <j>  =  60° 


ARM  POSITION 

MAGNITUDE 

SUBJECT 
FS  1 

SUBJECT 
FS  2 

SUBJECT 
FS  3 

FORCE  APPLICATION 

ON  THE  ELBOW 

0  =  30° 

FORCE 

(N) 

87.00 

107.62 

110.12 

128.26 

99.11 

110.32 

MOMENT 

(N-m) 

25.60 

31.20 

34.29 

40.10 

27.67 

32  26 

FORCE  APPLICATION 

ON  THE  WRIST 

0  =  30° 

FORCE 

(N) 

60.43 

66.36 

53.69 

64.96 

44.48 

46.39 

MOMENT 

(N-m) 

29.03 

31.36 

26.49 

32.64 

20.76 

22.41 

FORCE  APPLICATION 

ON  THE  ELBOW 

0  =  60° 

FORCE 

(N) 

50.67 

62.50 

106.58 

115.55 

89.47 

93.84 

MOMENT 

(N-m) 

17.20 

20.25 

30.76 

33.65 

24.52 

26.18 

FORCE  APPLICATION 

ON  THE  WRIST 

0  =  60° 

FORCE 

(N) 

33.01 

39.81 

56.45 

58.22 

43.46 

46.08 

MOMENT 

(N-m) 

16.63 

19.56 

28.75 

30.79 

23.88 

25.26 

FORCE  APPLICATION 

ON  THE  ELBOW 

0  =  90° 

FORCE 

(N) 

52.51 

60.87 

95.78 

110.53 

85.49 

87.08 

MOMENT 

(N-m) 

13.46  | 

18.46 

22.29 

23.41 

21.49 

22.88 

FORCE  APPLICATION 

ON  THE  WRIST 

0  =  90° 

FORCE 

(N) 

35.38 

38.09 

56.41 

59.69 

44.70 

45.27 

MOMENT 

(N-m) 

16.70 

19.80 

25.02 

26.52 

20.66 

21.70 

FORCE  APPLICATION 
n\t  tub  n  i  d  nw 

FORCE 

(N) 

85.92 

126.16 

116.93 

120.65 

91.87 

97.20 

OTr  1 1 1L  L  LDvll 

0  =  120° 

MOMENT 

(N-m) 

27.66 

38.02 

28.85 

30.01 

23.09 

24.43 

FORCE  APPLICATION 

FORCE 

(N) 

32.22 

34.90 

55.02 

57.79 

38.90 

48.45 

UN  Inc  nfKlol 

0  *  120° 

MOMENT 

(N-m) 

18.07 

19.18 

27.41 

28.60 

21.77 

26.49 

TABLE  XI 


MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  S  MOMENT  VECTORS 
AT  THE  SHOULDER  JOINTS  UF  THREE  MALE  SUBJECTS 
FOR  THE  DISLODGED  POSITIONS  OF  THE  ARM  WITH  <f>  *  90° 


ARM  POSITION 

MAGNITUDE 

SUBJECT 
MS  1 

SUBJECT 
MS  2 

SUBJECT 
MS  3 

FORCE  APPLICATION 
ON  THE  ELBOW 

0  =  30° 

FORCE 

CN) 

242.17 

2S8.91 

271.56 

298.58 

ISO. 16 
162.89 

MOMENT 

(N-m) 

82.09 

86.48 

77.14 

81.24 

56.40 

58.53 

FORCE  APPLICATION 
ON  THE  WRIST 

0  =  30° 

FORCE 

(N) 

128.28 

159.72 

154.37 

159.21 

79.44 

107.75 

MOMENT 

(N-m) 

65.18 

84.07 

71.70 

76.98 

43.34 

58.24 

FORCE  APPLICATION 
ON  THE  ELBOW 

0  =  60° 

FORCE 

(N) 

272.55 

282.05 

241.18 

242.61 

162.34 

189.96 

MOMENT 

(N-m) 

88.52 

91.36 

65.35 

66.48 

56.71 

70.75 

FORCE  APPLICATION 
ON  THE  WRIST 

9  =  60° 

FORCE 

(N) 

115.30 

138.25 

139.65 

145.75 

76.88 

86.38 

MOMENT 

(N-m) 

54.39 

65.97 

65.10 

67.60 

44.46 

46.05 

FORCE  APPLICATION 
ON  THE  ELBOW 

9  =  90° 

FORCE 

(N) 

274.79 

288.15 

239.06 

244.58 

145.39 

156.17 

MOMENT 

(N-m) 

77.34 

81.76 

64.72 

66.24 

53.02 

57.30 

FORCE  APPLICATION 
ON  THE  WRIST 

9  =  90° 

FORCE 

(N) 

171.14 

178.44 

134.31 

138.19 

81.39 

91.84 

MOMENT 

(N-m) 

77,34 

86.27 

61.70 

66.70 

47.96 

54.50 

FORCE  APPLICATION 
ON  THE  ELBOW 

0  =  120° 

FORCE 

(N) 

153.98 

168.87 

226.00 

238.61 

109.79 

113.48 

MOMENT 

(N-m) 

47.93 

62.96 

49.83 

57.34 

37.48 

43.00 

FORCE  APPLICATION 
ON  THE  WRIST 

0  a  120° 

FORCE 

(N) 

127.83 

132.02 

134.38 

145.78 

52.80 

54.11 

MOMENT 

(N-m) 

64.29 

65.25 

63.12 

71.79 

31.18 

32.31 
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MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  §  MOMENT  VECTORS 
AT  THE  SHOULDER  JOINTS  OF  THREE  FEMALE  SUBJECTS 
FOR  THE  DISLODGED  POSITIONS  OF  THE  ARM  WITH  $  =  90° 


ARM  POSITION 

MAGNITUDE 

SUBJECT 
FS  1 

SUBJECT 
FS  2 

SUBJECT 
FS  3 

FORCE  APPLICATION 

ON  THE  ELBOW 

0  =  30° 

FORCE 

(N) 

74.77 

81.77 

110.96 

114.49 

108.56 

115.19 

MOMENT 

(N-m) 

20.47 

24.59 

34.08 

36.15 

30.39 

32.38 

FORCE  APPLICATION 

ON  THE  WRIST 

6  =  30° 

FORCE 

(N) 

37.54 

38.80 

41.82 

45.62 

50.59 

53.53 

MOMENT 

(N-m) 

20.98 

21.74 

17.38 

19  .53 

24.86 

26.48 

FORCE  APPLICATION 

ON  THE  ELBOW 

9  =  60° 

FORCE 

(N) 

82.00 

113.11 

105.37 

134.07 

93.52 

111.81 

MOMENT 
(N  -m) 

21.44 

30.53 

31.32 

36.17 

28.31 

33.00 

FORCE  APPLICATION 

ON  THE  WRIST 

9*6  0° 

FORCE 

(N) 

39.97 

47.76 

57.53 

61.99 

47.47 

52.28 

MOMENT 

(N-m) 

20.02 

22.97 

30.81 

32.90 

24.25 

26.53 

FORCE  APPLICATION 

ON  THE  ELBOW 

0  =  90° 

FORCE 

(N) 

59.69 

70.72 

105.33 

118.81 

82.48 

93.33 

MOMENT 

(N-m) 

14.88 

17.94 

30.12 

32.07 

24.59 

27.68 

FORCE  APPLICATION 

ON  THE  WRIST 

0  =  90° 

FORCE 

(N) 

38.48 

42.89 

49.64 

51.27 

48.48 

56.44 

MOMENT 

(N-m) 

18.58 

19.41 

25.37 

26.17 

23.59 

27.08 

FORCE  APPLICATION 

ON  THE  ELBOW 

0  =  120° 

FORCE 

(N) 

61.91 

68.25 

96.00 

105.48 

62.37 

70.15 

MOMENT 

(N-m) 

30.90 

33.89 

45,99 

51.18 

19.81 

21.36 

FORCE  APPLICATION 

ON  THE  WRIST 

9  =  120° 

FORCE 

(N) 

36.28 

46.82 

38.44 

39.30 

23.74 

26.49 

MOMENT 

(N-m) 

17.45 

20.53 

23.33 

24.24 

13.65 

15.01 

142 


TABLE  XIII 


MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  $  MOMENT  VECTORS 
AT  THE  SHOULDER  JOINTS  OF  THREE  MALE  SUBJECTS 
FOR  THE  DISLODGED  POSITIONS  OF  THE  ARM  WITH  $  =  120° 


ARM  POSITION 

MAGNITUDE 

SUBJECT 
MS  1 

SUBJECT 
MS  2 

SUBJECT 
MS  3 

FORCE  APPLICATION 

ON  THE  ELBOW 

0  =  30° 

FORCE 

(N) 

232.37 

237.60 

311.03 

313.11 

187.99 

189.25 

MOMENT 

(N-m) 

74.79 

79.90 

92.98 

93.13 

68.50 

68.73 

FORCE  APPLICATION 

ON  THE  WRIST 

0  =  30° 

FORCE 

(N) 

136.84 

149.78 

138.64 

163.92 

89.37 

102.06 

MOMENT 

(N-m) 

65.56 

72.25 

67.48 

77.29 

48.52 

55.96 

FORCE  APPLICATION 

ON  THE  ELBOW 

0  =  60° 

FORCE 

(N) 

288.36 

293.48 

289.94 

305.39 

105.18 

129.64 

MOMENT 

(N-m) 

94.11 

95.55 

80.83 

87.81 

36.33 

43.55 

FORCE  APPLICATION 

ON  THE  WRIST 

0  =  60° 

FORCE 

(N) 

176.40 

181.44 

154.86 

171.18 

68.49 

76.49 

MOMENT 

(N-m) 

84.98 

87.58 

77.47 

81.65 

38.80 

43.29 

FORCE  APPLICATION 

ON  THE  ELBOW 

0  *  90° 

FORCE 

(N) 

282.09 

293.72 

237.45 

288.61 

130.68 

138.09 

MOMENT 

(N-m) 

92.61 

95.88 

68.43 

81.46 

50.85 

52.07 

FORCE  APPLICATION 

ON  THE  WRIST 

0  =  90° 

FORCE 

(N) 

157.57 

172.01 

135.00 

155.79 

57.08 

56.33 

MOMENT 

(N-m) 

75.89 

82.02 

65.14 

71.28 

32.07 

31.82 

FORCE  APPLICATION 

ON  THE  ELBOW 

0  =  120° 

FORCE 

(N) 

242.36 

255.76 

248.94 

258.96 

127.24 

130.89 

MOMENT 

(N-m) 

70.87 

84.12 

65.25 

68.23 

49.42 

50.58 

FORCE  APPLICATION 

ON  THE  WRIST 

9  =  120° 

FORCE 

(N) 

131.63 

143.31 

121.65 

123.97 

51.67 

54.42 

MOMENT 

(N-m) 

66.34 

68.84 

53.29 

54.88 

29.18 

30.58 
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TABLE  XIV 


MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  §  MOMENT  VECTORS 
AT  THE  SHOULDER  JOINTS  OF  THREE  FEMALE  SyB^ECTS 
FOR  THE  DISLODGED  POSITIONS  OF  THE  ARM  WITH  $  =  120° 


ARM  POSITION 

MAGNITUDE 

SUBJECT 

FS  1 

SUBJECT 

FS  2 

SUBJECT 

FS  3 

FORCE  APPLICATION 
ON  THE  ELBOW 

0  =  30° 

FORCE 

(N) 

123.64 

137.82 

96.08 

96.34 

159.53 

161.00 

MOMENT 

(N-m) 

38.07 

42.29 

31.20 

31.22 

49.07 

49.85 

FORCE  APPLICATION 
ON  THE  WRIST 

6  =  30° 

FORCE 

(N) 

65.77 

68.25 

56.79 

59.33 

88.96 

90.48 

MOMENT 

(N-m) 

31.  IS 
33.89 

29.73 

31.10 

45.04 

46.93 

FORCE  APPLICATION 
ON  THE  ELBOW 

0  =  60° 

FORCE 

(N) 

137.73 

163.10 

80.46 

87.29 

141.90 

160.51 

MOMENT 

(N-m) 

42.51 

49.77 

25.47 

26.66 

39.06 

44.16 

FORCE  APPLICATION 
ON  THE  WRIST 

0  =  60° 

FORCE 

(N) 

70.91 

71.63 

51.93 

52.81 

87.19 

88.96 

MOMENT 

(N-m) 

35.61 

36.37 

26.91 

27.87 

43.07 

43.47 

FORCE  APPLICATION 
ON  THE  ELBOW 

0  =  90° 

FORCE 

(N) 

141.04 

142.08 

77.37 

89.85 

113.56 

124.98 

MOMENT 

(N-m) 

42.58 

43.02 

23.28 

29.42 

31.23 

35.45 

FORCE  APPLICATION 
ON  THE  WRIST 

0  *  90° 

FORCE 

(N) 

63.86 

67.13 

40.11 

46.64 

58.16 

76.43 

MOMENT 

(N-m) 

31.53 

33.31 

20.62 

22.50 

29.05 

34.39 

FORCE  APPLICATION 
ON  THE  ELBOW 

0  =  120° 

FORCE 

(N) 

110.62 

119.85 

76.51 

80.79 

107.82 

112.13 

MOMENT 

(N-m) 

34.42 

37.13 

29.04 

33.68 

28.87 

29.81 

FORCE  APPLICATION 
ON  THE  WRIST 

0  -  120° 

FORCE 

(N) 

70.67 

77.66 

37.76 

44.07 

54.33 

58.01 

MOMENT 

(N-m) 

32.81 

34.73 

18.30 

21.32 

25.01 
27. 2S 
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TABLE  XV 

MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  AND  MOMENT  VECTORS 
AT  THE  HIP  JOINTS  OF  THREE  MALE  SUBJECTS  FOR  THE  INDICATED 
POSITIONS  OF  THE  LEG  WITH  6=0,  KNEE  LOCKED. 


LEG  POSITION 

MAGNITUDE 

SUBJECT 

MS  1 

SUBJECT 

MS  2 

SUBJECT 

MS  3 

FORCE  APPLIED 

ON  THE  KNEE 
<j>  =  30° 

FORCE 

(N) 

363.60 

360.47 

355.90 

326.85 

269.78 

251.76 

MOMENT 

(N-m) 

171.31 

171.94 

185.96 

174.60 

134.68 

125.10 

FORCE  APPLIED 

ON  THE  ANKLE 
<f>  =  30° 

FORCE 

(N) 

200.20 

182.44 

299.60 

274.15 

72.15 

71.90 

MOMENT 

(N-m) 

157.24 

137.62 

225.64 

208.13 

36.19 

28.31 

FORCE  APPLIED 

ON  THE  KNEE 
<j>  =  LIMIT 

FORCE 

(N) 

395.50 

392.38 

315.84 

305.42 

263.29 

228.96 

MOMENT 

(N-m) 

188.34 

189.84 

158.56 

153.90 

101.65 

84.40 

FORCE  APPLIED 

ON  THE  ANKLE 
<t>  =  LIMIT 

FORCE 

(N) 

233.98 

231.03 

248.35 

225.70 

236.88 

224.02 

MOMENT 

(N-m) 

183.13 

181.84 

183.39 

164.16 

142.20 

128.11 

TABLE  XVI 

MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  AND  MOMENT  VECTORS 
AT  THE  HIP  JOINTS  OF  THREE  FEMALE  SUBJECTS  FOR  THE  INDICATED 
POSITIONS  OF  THE  LEG  WITH  0=0,  KNEE  LOCKED. 


LEG  POSITION 

MAGNITUDE 

SUBJECT 

FS  1 

SUBJECT 

FS  2 

SUBJECT 

FS  3 

FORCE  APPLIED 

ON  THE  KNEE 
=  30° 

FORCE 

(N) 

187.70 

179.53 

210.83 

180.61 

250.16 
244 . 26 

MOMENT 

(N-m) 

84.04 

81.44 

76.28 

64.78 

120.61 

117.65 

FORCE  APPLIED 

ON  THE  ANKLE 
<p  =  30° 

FORCE 

CN) 

87.77 

66.52 

135.85 

131.49 

102.29 

101.66 

MOMENT 

(N-m) 

68.22 

53.18 

44.87 

39.73 

71.85 

66.19 

FORCE  APPLIED 

ON  THE  KNEE 
<j>  =  60° 

FORCE 

(N) 

202.43 

171.53 

255.41 

219.00 

229.59 

227.35 

MOMENT 

(N-m) 

92.92 

76.92 

79.32 

75.77 

112.51 

110.22 

FORCE  APPLIED 

ON  THE  ANKLE 
<p  =  60° 

FORCE 

(N) 

92.83 

82.05 

159.06 

152.79 

139.46 

133.20 

MOMENT 

(N-m) 

69.58 

62.51 

38.13 

30.38 

111.08 

105.28 

TABLE  XVII 


MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  AND  MOMENT  VECTORS 
AT  THE  HIP  JOINTS  OF  THREE  MALE  SUBJECTS  FOR  THE  INDICATED 
POSITIONS  OF  THE  LEG  WITH  6  =  40°,  KNEE  LOCKED. 


LEG  POSITION 

MAGNITUDE 

SUBJECT 

MS  1 

SUBJECT 

MS  2 

SUBJECT 

MS  3 

FORCE  APPLIED 

ON  THE  KNEE 
<j>  =  30° 

FORCE 

(N) 

280.75 

279.54 

354.66 

290.08 

273.60 

259.66 

MOMENT 

(N-m) 

126.14 

112.62 

184.87 

155.78 

169.88 

166.26 

FORCE  APPLIED 

ON  THE  ANKLE 

4  «  30° 

FORCE 

(N) 

221.50 

218.16 

271.14 

258.20 

146.71 

137.21 

MOMENT 

(N-m) 

162.71 

157.10 

204.23 

196.21 

124.76 

92.74 

FORCE  APPLIED 

ON  THE  KNEE 
<f>  =  LIMIT 

FORCE 

(N) 

328.91 

344.12 

297.32 

283.29 

322.74 

314.97 

MOMENT 

(N-m) 

162.98 
164. S8 

147.71 

141.79 

159.79 

157.63 

FORCE  APPLIED 

ON  THE  ANKLE 

4>  =  LIMIT 

FORCE 

(N) 

282.86 

281.36 

204.69 

179.79 

126.13 

122.85 

MOMENT 

(N-m) 

221.30 

221.20 

151.02 

126.05 

85.70 

83.42 
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TABLE  XVIII 

MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  AND  MOMENT  VECTORS 
AT  THE  HIP  JOINTS  OF  THREE  MALE  SUBJECTS  FOR  THE  INDICATED 
POSITIONS  OF  THE  LEG  WITH  0  =  40°,  KNEE  FLEXED. 


LEG  POSITION 

MAGNITUDE 

SUBJECT 

MS  1 

SUBJECT 

MS  2 

SUBJECT 

MS  3 

FORCE  APPLIED 

ON  THE  KNEE 
<j>  =  30° 

FORCE 

(N) 

235.97 

233.59 

261.16 

252.96 

261.98 
249 , 76 

MOMENT 

(N-m) 

115.01 

113.38 

126.40 

123.86 

137.74 

130.77 

FORCE  APPLIED 

ON  THE  ANKLE 
<t>  *  30a 

FORCE 

(N) 

150.60 

150.19 

212.05 

211.80 

113.12 

106.22 

MOMENT 

(N-m) 

93.73 

93.14 

158.23 

196.02 

32.88 

31.09 

FORCE  APPLIED 

ON  THE  KNEE 
=  LIMIT 

FORCE 

(N) 

352.08 

349.02 

344.08 

309.57 

301.40 

287.37 

MOMENT 

(N-m) 

169.38 

168.5? 

165. S6 
149.20 

148.51 

143.03 

FORCE  APPLIED 

FORCE 

(N) 

212.26 

196.85 

188.14 

186.31 

119.06 

117.53 

FORCE  APPLIED 
ON  THE  ANKLE 
<j>  =  LIMIT 


MOMENT 

(N-m) 


117.41 

104.79 


142.24 

146.95 


79.46 

76.34 


TABLE  XIX 


MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  AND  MOMENT  VECTORS 
AT  THE  HIP  JOINTS  OF  THREE  FEMALE  SUBJECTS  FOR  THE  INDICATED 
POSITIONS  OF  THE  LEG  WITH  9  =  40°,  KNEE  LOCKED. 


LEG  POSITION 

MAGNITUDE 

SUBJECT 

FS  1 

SUBJECT 

FS  2 

SUBJECT 

FS  3 

FORCE  APPLIED 

ON  THE  KNEE 
<f>  =  30° 

FORCE 

(N) 

187.70 

182.98 

142.93 

138.58 

230.44 

217.97 

MOMENT 

(N-m) 

84.04 

82.67 

69.34 

67.21 

118.87 

113.14 

FORCE  APPLIED 

ON  THE  ANKLE 
d>  =  30° 

FORCE 

(N) 

90.84 

87.77 

72.98 

63.82 

114.32 

107.73 

MOMENT 

(N-m) 

70.69 

68.22 

55.46 

50.28 

73.25 

60.54 

FORCE  APPLIED 

ON  THE  KNEE 
<f>  =  60° 

FORCE 

(N) 

202.43 

200.16 

178.80 

173.00 

215.65 

211.64 

MOMENT 

(N-m) 

92.92 

91.02 

87.92 

78.67 

102.08 

100.42 

FORCE  APPLIED 

ON  THE  ANKLE 

4  =  60° 

FORCE 

(N) 

92.83 

92.20 

54.67 

54.35 

147.56 

135.67 

MOMENT 

(N-m) 

69.58 

68.88 

42.60 

42.38 

90.81 

8S.58 

FORCE  APPLIED 

ON  THE  KNEE 

4  =  LIMIT 

FORCE 

(N) 

248.66 

237.98 

214.24 

203.18 

240.61 

233.66 

MOMENT 

(N-m) 

111.52 

106.46 

104.75 

97.92 

114.86 
105. S6 

FORCE  APPLIED 

ON  THE  ANKLE 

4  *  LIMIT 

FORCE 

(N) 

133.51 

121.57 

107.64 

83.15 

147.95 

143.17 

MOMENT 

(N-m) 

97.69 

91.37 

85.58 

70.68 

87.78 

81.78 

TABLE  XX 


MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  AND  MOMENT  VECTORS 
AT  THE  HIP  JOINTS  OF  THREE  FEMALE  SUBJECTS  FOR  THE  INDICATED 
POSITIONS  OF  THE  LEG  WITH  9  =  40°,  KNEE  FLEXED. 


LEG  POSITION 

MAGNITUDE 

SUBJECT 

FS  1 

SUBJECT 

FS  2 

SUBJECT 

FS  3 

FORCE  APPLIED 

ON  THE  KNEE 

4>  =  30° 

FORCE 

(N) 

187.86 

168.84 

186.09 

185.33 

216.10 

213.08 

MOMENT 

(N-m) 

85.35 

75.95 

88.35 

88.31 

87.22 

85.78 

FORCE  APPLIED 

ON  THE  ANKLE 
<j>  =  30° 

FORCE 

(N) 

105.71 

102.58 

130.57 

124.07 

100.31 

93.17 

MOMENT 

(N-m) 

69.02 

68.49 

94.35 

81.98 

57.44 

54.95 

FORCE  APPLIED 

ON  THE  KNEE 

4  =  60° 

FORCE 

(N) 

225.16 

186.81 

175.74 

174.25 

228.90 

228.40 

MOMENT 

(N-m) 

103.07 

85.45 

68.47 

62.68 

103.51 

101.71 

FORCE  APPLIED 

ON  THE  ANKLE 

4  =  60° 

FORCE 

(N) 

92.49 

92.22 

119.73 

116.03 

121.16 

110.63 

MOMENT 

(N-m) 

62.72 

58.34 

86.29 

84.27 

72.39 

64.10 

FORCE  APPLIED 

ON  THE  KNEE 

4  =  LIMIT 

FORCE 

(N) 

237.39 

229.24 

147.70 

146.74 

217.49 

214.87 

MOMENT 

(N-m) 

108.24 

104.88 

59.85 

50.65 

97.13 

96.18 

FORCE  APPLIED 

ON  THE  ANKLE 

4  =*  LIMIT 

FORCE 

(N) 

115.53 

112.21 

95.89 

91.63 

146.48 

127.15 

MOMENT 

(N-m) 

79.57 

76.89 

65.11 

65.01 

88.64 

73.77 
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TABLE  XXI 


MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  AND  MOMENT  VECTORS 
AT  THE  HIP  JOINTS  OF  THREE  MALE  SUBJECTS  FOR  THE  INDICATED 
POSITIONS  OF  THE  LEG  WITH  9  =  90°,  KNEE  LOCKED. 


LEG  POSITION 

MAGNITUDE 

SUBJECT 

MS  1 

SUBJECT 

MS  2 

SUBJECT 

MS  3 

FORCE  APPLIED 

ON  THE  KNEE 

4>  =  0 

FORCE 

(N) 

179.58 

171.00 

104.88 

100.85 

156.93 

144.48 

MOMENT 

(N-m) 

76.36 

74.22 

47.17 

45.35 

76.97 

72.06 

FORCE  APPLIED 

ON  THE  ANKLE 

ip  =  0 

FORCE 

(N) 

122.76 

113.25 

108.63 

103.95 

68.70 

65.53 

MOMENT 

(N-m) 

97.03 

89.35 

81.76 

78.51 

65.35 

61.85 

FORCE  APPLIED 

FORCE 

(N) 

258.43 

255.37 

278.49 

259.65 

216.60 

204 . 70 

ON  THE  KNEE 

4)  =  30° 

MOMENT 

(N-m) 

123.37 

121.42 

141.21 

130.89 

100.90 

94.23 

FORCE  APPLIED 

ON  THE  ANKLE 
i  =  30° 

FORCE 

(N) 

206.81 

199.42 

219.08 

212.99 

120.17 

114.86 

MOMENT 

(N-m) 

158.40 

153.42 

164.78 

159.56 

104.78 

100.19 

FORCE  APPLIED 

ON  THE  KNEE 
<p  =  60° 

FORCE 

(N) 

241.66 

209.29 

293.26 

264.13 

262.80 

257.43 

MOMENT 

(N-m) 

111.01 

100.19 

143.70 

130.06 

128.19 

125.37 

FORCE  APPLIED 

ON  THE  ANKLE 
<p  =  60° 

FORCE 

(N) 

193.13 

191.78 

235.41 

172.72 

143.23 

122.51 

MOMENT 

(N-m) 

153.58 

152.41 

138.85 

107.27 

110.64 

98.22 

FORCE  APPLIED 

ON  THE  KNEE 

FORCE 

(N) 

363.81 

340.15 

404.89 

366.06 

240.94 

217.99 

<p  =  LIMIT 

MOMENT 

(N-m) 

151.83 

141.75 

196.25 

178.22 

118.43 

107.12 

FORCE  APPLIED 

ON  THE  ANKLE 
<(>  »  LIMIT 

FORCE 

(N) 

216.89 

202.71 

214.93 

188.56 

97.62 

86.39 

MOMENT 

(N-m) 

152. S6 
142.40 

170.00 

142.70 

82.52 

75.08 
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TABLE  XXII 


MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  AND  MOMENT  VECTORS 
AT  THE  HIP  JOINTS  OF  THREE  MALE  SUBJECTS  FOR  THE  INDICATED 
POSITIONS  OF  THE  LEG  WITH  9  =  90° ,  KNEE  FLEXED. 


LEG  POSITION 

MAGNITUDE 

SUBJECT 

MS  1 

SUBJECT 

MS  2 

SUBJECT 

MS  3 

FORCE  APPLIED 

ON  THE  KNEE 

<J>  =  0 

FORCE 

(N) 

267.13 

265.47 

225.98 

219.46 

137.94 

122.20 

MOMENT 

(N-m) 

126.33 

125.96 

112.07 

110.62 

65.03 

57.13 

FORCE  APPLIED 

ON  THE  ANKLE 

0  =  0 

FORCE 

(N) 

246.51 

235.51 

174.45 

139.92 

73.48 

69.44 

MOMENT 

(N-m) 

149.44 

146.99 

108.59 

92.42 

49.43 

47.89 

FORCE  APPLIED 

ON  THE  KNEE 

0  =  30° 

FORCE 

(N) 

301.30 

287.74 

262.33 

231.63 

182.97 

174.16 

MOMENT 

(N-m) 

141.86 

134.09 

135.88 

125.28 

83.35 

78.26 

FORCE  APPLIED 

ON  THE  ANKLE 

0  =  30° 

FORCE 

(N) 

242.00 

241.29 

194.50 

174.54 

93.01 

80.29 

MOMENT 

(N-m) 

138.08 

137.19 

124.76 

113.05 

55.12 

44.41 

FORCE  APPLIED 

ON  THE  KNEE 

0  =  60° 

FORCE 

(N) 

274.83 

261.89 

326.92 

312.84 

224.62 

219.62 

MOMENT 

(N-m) 

117.67 

112.03 

159.67 

151.90 

122.04 

113.92 

FORCE  APPLIED 

ON  THE  ANKLE 

0  =  60° 

FORCE 

(N) 

341.40 

317.56 

206.43 

196.64 

110.16 

105.27 

MOMENT 

(N-m) 

132.69 

125.01 

115.72 

109.67 

66.65 

63.88 

FORCE  APPLIED 

ON  THE  KNEE 

0  =  LIMIT 

FORCE 

(N) 

269.76 

255.58 

281.39 

276.75 

252.14 

232.98 

MOMENT 

(N-m) 

122.66 

117.02 

139.54 

138.47 

126.32 

111.52 

FORCE  APPLIED 

FORCE 

(N) 

212.60 

204.01 

227.62 

220.20 

86.18 

85.76 

FORCE  APPLIED 
ON  THE  ANKLE 
0  =  LIMIT 


MOMENT 

(N-m) 


108.11 

102.94 


130.84 

124.47 


56. 4S 
55.39 


TABLE  XXIII 


MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  AND  MOMENT  VECTORS 
AT  THE  HIP  JOINTS  OF  THREE  FEMALE  SUBJECTS  FOR  THE  INDICATED 
POSITIONS  OF  THE  LEG  WITH  9  =  90°,  KNEE  LOCKED. 


LEG  POSITION 

MAGNITUDE 

SUBJECT 

FS  1 

SUBJECT 

FS  2 

SUBJECT 

FS  3 

FORCE  APPLIED 

ON  THE  KNEE 
<fr  =  0 

FORCE 

(N) 

86.78 

86.70 

118.69 

113.95 

81.30 

80.98 

MOMENT 

(N-m) 

37.36 

37.23 

55.71 

54.06 

35.29 

34.82 

FORCE  APPLIED 

ON  THE  ANKLE 
=  0 

FORCE 

CN) 

53.37 

51.38 

49.31 

48.59 

48.04 

44.04 

MOMENT 

(N-m) 

42.28 

40.59 

40.02 

39.26 

38.35 

35.49 

FORCE  APPLIED 

ON  THE  KNEE 
<f>  =  30° 

FORCE 

(N) 

158.80 

149.29 

173.84 

172.29 

182.28 

181.51 

MOMENT 

(N-m) 

69.35 

64.69 

82.32 

81.42 

83.30 

79.89 

FORCE  APPLIED 

ON  THE  ANKLE 
<)>  =  30° 

FORCE 

(N) 

84.34 

68.08 

78.35 

72.15 

103.24 

97.51 

MOMENT 

(N-m) 

72.76 

59.13 

67.72 

62.19 

83.75 

79.66 

FORCE  APPLIED 

ON  THE  KNEE 

4>  =  60° 

FORCE 

(N) 

157.12 

153.24 

143.48 

132.43 

251.60 

241.19 

MOMENT 

(N-m) 

63.69 

61.91 

69.67 

67.03 

114.94 

109.64 

FORCE  APPLIED 

ON  THE  ANKLE 

4  =  60° 

FORCE 

(N) 

70.02 

68.02 

89.32 

83.98 

118.93 

102.39 

MOMENT 

(N-m) 

53.57 

52.49 

74.39 

70.45 

99.27 

84.42 

FORCE  APPLIED 

ON  THE  KNEE 
<f>  =  LIMIT 

FORCE 

(N) 

194.25 

187.30 

156.35 

153.27 

298.85 

278.59 

MOMENT 

(N-m) 

88.25 

84.68 

71.83 

71.46 

143.99 

131.10 

FORCE  APPLIED 

ON  THE  ANKLE 
<t>  =  LIMIT 

FORCE 

(N) 

122.73 

118.21 

77.48 

71.58 

150.75 

148.18 

MOMENT 

(N-m) 

81.24 

80.27 

58.69 

55.74 

115.65 

114.77 

154 


TABLE  XXIV 


MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  AND  MOMENT  VECTORS 
AT  THE  HIP  JOINTS  OF  THREE  FEMALE  SUBJECTS  FOR  THE  INDICATED 
POSITIONS  OF  THE  LEG  WITH  6  =  90°,  KNEE  FLEXED. 


LEG  POSITION 

MAGNITUDE 

SUBJECT 

FS  1 

SUBJECT 

FS  2 

SUBJECT 

FS  3 

FORCE  APPLIED 

ON  THE  KNEE 
$  a  0 

FORCE 

(N) 

152.27 

148.94 

152.48 

137.03 

89.66 

85.70 

MOMENT 

(N-m) 

65.22 

63.35 

74.04 

64.12 

38.21 

32.84 

FORCE  APPLIED 

ON  THE  ANKLE 
*  0 

FORCE 

(N) 

85.31 

75.21 

103.33 

94.27 

74.29 

60.87 

MOMENT 

(N-m) 

46.72 

38.77 

67.50 

61.47 

43.99 

36.20 

FORCE  APPLIED 

ON  THE  KNEE 

4>  =  30° 

FORCE 

(N) 

162.03 

159.67 

120.55 

120.37 

157.50 

152.97 

MOMENT 

(N-m) 

87.01 

84.58 

56.02 

55.50 

73.13 

70.08 

FORCE  APPLIED 

ON  THE  ANKLE 
<4  =  30° 

FORCE 

(N) 

104.54 

103.65 

81.15 

78.83 

132.52 

122.08 

MOMENT 

(N-m) 

67.25 

67.41 

64.62 

62.59 

78.08 

68.17 

FORCE  APPLIED 

ON  THE  KNEE 
<p  =  60° 

FORCE 

(N) 

248.47 

242.56 

160.26 

156.69 

209.12 

196.52 

MOMENT 

(N-m) 

109.61 

105.34 

75.57 

73.06 

95.97 

90.52 

FORCE  APPLIED 

ON  THE  ANKLE 
<j>  =  60° 

FORCE 

(N) 

132.46 

128.06 

86.49 

84.67 

150.92 

144.33 

MOMENT 

(N-m) 

79.41 

76.88 

62.12 

60.97 

99.27 

94.45 

FORCE  APPLIED 

ON  THE  KNEE 

4)  =  LIMIT 

FORCE 

(N) 

185.25 

183.18 

150.59 

149.47 

256.72 

234.51 

MOMENT 

(N-m) 

84.90 

83.31 

71.08 

71.14 

117.89 

112.69 

FORCE  APPLIED 

ON  THE  ANKLE 

4>  *  LIMIT 

FORCE 

(N) 

160.47 

149.07 

94.34 

78.24 

185.65 

174.91 

MOMENT 

(N-m) 

90.90 

83.66 

57.21 

43.79 

113. 4S 
103.65 

TABLE  XXV 


MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  AND  MOMENT  VECTORS 
AT  THE  HIP  JOINTS  OF  THREE  MALE  SUBJECTS  FOR  THE  INDICATED 
POSITIONS  OF  THE  LEG  WITH  6  =  120°,  KNEE  LOCKED. 


LEG  POSITION 

MAGNITUDE 

SUBJECT 

MS  1 

SUBJECT 

MS  2 

SUBJECT 

MS  3 

FORCE  APPLIED 

ON  THE  KNEE 
cp  =  0 

FORCE 

(N) 

95.27 

62.76 

119.22 

113.54 

90.25 

86.82 

MOMENT 

(N-m) 

43.24 

28.24 

56.50 

54.22 

37.54 

37.24 

FORCE  APPLIED 

ON  THE  ANKLE 

4  =  0 

FORCE 

(N) 

70.30 

66.54 

73.67 

70.27 

30.15 

29.86 

MOMENT 

(N-m) 

54.58 

51.71 

54.97 

51.58 

26.13 

25.86 

FORCE  APPLIED 

ON  THE  KNEE 

4  =  30° 

FORCE 

(N) 

211.95 

207.65 

256.31 

238.03 

145.24 

119.63 

MOMENT 

(N-m) 

97.90 

94.76 

122.13 

97.47 

65.94 

54.85 

FORCE  APPLIED 

ON  THE  ANKLE 
<p  =  30° 

FORCE 

(N) 

80.23 

75.92 

166.02 

151.33 

92.29 

90.18 

MOMENT 

(N-m) 

62.18 

58.54 

131.79 

109.71 

78.81 

76.91 

FORCE  APPLIED 

ON  THE  KNEE 

4  =  LIMIT 

FORCE 

(N) 

238.81 

231.80 

388.21 

343.50 

229.76 

219.12 

MOMENT 

(N-m) 

112.32 

109.51 

171.84 

148.21 

91.09 

87.01 

FORCE  APPLIED 

ON  THE  ANKLE 

4  =  LIMIT 

FORCE 

(N) 

157.92 

155.16 

177.56 

175.97 

101.06 

92.83 

MOMENT 

(N-m) 

124.91 

122.65 

138.69 

137.08 

75.50 

69.53 

156 


TABLE  XXVI 


MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  AND  MOMENT  VECTORS 
AT  THE  HIP  JOINTS  OF  THREE  MALE  SUBJECTS  FOR  THE  INDICATED 
POSITIONS  OF  THE  LEG  WITH  6  =  120°,  KNEE  FLEXED. 


LEG  POSITION 

MAGNITUDE 

SUBJECTS 

MS  1 

SUBJECTS 

MS  2 

SUBJECTS 

MS  3 

FORCE  APPLIED 

ON  THE  KNEE 

<j>  =  0 

FORCE 

(N) 

125.05 

118.38 

111.30 

100.73 

88.23 

86.67 

MOMENT 

(N-m) 

59.72 

57.10 

51.45 

45.72 

39.14 

38.30 

FORCE  APPLIED 

ON  THE  ANKLE 
$  =  0 

FORCE 

(N) 

79.23 

75.34 

95.84 

77.50 

63.49 

60.42 

MOMENT 

(N-m) 

47.00 

45.89 

51.11 

47.13 

34.88 

33.94 

FORCE  APPLIED 

ON  THE  KNEE 
i>  =  30° 

FORCE 

(N) 

168  79 
164.55 

214.14 

211.56 

153.79 

149.71 

MOMENT 

(N-m) 

77.20 

77.38 

92.18 

89.05 

61.31 

60.42 

FORCE  APPLIED 

ON  THE  ANKLE 
<J>  =  30° 

FORCE 

(N) 

106.08 

98.62 

127.08 

118.62 

226.19 

218.82 

MOMENT 

(N-m) 

66.65 

61.97 

55 . 62 
49.98 

97.45 

92.64 

FORCE  APPLIED 

ON  THE  KNEE 
<J>  =  60° 

FORCE 

(N) 

278.13 

275.60 

363.66 

348.63 

122.97 

112.15 

MOMENT 

(N-m) 

127.30 

126.00 

159.34 

156.51 

74.84 

66.94 

FORCE  APPLIED 

ON  THE  ANKLE 
<{>  =  60° 

FORCE 

(N) 

158.36 

151.98 

129.91 

123.21 

76.15 

69.57 

MOMENT 

(N-m) 

97.88 

93.00 

87.30 

80.78 

52.87 

47.41 

FORCE  APPLIED 

ON  THE  KNEE 
<p  =  LIMIT 

FORCE 

(N) 

282.46 

229.10 

373.14 

333.97 

246.03 

241.98 

MOMENT 

(N-m) 

128.14 

104.48 

170.30 

151.44 

129.97 

121.63 

FORCE  APPLIED 

ON  THE  ANKLE 
$  =  LIMIT 

FORCE 

(N) 

183.64 

158.55 

199.53 

186.21 

99.50 

95.45 

MOMENT 

(N-m) 

118.56 

101.99 

127.55 

118.86 

70.00 

67.86 

TABLE  XXVII 


MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  AND  MOMENT  VECTORS 
AT  THE  HIP  JOINTS  OF  THREE  FEMALE  SUBJECTS  FOR  THE  INDICATED 


POSITIONS 

OF  THE  LEG  WITH 

e  =  120°, 

KNEE  LOCKED. 

LEG  POSITION 

MAGNITUDE 

SUBJECT 

FS  I 

SUBJECT 

FS  2 

SUBJECT 

FS  3 

FORCE  APPLIED 

ON  THE  KNEE 

<)>  =  0 

FORCE 

(N) 

90.29 

83.41 

124.89 

107.25 

102.98 

94.30 

MOMENT 

(N-m) 

38.28 

37.58 

51.14 

43.39 

43.49 

38 . 30 

FORCE  APPLIED 

ON  THE  ANKLE 

<J>  =  0 

FORCE 

(N) 

38.74 

32.15 

31.38 

29.53 

48.66 

41.14 

MOMENT 

(N-m) 

29.46 

24.33 

25.28 

22.92 

35.55 

30.03 

FORCE  APPLIED 

ON  THE  KNEE 

4  =  30° 

FORCE 

(N) 

158.86 

152.75 

175.54 

166.62 

176.50 

158.17 

MOMENT 

(N-m) 

69.85 

64.97 

75.97 

71.84 

63.98 

56.21 

FORCE  APPLIED 

ON  THE  ANKLE 
<p  =  30° 

FORCE 

(N) 

67.08 

63.18 

79.88 

77.37 

81.36 

76.66 

MOMENT 

(N-m) 

50.91 

47.82 

62.75 

61.18 

64.57 

58.85 

FORCE  APPLIED 

ON  THE  KNEE 
<p  =  60° 

FORCE 

(N) 

194.23 

183.88 

164.48 

159.48 

300.22 

290.35 

MOMENT 

(N-m) 

65.68 

63.92 

79.70 

78.51 

132.96 

150.13 

FORCE  APPLIED 

ON  THE  ANKLE 
<t>  =  60° 

FORCE 

(N) 

107.63 

106.79 

85.54 

84.10 

216.26 

198.01 

MOMENT 

(N-m) 

67.16 

65.62 

69.20 

67.79 

131.19 

123.01 

FORCE  APPLIED 

ON  THE  KNEE 
$  =  LIMIT 

FORCE 

(N) 

237.61 

220.53 

192.01 

191.53 

528.22 

318.57 

MOMENT 

(N-m) 

80.38 

74.56 

90.85 

88.44 

109.14 

111.03 

FORCE  APPLIED 

ON  THE  ANKLE 
<p  *  LIMIT 

FORCE 

(N) 

136.84 

134.97 

94.45 

91.99 

196.53 

185.14 

MOMENT 

(N-m) 

107.59 

105.50 

78.67 

76.92 

121.08 

120.42 

158 


TABLE  XXVIII 


MAXIMUM  MAGNITUDES  OF  THE  ACTIVE  MUSCLE  FORCE  AND  MOMENT  VECTORS 
AT  THE  HIP  JOINTS  OF  THREE  FEMALE  SUBJECTS  FOR  THE  INDICATED 
POSITIONS  OF  THE  LEG  WITH  0  =  120°,  KNEE  FLEXED. 


LEG  POSITION 

MAGNITUDE 

SUBJECT 

FS  1 

SUBJECT 

FS  2 

SUBJECT 

FS  3 

FORCE  APPLIED 

ON  THE  KNEE 

4>  =  0 

FORCE 

(N) 

73.04 

63.44 

128.22 

121.21 

118.36 

105.76 

MOMENT 

(N-m) 

26.83 

23.82 

48.71 

47.27 

44.83 

38.11 

FORCE  APPLIED 

ON  THE  ANKLE 

<j>  =  0 

FORCE 

(N) 

65.91 

64.64 

85.04 

76.72 

67.72 

63.20 

MOMENT 

(N-m) 

43.00 

42.14 

26.13 

23.62 

22.29 

21.49 

FORCE  APPLIED 

ON  THE  KNEE 
<P  =  30° 

FORCE 

(N) 

148.41 

142.21 

168.18 

163.47 

166.79 

164.85 

MOMENT 

(N-m) 

62.48 

59.89 

81.61 

79.10 

53.31 

51.46 

FORCE  APPLIED 

ON  THE  ANKLE 
<f>  =  30° 

FORCE 

(N) 

86.63 

80.48 

92.87 

87.28 

159.17 

139.16 

MOMENT 

(N-m) 

29.76 

28.72 

55.34 

52.83 

88.94 

79.16 

FORCE  APPLIED 

FORCE 

(N) 

208.35 

195.93 

161.91 

155.49 

219.02 

215.78 

UN  lllu  NNLL 

4  =  60° 

MOMENT 

(N-m) 

92.03 

86.96 

68.73 

65.66 

64.00 

64.38 

FORCE  APPLIED 

ON  THE  ANKLE 

4  =  60° 

FORCE 

(N) 

133.83 

122.38 

85.62 

84.89 

209.98 

170.52 

MOMENT 

(N-m) 

74.55 

72.80 

52.75 

51.65 

84.09 

62.67 

FORCE  APPLIED 

f\M  TUC  L'MEE 

FORCE 

(N) 

237.32 

218.11 

194.45 

187.26 

250.91 

249.38 

UN  IMl  Mill 

4  =  LIMIT 

MOMENT 

(N-m) 

107.06 

98.91 

86.77 

85.28 

111.87 

114.84 

FORCE  APPLIED 

ON  THE  ANKLE 

4  -  LIMIT 

FORCE 

(N) 

146.90 

146.00 

81.03 

75.25 

216.37 

196.18 

MOMENT 

(N-m) 

99.37 

99.00 

41.14 

40.78 

127.11 

114.03 
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